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The present paper treats the symmetrical and charged pseudoscalar theories of the meson 
field, using the strong coupling approximation; it restricts itself to the case of a single source. 
The energy levels of the excited states of the heavy particle and the scattering cross section 
for free mesons are computed by wave mechanical methods. An expression is also obtained 
for the magnetic moment of the proton or neutron. While the scattering cross section can, 
with reasonable assumptions, be brought into agreement with experimental values, the results 
for the magnetic moment are qualitatively at variance with the known values in that equal 
and opposite moments are predicted for proton and neutron. 


1. INTRODUCTION 


HE perturbation theoretic treatment of the 

coupling of the meson field to the heavy 
particle, based on the assumption that the 
coupling is weak, encounters several funda- 
mental difficulties in its application. The di- 
vergences which arise from the treatment of the 
heavy particle as an infinite point mass have 
been in some instances arbitrarily removed by 
cut-off methods, but it can be shown that such 
methods cannot be consistently formulated in 
such a way as to satisfy the criterion for weak 
coupling.' Moreover, calculations of the scatter- 
ing cross section for mesons in this theory lead 
to values which are generally much too high to 
agree with cross sections observed for cosmic-ray 
mesons. The results vary somewhat depending 


1 This point will be more fully discussed in a subsequent 

per in which the two-source problem will be considered. 

e are indebted to Professor R. Serber for valuable 
discussions on this point as well as on several others 
considered in this paper. 


on the particular type of meson field assumed. 
If we assume that the spin of the meson is zero* 
(charged or symmetrical pseudoscalar), pertur- 
bation theory leads to the following result* for 
the scattering of a meson by a nucleon.‘ 


dq=gix*(p/x)*(1+cos for hp>Mc, 
for hp<Mc, 


Here g is the dimensionless constant which 
expresses the magnitude of the coupling ; «= yuc/h 
where yu is the rest mass of the meson; p is the 
momentum of the incident meson divided by h; 
E the total energy divided by hc; M the rest 
mass of the nucleon; and @ the scattering angle. 
These cross sections refer to processes in which 


2 This is one of the two possibilities left open by a | 
consideration of the electromagnetic radiation processes 
of the meson. See R. F. Christy and S. Kusaka, Phys. 
Rev. 59, 405 and 414 (1941). The other possible value 
of the spin (4) will be discussed by J. R. Oppenheimer 
and E. Nelson in a forthcoming publication. 

3H. Yukawa and Y. Tanikawa, Proc. Phys. Math. Soc. 
Japan 23, 445 (1941). 

4 Nucleon is equivalent to ‘‘proton-neutron.” 
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the meson is scattered with the original charge. 
In the case of the symmetrical theory there 
exists also a scattering process in which the 
meson changes its charge. If a negative (positive) 
meson collides with a proton (neutron) it can 
be scattered as a neutral meson, while the 
nucleon changes into a neutron (proton). Ac- 
cording to computations of F. Adler, not yet 
published, the corresponding cross section is 


dq =g'x~*(p/x)?(1—cos 6/2)*d2 for hp>Mc, 
dq = sin? 6dQ for hp<Mc, ~~. 


A numerical estimate of the total cross section 
depends upon an assignment of the magnitude 
of g. This is ordinarily done by considering the 
nuclear forces predicted by the theory and 
adjusting g to fit the properties of the deuteron. 
This involves a cut-off or other readjustment of 
the radial dependence of the interaction. As 
remarked above, such a procedure cannot be 
carried through in a manner consistent with the 
weak coupling hypothesis. Nevertheless, the 
value of g so determined (g?~0.1) will be used 
to give an order of magnitude estimate of the 
scattering cross section. We obtain a total cross 
section ~2X10-** cm?, as compared with the 
upper limit of 510-*% cm? determined from 
experiments on cosmic-ray mesons.$ 

Two apparently different theories have been 
proposed to explain the smallness of the meson 
scattering cross section. Heisenberg pointed out 
that the reaction of the eigenfield of the nucleon 
to the motion of its spin—and particularly the 
terms proportional to 1/a which one can con- 
sider as an inertia of the spin—are considerable.® 
As an example, he computed by classical 
methods, the scattering in a neutral theory with 
spin 1 mesons; he found a scattering cross 
section proportional to a*(p/£)! if 

The second theory is due to Heitler® and 
Bhabha® and based 6n the assumption of the 


5 For instance, see R. P. Shutt, Phys. Rev. 61, 6 (1942). 
Older literature is reviewed. 

‘The quantity @ measures the radius of the nucleon. 
A precise definition is given by (4). 

7W. Heisenberg, Zeits. f. Physik 113, 61 (1939). The 
problem of interaction of a dipole with its own electro- 
magnetic or meson field is treated in detail by H. J. 
Bhabha and H. C. Corben, Proc. Roy. Soc. A178, 273 
(1941) and H. J. Bhabha, Proc. Roy. Soc. A178, 314 
(1941). A classical but relativistically invariant method 
is used. The terms in the Hamiltonian proportional to 
1/a are replaced by terms containing undetermined 
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existence of excited states (isobars) of the nucleon 
with higher values of the charge (Bhabha) and 
of the spin (Heitler). In this theory the energy 
difference AE between consecutive isobaric states 
is not derived, but arbitrarily assumed. On the 
basis of these assumptions, one finds a con- 
siderably reduced value for the cross section for 
scattering of a meson by a nucleon (of infinite 
mass) at rest: 


= q0(AE/E)’. 


E is the meson energy and gp is the total cross 
section (given above) for scattering as calculated 
by perturbation theory. 

It was shown by Oppenheimer and Schwinger’? 
that the two explanations are in reality one and 
the same. One can indeed derive from Heisen- 
berg’s assumptions the existence of excited 
states of the nucleon, corresponding to the degree 
of freedom associated with the reaction of the 
nucleon spin. It is easy to see that if we assume 
this excitation energy proportional to 


AE~ 
and combine this with Heisenberg’s result 
qi~a*(p/E)* 


we obtain Heitler and Bhabha’s formula for 
qi/qo. The above formula for AE, presumably 
valid for vector mesons, is identical with the 
corresponding formula that will be rigorously 
derived for the pseudoscalar case (see below). 


2. RESULTS AND CONCLUSIONS 


We give here the principal results of this 
investigation, including the excitation of the 
isobaric states, the cross section for meson 
scattering, and the resultant magnetic moments 
of the proton and neutron. As the condition for 
the validity of the strong coupling approximation 
we obtain g>xa, in the case of small source size 


constants. Such a procedure may approximate an exact 
treatment more closely than the one used in this paper, 
but it is not suitable for a quantum-mechanical treatment 
of the problem. 

8’ Compare the detailed paper of W. Heitler and S. T 
Ma, Proc. Roy. Soc. A176, 368 (1940). 

9H. J. Bhabha, Phys. Rev. 59, 100 (1941). We refer to 
this paper for the discussion of the influence of the Coulomb 
forces on the probability of the generation of a doubly 
charged proton. 

10]. R. Oppenheimer and J. Schwinger, Phys. Rev. 60, 
150 (1941). 
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(xa<1). For the energy of excitation of the 
isobars we find , 


AE 3 3]. 

j(j+1) in the symmetrical theory, 
uc? 2g? 4 

AE 5]. 

— 2j(j+1) —n*?—-— |in the charged theory. 
uc? 2¢g 4 


Here the half-integer j is the total angular 
momentum, and the half-integer quantum 
number 1 is the electric charge (in units e) minus 
one-half. In both cases, the value of 7 is restricted 
to —j=n=j, otherwise the isobar energy in the 
first case is independent of n. The definition of 
the radius a of the nucleon is given by (4). The 
proton and neutron are identified as the two 
lowest states of this system, corresponding to 
j=}, n= +}. Higher states will be stable against 
“meson decay’’ to these ground states only if 
AE<uc? or g?/xa>1. 

For the total scattering cross section (elastic 
plus inelastic scattering with and without change 
of the charge) we find under the conditions 
g>«ka, ea<1, pa<1: 


dq = }(p/E)*a*(1+cos* 6)dQ, 
q=4n(p/E)*‘a? in the symmetrical theory; 
dq = 4(p/E)‘a*(1+3 cos? @)dQ, 
q=6n(p/E)*a? in the charged theory. 


As has already been shown by Oppenheimer and 
Schwinger, one gets agreement with the experi- 
mental cross section for meson scattering, in the 
symmetrical theory for instance, by assuming 
a~h/Mc=2X10-" cm which means xa~0.1. 
Referring to the condition for strong coupling, 
g>>xa, we see that values of g?>0.01 satisfy this 
condition. 

The magnetic moment of the system nucleon 
+meson field is found to be 

10g? 1 


+]— proton magnetons 
E ka 6 


the + symbol referring to proton and neutron, 
respectively. As regards order of magnitude; this 
is not inconsistent with the known magnitudes. 
A moment of 1.93 proton magnetons, such as 
has been observed for the neutron, can be 
obtained by taking g?/xa=6.3; using xa=0.1 
we find g?=0.63, well into the strong coupling 


domain. However, the prediction that the 
neutron and proton moments are equal and 
opposite is in contradiction with experience. 
The second term in the bracket above represents 
the contribution of the nucleon itself. This 
contribution is somewhat uncertain because of 
our ignorance of the fundamental properties of 
the ‘‘bare’’ proton or neutron. Consequently 
the discrepancy cannot be taken as an unam- 
biguous disproof of the strong coupling hy- 
pothesis. However, the hope that the anomalous 
moments would follow in a simple way from the 
theory has not been justified. 

It is to be noted that the results obtained 
above are quantitatively different from the 
results of the application of semi-classical meth- 
ods employed by Oppenheimer and Schwinger 
and described in detail by Schwinger." The 
latter method treats the spin @ and the isotopic 
spin + as classical unit vectors. In the charged 
scalar and neutral pseudoscalar theories, it gave 
a quantitative agreement with the wave me- 
chanical results in the limit of strong coupling. 
This is, however, not true in the present case; 
in particular the restriction —j=n3j does not 
appear. The discrepancy cannot be considered 
surprising in view of the neglect, in the semi- 
classical method, of the commutation relations 
between the various components of + and @. 
One sees in the wave mechanical treatment that 
these relations play an essential part in deter- 
mining the minimum eigenvalue of the inter- 
action energy (see Section 4). Once this part of 
the problem has been solved, however, there is 
no difference between the two methods of 
treatment. 

The remaining sections are devoted to the 
derivation of the results quoted above. In 
Sections 3 to 6 we introduce such new field 
variables for the symmetrical pseudoscalar field 
as to express the energy in terms of unbound 
mesons and mesons bound to the nucleon. In 6 
we find the dependence of the energy on the 
charge and angular momentum of the bound 
meson cloud, and discuss in detail the conditions 
for validity of the various strong coupling 
approximations introduced. In 7 we calculate 
the scattering cross sections. Section 8 contains a 


"J. Schwinger, publication in preparation. 


= 


specialization of all these calculations as applied 
to the charged pseudoscalar case. Section 9 
consists of the magnetic moment computation 
and an appendix describes an alternate method 
for expanding field quantities, with the use of 
Euler angles. 


3. SPLIT OF THE FIELD INTO A ZERO STATE 
AND UPPER STATES 


In the symmetrical pseudoscalar theory it is 
convenient to describe the field by three real 
quantities g_(x), a=1, 2, 3. The Hamiltonian H 
consists of the part H of free particles and the 
interaction energy H,. The first part is given by” 


Hy= (Vga)? +k pa? |d V 
a=1 


3 
= [re ¢a(—A+k) ga (1) 
a=1 
where the 7, are the momenta conjugate to ga 
and fulfill, at a given instant of time, the com- 
mutation rules 


iL wa(X), ]= 5ag5(x—x’). (2) 
In the interaction energy we suppose the heavy 
particle, which we shall call nucleon, to be at 
rest and characterized by a spherically sym- 
metrical source function K(x) =K(|x|) which is 
normalized according to 


V=1 (3) 


and determines a radius a of the nucleon: 
ote f f (4) 


where here and in the following, r has the meaning 
(5) 


Furthermore, the nucleon is capable of existing 
in two states corresponding to proton and 
neutron, and moreover to two states correspond- 
ing to different directions of the spin. The first 
degree of freedom is described by the isotopic 
spin matrices 71, T2, 73, OF Ta With a=1, 2, 3; the 


r=|x—x’|. 


122 We always use here natural units where h=c=1 or, 
in other words, energies are divided by hc, momenta by h, 
angular momenta by &, electric charges by (hc)*; x is the 
rest mass u of the free meson divided by he. 
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second by the spin matrices ¢,, ¢,, ¢:, or «. Both 
of them fulfill the same kind of relations 


(7) 
where the --- denotes similar relations derived 


by cyclic permutations. The interaction energy 
is then given by 


= (g/x)(4r)!/v2 [VK (8) 


The factor (47)! is introduced in order to measure 
g in ordinary units, not in Heaviside units; the 
factor 1/v2 in order to bring it into accordance 
with the notation in the theories which introduce 
charged particles only; and the factor 1/x in 
order to make g dimensionless." The total charge 
of the system meson field+nucleon, measured in 
the unit e of the electron charge is given by 

c= (9) 
The symmetrical theories are characterized by 
the fact that the charge is only one component 
of the isotopic spin Tas= —T ga (a, B=1, 2, 3) 
which is a more general integral of motion and 
given by 


with — 721=73, The charge is then the 
12 component of the isotopic spin, e= 712+}. 

The form (8) of the interaction energy suggests 
the definition" 


— (4x)! f K(x) V 
OX; 
=f alV, (11) 
’ 


which gives 


K a,k 


Moreover, we split the field g(x) into two parts 


18 The connection of our constant g with the constant 
gy introduced by Yukawa is 
g=fy(hc) 4. 
4 Small roman indices 7, k, --- run from 1 to 3 and 
denote vector components in the ordinary space. 
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of which one, ¢a’(x), is orthogonal to the gradient 
of the source function 


f ¢a'VKdV =0 (12) 


and the other proportional to the gradient of the 
potential X(x) which the source function gener- 
ates according to 


(—A+x?)X=4rK, 
= f (13) 


If we define 
OX 0K 
= f —dV 
OX; OX; 
0K dK 
= —w—dVdV’, (14) 
J OX; F OX, 
X(x) OK 
&(x) =——, bu= (15) 
I Ox; OX, 


then it follows from (11) and (12) that 


0 
hence 


= ga(x) —VE- f VKeadV. (17) 


The circumstance that we choose the first part of 
the field as a linear function of the 0X/dx, (and 
not of the 0K /dx;) has the consequence that in 


virtue of (12) no cross terms between eas and 
¢a (x) occur in the potential energy //». Indeed, 
one has, using (13), (15), 
0 
(Gar)? 


BE 
+32 f ga'(—A+k*) V. (18) 


The conjugate decomposition of the momentum 
m(x) is given by 


aK 

k OX, 

f x.'VidV=0, (12') 

wen f (at/dx,)medV, 

we! = f VirdV, —(17’) 


which leads, in accordance with (2), to the 


commutation rules 
iL mes, ]= bap dik, (19) 
(20) 


iL ma'(X1), op’ (X2) ] 
= bas(6(X1—X2) —ViK-Vok). (21) 


The latter relation is also in direct agreement 
with the orthogonality conditions (12), (12). 
The isotopic spin (10) decomposes simply into 


f ma’ — ta’ )dV 


0 
The kinetic energy however gets cross terms in 


0 
Tai and m,_'(x), and if we use the abbreviation 


0K dK 
——dV 


Ox; 
or 


f (VK)%4V, (23) 


it is given by 


f roid V=3N (mas) 
us ak 
+4) Ex j= (x)d VV 


445 f (24) 


One gets the total Hamiltonian by collecting (24), 
(18), (8’): 
0 0 
H=3N (mar)? +3 
a,k a,k 


a,k 


+45 f 
+45 f (I) 


Just as for the isotopic spin there exists in virtue 
of the spherical symmetry of the source function 
K(x) the integral of the angular momentum 
—Li; given by 


xf ( : eal V+4 (10’) 


)) 
e 
n 
e 

it 
) 
) 
e 
| 
) 
rt 


It decomposes as follows: 
tout La ( - 


OX, Ox) 


fr ¢ga'dV. (22') 
Ox; 
The orthogonality relations (12), (12’) do not 
make the cross terms disappear in this case. We 
emphasize, however, the simple form of the 
angular momentum of the zero field, given by 
the first sum in (22’). 


4. THE EIGENVALUES OF THE INTERACTION 
ENERGY. NEW VARIABLES FOR 
THE ZERO STATE 

While in the usual treatment the interaction 
energy H, is considered as a small perturbation 
(weak coupling), we consider in this paper the 
opposite limit of strong coupling. In this case one 
has first to investigate the eigenvalues of the 
interaction energy H, given by (8’), and then to 
retain only the lowest of them under the assump- 
tion, whose validity will be investigated, that the 
next higher level of 77; will not appreciably 
perturb a system in this lowest state. 

In order to determine the four eigenvalues of 
the form 

a,k 

in which the 9 coefficients a were real, we have 
to bear in mind that we are free to subject the 
7's and o’s to independent real orthogonal 
transformations 


3 3 
Taz >, tr Bra, (25) 
r=1 s=1 
where B and A fulfill orthogonality conditions 
which we can write in matrix notation 


AA=AA=1, BB=BB=1, (26) 


where A and B are the transposed matrices 
(An=Arse, Bar=Bra). The eigenvalues of our 
form are therefore the same as those of the other 


form 
(27) 


r,s 
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We shall prove now that we can choose the or- 
thogonal matrices A, B in such a way that the 
new matrix Byg®A becomes diagonal: 


(B deo (28) 


or, according to (26), 


(29) 


For this purpose we consider the two sym- 
metrical forms 


and 


, 0 0 
C’=¢0 or Cap=Cpa= Li Garver. 
We choose A as the orthogonal matrix which 
transforms C to principal axes (which is always 
possible because C is symmetric) and denote the 


eigenvalues of C by (Q,)*. These are certainly 
never negative because the form 


Xe = Le (dx Parte)? 
i,k 


is definite. This means that we have 
CA=AQ? or ACA=(Q?. (31) 


We note that the matrix A is unique except for a 
sign, and except for the case where some of the 
Q,’s are zero or some of them are equal to each 
other. We now define the matrix B by 


vA = BQ, (32) 


where Q, is defined as the square root of the 
eigenvalue Q,? of C, the sign of which we shall 
choose later. Multiplication by ¢ from the left 
gives, by virtue of (30), (31) 


AQ= 
AQ=eB, AAQ=Q=ABB, 
and by transposing the matrices, Q being diagonal 
BeA=Q 


which proves Eq. (28). Now multiplying (32) by 
B from the left, we get 


BeA =Q=BBO. 


Therefore BB =1, which proves the orthogonality 
of B. (For the case of a multiple eigenvalue Q, 
the orthogonality BB=1 does not follow, but 
can easily be achieved.) Multiplying AQ= 


or 


A 


b 


al 


- 


it 
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by Bg from the left, we get 
BC'=@B or BC'B=@, 


which is analogous to (31) and shows that C’ has 
the same eigenvalues as C. 

If the eigenvalue Q, is zero, the row B,. of B 
is not defined by (32) and the preceding proof is 
incomplete. Let us assume first that there is a 
simple eigenvalue Q;=0, Q:+0, Q2#0. From 
Eq. (32) we get for r=3 


The equation AQ=@B is in this case not a con- 
sequence of the proved equation AQ?= gBQ but 
can be used for r=3 as a definition of the third 
row of the matrix B, namely 


= 0. 
This can always be fulfilled, the determinant of 


(33) 


da being zero in this case. There is still an arbi- 
trary factor independent of a@ undetermined in 
Bz. This factor can be adjusted to satisfy the 
equation Q=BBQ; it can simultaneously be 
made to satisfy the equation 


which does not follow directly from the above. 


Combining Q=BBQ and ¥.Bse=1, we obtain 
again BB=1. 
In Section 8 we shall meet the particular case 


Q:=0, Q:¥0, Q2¥0, 


In this case one has B3.=0 for a=1, 2, which 
follows from 


3a >= 0. 


Moreover, (32) leads to Bas=0 for a=1, 2. 
Finally, B33 is fixed by BB=1 and is equal to 
one. The latter condition also expresses the 
orthogonality of the two-row and two-column 
matrix B,. (r=1, 2; a=1, 2). 

In a similar way the case of a multiple eigen- 
value zero of C and C’ can be treated by con- 
sidering all equations in the preceding proof 
which cannot be directly derived as possible 
additional defining postulates. 

In order to have the determinants of A and B 
both +1, which is essential for the validity of the 


correct algebraic relations for the r,’, o,’ we get 
from (32) the condition 


det. ¢=Q10203. (34) 
Apart from this condition, the signs of the Q,’s 
are arbitrary and can be fixed by definition. This 
means that to every combination of signs of the 
Q, which fulfills (34) there exist, for a given 


Pe matrices A and B which are orthogonal in 
the usual sense and which fulfill (28). We can for 
instance define all signs of the Q,’s as positive 
(negative) if det. ¢ is positive (negative). 

Now inserting (28) in (27) we get 


a,k r 


The eigenvalues of every 7,’, ,’ are +1, but the 
product of the three matrices 7,’¢,’ is always —1. 
Therefore, the eigenvalues of 17, for a given field 


Gar are 
(g/xv2)>, with (36) 
E,= —(g/«v2)(Qit+Q2+Qs), 
E2= (g/«v2)(—Qi+Q2+Qs), 
E3=(g/xv2)(Qi—Q2+Qs), 


The sum of all four E’s is zero, but if E is an 
eigenvalue in general —£ is not an eigenvalue 
unless at least one of the three Q,’s is zero. We 
repeat that the Q,’s are the square roots of the 
eigenvalues of the symmetrical matrix C (or C’). 
From (31) we get immediately by evaluating the 
trace 


(36’) 


(var)? (37) 


r=1 k 
The sum of the Q,’s themselves, however, cannot 
be expressed rationally by the field i 


We now introduce into the kinetic energy of 
the zero field the variables corresponding to Q,, 
Ax, Bra. Instead of expressing the latter variables 
by Euler angles and using their conjugate 
momenta, we can also use the components of the 
angular momentum 


and the isotopic spin 
0 0 0 ’ 
Tap = Lit (Par Tak). (220) 


1e 
*h 
ys 
1e 
| 
1e 
+h | 
i 
1e 
ft 
_| 
), 
_| 
B 


As operators they correspond to infinitesimal ro- 
tations of the ordinary space and of the isotopic 
spin space, and fulfill commutation relations with 
the A;,, and B,, analogous to those with the 


namely 
A = ir — Bit 

(Lis, (38) 

Aur]=0, (38’) 


0 0 
while the components of Ly, or Tag fulfill with 
each other the well-known commutation relations 


Lim] = 
(39) 
or, with the vector notation L,°, L2°, L;° for 
Les, Ln, Liz 
iL L1°, L2° |= —L;°, ---. (39) 
As scalars the Q,’s commute with both the 
Liz and the 
[Lit, Q-]=[Tas, Q,]=0. (40) 
In the same way we have 
i[ Tos, = as 
+ — Say (39) 
the and the commute. 
It is useful to introduce the components of 


(L°), (T°) with respect to the axes defined by 
(Arr), (Bra) which we define by 


i,k 


i,k 


(41) 
a,B a,B 


We note that in spite of the non-commutativity 


of Ly with the A;, and of Ts with the B,, the 
Lo” and 7, are Hermitian operators like the 

0 0 
Ly and Tag as a consequence of the equalities 
indicated in (41) which follow from (38), (38’). 
The total square of Lo and 7) is given by 


i<k 
(Tos)*=E (To")?. (42’) 


The commutation relations of the Ly” with the 
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A,, and of To" with the B,. are given by 
Ane] cs — kr, (43) 
iLT0", Bia 5stBra, (43’) 
and of the Lo” or the To" with each other 
Lo", Lo] =6,.L0™ 


iLTo", To ]=6.T0™ 
+ To" bru To". (44’) 
The signs in the relations (43), (44) or (43), 
(44’) are just opposite to those in the relations 
(38), (39) or (38’), (39’’). 
Inserting (29), (220), and (22)’) into (41) we get 


(45) 


ka 
From this we get the important expression 
1 To" To" 


Finally, we get the variable canonically conju- 
gate to Q,, which we call P,, by putting for r=s, 


= P,. (47) 


for (46) 


Therefore, we have finally 
Tes > AirBsa 


1 Lo" T 1 Lot* — To" 
x 0 + 0 0 0 
where the second term is to be taken for 7 #s only. 
Assuming 


Ax. |=(P,, B..)= iP. Lo**] 
one can verify the canonical commutation rela- 


tions (19) for the m2; and gs with the help of (43) 
(44), (43’), (44’). It may be noted that the 
operator P, defined by (47) is not Hermitian and 
that P, in (48) has to be replaced by P,* if the 
factor Aj,B.a is put on the right side of the 
bracket. For the kinetic energy of the zero field 
(apart from the factor N/2) we get from (48), 
taking into account the anti-symmetry in 7, s of 


» (48) 


8) 
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(Lo — To"*)/Q,+Q,; and the symmetry in r,s of 
(Lo”* + To") /Q; Q; 


0 +7 
1 — To")? 
(© 


r,8 (Q,+Q.)? 


We conclude from (46) that if @,=Q, then it 
must also be true that Lo"*= — 79". In view of the 
importance of this case it is sometimes convenient 
to use different coordinates. Indeed, if all three 
Q,’s are equal, which will be proved to be true in 
good approximation in the case of large coupling, 
according to (29) only the orthogonal matrix 


Cak = AirBra (51) 
which fulfills 
Dia Cailak=Sixy Dox (51a) 


is uniquely determined by the field gx and not 
the A,, and B,. themselves. Therefore, besides 
the three independent variables e,,, we can intro- 
duce the six symmetrical quantities 


0 0 
das = Pakepk = B,aB,sQ, = Gea (52) 
to describe the field 
= Ds (53) 


The components 


(54) 


r,8 i,k 


0 
commute with the qg,s and fulfill similar commu- 
tation relations with each other and with the e,; 
as do the Lo” with each other and with the Ax, 


iL eyx_] = — 53y€aks (55) 
iL Lo, Lo” |= 

+ — — (56) 

Finally, we search for symmetrical quantities 


Pas= psa Which are conjugate to the gag in the 
sense that they fulfill the commutation relations 


iL pes, (57) 


(This means, for instance, that i[pi1°, qi1°]=1, 
iL p12", ]=4; this is more convenient than to 
have the value 1 for the latter bracket.) This is 
just the case if, with the help of the first two 


terms in the bracket in (48), we put 
1 To" 


~ B,aBrgP, +>’ (58) 
r¥s Q, 
From as) we find also 
(Cart n+ (59) 


The remark made about the non-Hermitian 
character of the P, in (47) applies also to the 


Pe defined here. One can prove (57) with the 
help of (54), (43), (43’), and moreover 


[pes, Lo ]=0. (60) 


We have therefore three kinds of variables for the 


zero field: 
0 0 


(1) Pak, 

(2) Q,, Aur, P., Lo", To", 
0 0 

(3) das Cak; Pab=Psa, Lo®®. 


While we were able to get an explicit expression 
for the kinetic energy in terms of the variables 
(2), this is not rigorously possible with the 
variables (3). From (48), ~ we get only 

ot 0 — To")? 


abPa 61 


In rel next section we shall see, however, that in 
the case of small deviations of the Q, from a large 
common value Q,°= D, further simplifications are 
possible, and that the variables (3) are the most 
convenient ones in this case. 


5. MINIMUM OF THE POTENTIAL ENERGY. 
STRONG COUPLING APPROXIMATION 


We determine first the minimum of the total 
potential energy of the zero field which is, ac- 
to (18), (37), equal to 


2] r=1 
(€1€2¢3= —1). (62) 
The minimum corresponds to 
g 
a=@=eg=—-1, 


3 
Enin= (63) 


According to (53) the corresponding field is given 
by 


= Dear (64) 


3) 
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or, in vector notation 


¢a° = Deg. (64a) 


This field corresponds to the smallest eigenvalue 
of the interaction energy, determined by (36), 
and to the absence of free mesons [we neglected 
but still contains the g numbers e,. 
Moreover, it corresponds to the state of the r 
and o where 

TaS = —3. 

a,k 
By a suitable S transformation of the Hamil- 
tonian (see Appendix) we can bring this to the 
form 


'= TaVa- (65) 


a=1 

The S transformation is not uniquely deter- 
mined, just as is true for the A and B matrices; 
it is, for instance, sufficient in this case to trans- 
form one system of the matrices 7, or o; alone. 
In view of the application to the charged pseudo- 
scalar theory, which is discussed in Section 6, we 
decided to leave the 7.’s unchanged and to 
transform the o;’s, 


St_57! =Ta; S(d: (65a) 
This is equivalent to 
on! => = CakTa: (65b) 


The state in question can be described as the 
singlet state inasmuch as the composition of the 
T and og is analogous to the composition of two 
spins of: equal magnitude. This state is not 
degenerate according to (36’), while the three 
other eigenvalues of the interaction energy are 
equal : 


E,= —(g/xv2)3D, Ex=E3=Es=(g/xv2)D. (66) 


In the following we assume that it is legitimate 
to restrict our attention to the lowest, or singlet 
state. This implies a condition on the magnitude 
of such terms in the Hamiltonian as couple the 
lowest to the three upper states of the inter- 
action energy. 

For the eigenvalue £;, the following expres- 
sions are diagonal: 


(67) 


15 It will be shown in Section 7 that the field gy’ describes 
free mesons. 


7101 = T202=7303= —1; 


On the other hand, the matrices 7,, ¢,, and 
7192— 7201, *** are not diagonal. We shall simply 
neglect terms in the energy in which these 
matrices appear as a factor. In Section 6, we 
shall discuss the region of validity of this ap- 
proximation. 

As is shown in the Appendix, the S trans- 
formation changes the angular momentum L 
according to 


Ya (68) 


As a consequence of this and (54) we have the 


relation 
SL 5% S-! = Lo? — (69) 


One checks directly that L’ commutes with o’ if 
L commutes with @ using (55). As the o. can be 
disregarded in our strong coupling approximation 
if they do not occur multiplied by the 7,’s, the 
L’ have in this approximation the same eigen- 
values as the L, while the total sum L” [see (42) ] 
has eigenvalues } larger than L*. Moreover, the 
eigenfunction W’=SwW is double-valued in the 
angles which enter in S and in which W is one- 
valued. Therefore (L’)? has the eigenvalues 
3+ j(j+1) with half-integer 7. 


We now split the variables Ges defined by (52) 
into two parts according to (64) 
(70) 


of which we consider the second part as small. 
With D given by (63) 


1 0 2 4 0 
This leads, by the S transformation (65) which 


0 
does not change the i but does change the max, 
to 


, 1 2 
Epot = (gaa) 
as 


= 


«v2 a,B 
Using (67) we have 
1 0.2 g 0 
pot 


which simply reduces to 


Sgt 
Kaw 


71 
4x2 71) 


aks 
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We now introduce (64) into (229), which means 
that we disregard os and get 
Lie (Caitak— 
and from (54) 


It therefore follows from (59) that 
1 
Tat = en( 


(72) 


Inserting es = De, one can finally check directly 
= (73) 


where the disregarded terms are of the relative 
order of magnitude D~ in comparison with the 
main term. This is in agreement with the result 
of the last section. If we introduce T»’*= —Lo"* 
and Q,+Q,=2D in the last term in (48) we get 
(72) directly. The relation (73) has the important 
consequence that the total square of the isotopic 
spin is equal to the total square of the angular 
momentum; after the S transformation, both 
correspond to half-integer quantum numbers. 
Therefore the quantum number 1, which is equal 
to the charge number minus } [see Eqs. (9), 
(10) ], always has for the zero field (that means 
in absence of free mesons) an absolute value 
smaller than the total angular momentum 
quantum number j: 


(74) 
This is of course due to the fact that the r, and 
oq are in the ‘‘singlet state’’ defined by (67) ; this 


also enables us to disregard the terms linear in 7 
and o as already mentioned. 


t 
If we introduce (pas)u = 3(Pas+ Pas), the ( 
denoting Hermitization, we get 


1 
=> +—(L*%eg;.) 4. (72’) 
8 2D 


This leads to the following expression for the 
kinetic energy of the zero field (neglecting an 
additive constant in the term proportional to 
D-): 
N o2 N o 2 1 
— at) = — —T,*}, (72a 


where the square of the total angular momentum 
is given by (42). A result identical with (72a) is 
obtained from (61) if we insert in the last term 
Q0,+Q0,=2D, —To'*=Lo". We note that the 
only result of the S transformation is that j gets 
half-integral values, the constant } added to Lo? 
being negligible for our applications. 

Inserting (71), (72), (72a) into the Hamiltonian 
(1), we get the ‘‘reduced” Hamiltonian of the 
strong approximation 


3g 


21 a,B a, B,k 


1 
x f > fra V 
OX), 2 a 
1 
fi 


1 0K 

+(4r)! f —ma'(x)dV. (IT) 
a, B,k 2D Ox, 

Inspection with respect to the approximation 

contained in (72) shows that the disregarded 

terms in H/ are of the following type: terms of 


the order D~* independent of Pal terms of the 


order D-? linear in the Li, and in the qu; and 
terms of higher order than D-*. We shall see in 
the next two sections that such terms are of no 
importance for the questions which are treated 
in this paper. The last term in (II), however, 
will be essential. 


6. THE EXCITED STATES (ISOBARS) OF 
THE NUCLEON 
We shall now treat two different kinds of 
problems. In the first kind we investigate the 
states of the nucleon in the absence of free 
mesons, while in the second kind we treat the 
scattering of free mesons by the nucleon and its 
proper field. It will turn out in the next section 


that the field which is described by Pas, Ta (x); 


Pi ¢a'(x) represents free mesons, which are 
partly undisturbed plane waves, partly scattered 
waves. On the other hand the problem described 


by @ak, Fae is the problem of the possible states 
of the free nucleon in the absence of free mesons. 


nd 
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For the latter we have to retain the terms up to 
the order D~? but we can omit terms of the order 


D~* which are independent of oom or terms of the 
order D-? linear in Fae which are at the same time 


linear in qu or bilinear in and gp’ (xe). 
The latter terms do not have an expectation 
value of the order D~? in the state where the cor- 
responding waves are unexcited, and the terms 


independent of Pd do not influence the energy 
difference between the excited states of the 
nucleon and the ground state (isobar separation). 
These more complicated terms, however, give 
rise to a correction of higher order to the scatter- 
ing of free mesons, which question is not treated 
in this paper. In the following section we treat 
the meson scattering only as a problem of the 
order D® and disregard all terms of higher order, 
while in this section we treat the problem of the 
nucleon in the absence of free mesons up to the 
order D~. 

In order to get rid of the last terms in (II) we 
have to shift the zero point of ,’(x) but in such 
a way that the orthogonality relation (12’), 
namely =0 is not violated. The 
suitable assumption is 


4x)! 
Ta’ (x) = (x) — = = 
75) 
S (d&/dx)*dV 


which agrees with the orthogonality condition 
(12’) because of (15). Indeed, we have, again 
using (12’), 


= 


OX, 


(4x)! 0K , 


Taking into account (15) and the definition (23) 
of N we find 


OK 3 
Ox S (0&/dx)*d V 


16We omit the index on the variable x if an integral 
has the same value for all three values of this index. 


The term with N just gives the second term 
in (II); therefore we get finally, inserting 
&(x) = X(x)/J and the value (63) of D 

4n 


~ 


(1), 2 


H,=3N = (qas) 


OK 


a, 


3g 


(76) 


+4 


HE (77) 


In the third term of the shift of 2.’(x) does not 
give any contribution, because 


Lo** pag =0, 
a,B 


Lo“ being antisymmetric, pas symmetric in a@ 
and 8. 

Before we proceed, it is necessary to complete 
the discussion of the shift (75) of 2.’(x). As this 
shift depends on Lo and eg; it is a g number and 
necessitates corrections to Lo and eg, them- 
selves to insure that they commute with the new 
Ta''(x) and at the same time maintain their own 
commutation relations. This goal is reached by 
the canonical transformation 


F'e-U = F’] 


[U, (78) 


We can substitute for F the old variables 7,’(x), 
esx, for F’ the new variables Lo*®’, 
es.’ if we choose for U 


2D ape 
x| ar] (79) 


which again gives the shift (75) for the first-order 

correction to mq'(x) as a result of (21). 
Fortunately it is not necessary to give here 

the explicit formulas for the corrections to Ly?*’, 


ng 


6) 


7) 


ot 
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and the second-order corrections to (x) 
because the additional terms in the Hamiltonian 
of the order D~? are either independent of Lo** 
or linear in Lo** and at the same time linear in 
¢a'(x) or bilinear in ma’(x), gs’(x). The expecta- 
tion value of the terms linear in Lo%* is zero in 
the stationary states of the Hamiltonian (76), 


(77). It may be noted that io is only an approx- 
imate integral of the motion, but if we replace it 
in (76) by the exact integral (22’) of the total 
angular momentum, the additional terms in the 
Hamiltonian are again of the type which is dis- 
regarded here. 

Because the eigenvalues of Lo? are j(j+1) 
with half-integer 7 after the transformation (54), 
we get directly from (76) for the energy of excita- 
tion of the levels of the nucleon above the ground 
level 


3 x? 


(j+1)—3}. (00 


The charge quantum number nm does not enter 
but is restricted by (74) to values less than or 
equal to j. It is characteristic of the symmetrical 
pseudoscalar theory that the system is degen- 
erate and the energy of all states with the same 
j and different charges has the same value. For 
the ground state, j=}, the number m has the 
two possible values n=—} and n=}, which 
correspond to neutron and proton. 
The integral 


1 1 
fxaxav 
4r 


2 
= [xKav-= far, 


[see (13)] and can be evaluated in the two 
limiting cases of small sources xa<1 and large 
sources xa>1 as follows. Consider the expression 


f XKdV= f f 


In the first limit we may substitute 1 for e~*’, the 
error being of the relative order xa. This means, 
according to (4), that we put f/XKdV=a~. In 
this limit the second integral (x?/42) can 
be disregarded because it becomes independent 


of a if xa is small. Therefore we get 


1 1, 
V=-, 
a 


AE 3xa 
for 


pc? 
The expression (AE/yc*), » being the rest mass 
of the meson, refers to ordinary units of the 
energy and it is equal to AE/x« in our units. 

For large sources xa>1 it is possible to inte- 
grate (13) by the following development 


4a i\- 
K K 


which corresponds to a development in increasing 
powers of (xa)~' in the resulting integrals. The 
leading term gives [see (23) ] 


1 4n 3 
av=— (6K) vide 


Hence we get, using N=aa~* (with a numerical 
factor a depending on the shape of the source 
function) 


f (VX)%dV 


AE 1 


We have now to consider the various terms 
neglected and the conditions which are implied 
by such procedures. In the first place, we ne- 
glected terms in the Hamiltonian which had as 
factors 7,, ¢,, since these had no diagonal matrix 
elements for the system of states chosen, but 
introduced coupling between them. In order of 
magnitude, such terms are given by 


1 
av]. 


For small sources (xa<1) we get from (80a) 
h~ 


The energy of separation between the states of 
(62) is of the order g*J/x*. An evaluation of the 


{j(j+1)—3} for xa>>1. (80b) 


— 
) | 
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integral J, defined by (14), shows that in the 

limit xa<1, J~1/a*. Consequently, if we demand 

that h be small compared with the energy of 
separation we get 

xa/g<1. (81) 

The second approximation was introduced in 


(70) and (72) when das was treated as small 
compared to D; we restricted ourselves at this 
point to small oscillations of the field about its 
position of equilibrium. This restriction is not 
necessarily essential to the strong coupling hy- 
pothesis, but was introduced for mathematical 
simplicity. If it failed, the isobar energy (80) 
would have had a more complicated dependence 
on the angular momentum quantum number j. 
We estimate the validity of the approximation 


by computing ((qaa)?)n/D?; we do this for the 
case where no free mesons are present, so that 


the only contribution to (gap?) comes from 
the zero-point oscillations of the field. We find 


for small sources 
for xa<1. 


In the same limit D=gI/x~g/xa*. Once again, 
the criterion proves to be 


xa/g<1. 


Finally, we consider the unitary transforma- 
tion (78). To the first order, its effect was simply 
to shift 2.'(x) according to (75) in such a way as 
to separate the Hamiltonian completely into a 
bound part and a part (77) associated with 
unbound mesons. However, higher order terms in 
(78) spoil this separation; they give a depend- 
ence of the isobar energy on the distribution of 
free mesotrons, or conversely a dependence of 
the scattering of mesotrons on the isobaric state. 
They make possible such processes as the transi- 
tion from one isobaric state to another with the 
inelastic scattering of mesons. If we demand that 
the ratio of these terms to the isobar separation 
be small (i.e., that the half-width of the isobaric 
states be small compared to their separation), we 
arrive at the condition 


f / f (*)avec (82) 


Again we substitute for ¢,’ values corresponding 


to the zero-point oscillation of the unbound 
meson field. Evaluating the expressions in the 
small source limit, we find!” 


(xa)? In xa 

g? 
For xa <1, this condition is certainly fulfilled for 
values of g large enough to satisfy (81). For 


comparison we give the preceding expressions in 
the large source limit (xa>1): 


I~1/x*a', Dr~g/x'a', 
xa’. 


Condition (81) for the first and second neglects 
is replaced by 


K1. (83) 


(xa)*/g? <1. (81a) 


Inequality (81a) also expresses the validity of 
(82) for this limit. 


7. THE SCATTERING OF FREE MESONS 


As mentioned already, in this paper we treat 
the scattering with disregard of all terms of 
order D~ which means that we do not distinguish 
between elastic and inelastic scattering. Our 
treatment of the problem therefore will give us 
the total scattering cross section. Omitting the 
irrelevant self-energy constant, we can use here 
the Hamiltonian (77) directly, where the differ- 
ence between 7,/(x) and z,’’(x) can be disre- 
garded. The quantities e, or é., can be treated 
here as c numbers for they commute with all 
other observables in the Hamiltonian.'® It is 
convenient to go back from the variables in (77) 
with the help of (16), (16’), (70), (72), to the 
functions 


OK 
Pa(x) (47)! Pagesi——+ Ta (x) 
OX; 


1 0K 


(84) 
2D Bk Ox; 
= —1/(4)! Cax(OE/Ox,). (85) 


17 The condition from which (82) was derived also 
restricts the values of 7 for which the approximation 


holds. We find ]7<g/xa. 
18 Tf, on the other hand, one uses a matrix representation 


0 0 
where 7\2, Liz, and Le? are diagonal, the éax are matrices 
which are a generalization of the rotator matrices already 
used by Heitler and Ma (reference 8). 
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We see that pa(x), ga(x) represents a field of 
free mesons, the meson eigenfield of the nucleon 
being subtracted. It fulfills the commutation 
relations [see (21), (55) ] 


i[Pa(X1), 9a(X2) ]= 5ag5(X1— Xz) 


OK @ OK @ 
~ Cak€pi—— =. (86) 


OX2x OX1; OX2% 

These relations camel to the fact that the 
symmetry of the as and Pas makes the field 
Pa(X), Ya(x) fulfill the subsidiary conditions [com- 
pare (12), (12’)] 

— ppd f v| =0, (87a) 

OX; OX). 


k 


jeu f mt 0. (87b) 


On the other hand, the Hamiltonian (77) 
becomes formally a Hamiltonian for free plane 
waves 


f f V. (88) 


The fact that the equations of motion differ from 
the wave equation of free particles is entirely 
due to the symmetry conditions or to the 
difference of the commutation relations (86) 
from the usual canonical ones. 

From (87a, b) and (88) we get, by the method 
of Lagrangian multipliers which have here the 
form of antisymmetric constants Aas = —Aga and 


= —Aga, the equations of motion 
Ga=Pa— Axx), 
Bk 
— Pa= (—A+k?)ga— /Ax;). 
8k 
With the help of the time derivatives of (87a), 


(87b) in which one has to insert the expressions 
for Ga; Pa One gets 


aK 
Aap = 3 f (Palsx ppax)—d V, Aas = 0. 
k OX: 


Therefore we get the equations of motion in 


1° The analogous result in the charged scalar theory was 
first obtained by J. Schwinger. 


the form 
a 
OK 
x f (Peta: (89a) 
(89b) 


The same equations follow from (86) with the 
help of the relation f=iLH,f] for f=p,. and 
f= Qa- 
If we first consider ga, Pa as Classical quantities, 
we can integrate (89a, b) by putting 


Pa= tla, da = 3 


OK 
x f (90) 


where ue has, according to (89b), to fulfill the 
wave equation 


— ua =— 
2I 


0K 
x f (91) 
Ox; 
and the subsidiary condition 
0g 
> Jen f | V| =0. (91a) 
k OX, Ox, 


This condition is, however, not independent of 
the wave Eq. (91), for the second time derivative 
of (91a) vanishes as a consequence of (91). 
The subsidiary condition essentially excludes 
only the static solution 


k 
The quantum-theoretical formalism leads to 
an analogous result. Let u.‘" be a normalized 
orthogonal set of functions (c numbers), which 
fulfill (with the usual boundary conditions in a 

hole) 
aK 


2I 


0K 
X f(a (91b) 
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and the subsidiary condition (91a). We put 


4 
= A (=) (92) 
0 
=2 A,(2wn) (n) _ >> 
4 f ~<a}, 
Ox; 
Ain=Ant, =ug™*, win=wn>0. (92a) 
If we assume 
[An, Am*]=[An, A—m]=68nm, (93) 


we get just the commutation relation (86), the 
system u,‘” being incomplete because, as a result 
of (91a), the static solution 
is missing. Moreover, the Hamiltonian (77) 
assumes the normal form 


1 n(AntAnt}). 


Our scattering problem is therefore reduced 
to the determination of the solution of (91) 
for, aside from a factor exp (—iwt), of (91b)] 
which describes a plane wave incident in a given 
direction and an outgoing scattered wave. 
Corresponding to the three pseudoscalar fields 
Ua (a=1, 2, 3) we shall get three linear inde- 
pendent solutions or of this 


(94) 


problem, which correspond to an incident 


particle with p=1, p=2, or p=3. 
Instead of this we can also use the solutions 


Uat =(1/V2) (ua? +ita™), 
Ua = (1/V2) (ua —ita™), ta, 


which correspond to an incident positive or 
negative charged or neutral particle. 
We now solve (91b) by the assumption 


Ua” (x) = bay exp [ik(n-x) ] 


i,k 


x f (95) 
Here p=1, 2, 3 enumerates the three different 
solutions, n is a unit vector in the direction of 
the incident wave, k = (w?—x*)!, and the constant 
A has to be determined. The first part, the 
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plane wave, fulfills the wave equation for free 
particles without the right side of (91) while the 
scattered wave and therefore the whole ex- 
pression (95) fulfills the equation 


ane (A+k?)u, 
= 4rA Zz. (SapMk i) (6K /dx,) 


i,k 
In order to check (91b) we introduce the ab- 
breviations 


(95’) 


fro exp [ik(n-x) ]dV 


sinkR 
K(x)——dV=F, (96) 
kR 
(where R= |x|, r= |x—x’|) 
0K(x) cos kr K(x’) 
ff —dVdV'=6x%J. (97) 
Ox; r Ox)” 
Further we have to evaluate the integral 
0K in kr 0K(x’ 
f f (x) sin kr 0K (x vay" 
Ox; r Ox," 
f f \dVaV" 


But as a consequence of the fact that sin kr/kr 
has no singularity for x=0, (sin kr/r)K(x’)dV’ 
is a solution of the wave equation (A+?)u=0, 
and therefore 


f (sin kr/r)K(x')d V’ =(sin RR/R)C. 
Moreover, C=F because for R=0 the left side 


of the equation is equal to RF, the right side 
equal to kC. Therefore 


0K (x) sin kr 0K(x’) 
ff dVdV'=6,—F?, (98) 
Ox; r Ox,” 3 
and finally 
OK kr OK (x’ 

ff (x) 

Ox; 
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With the assumption (95) we get 
0K 
k OX), 
= Le — 5 
The differential Eq. (89b) therefore gives 
(A+k*)u LX (Mr bap — Core 
a = pkOailli ax, 
x | —ikF+2A (4+5F) 
and the comparison with (95’) gives 
ikF 
; (100) 
ikF/2 


ik 
J-I 
Introducing 


tan 


A = (3i/2k?F) sin de®. (100’) 
We get the asymptotic behavior of the scattered 
wave for very large r in the direction n, of the 
scattered wave by putting r= R—ck(n,-x), 
V=-—v7kn,, 
Jecatte ~ tk FA [6,(n-n,) 
—(€q-n)(e,-m,) J(e*"/R). (102) 
We now compute the square of this amplitude 
for a given n, and sum over a, 
Dea (x) =k? *[(n-n,)? 
—2(n-n,)(e,-n)(e,-m.) +(e,-n,)*]. 
Further, we take the average over-all directions 
of the axis system e, and get for the right side,”° 
with cos @=(n-n,), @ being the scattering angle 
(kFA)*}(1+ cos? 6). 


*° This is in the approximation where the energy loss of 
the meson in the scattering is disregarded, the same as 
averaging over-all directions of the angular momentum 


0 
Lix in the initial and in the final states. 


This is the cross section dg for the solid angle 
dQ= 2x sin 6d6; inserting (100’) we get 
4 cos? (103) 
== —(1+ cos 
4 


and the total cross section 
q= (sin? 5/k?)4x. (104) 


On the other hand, if we integrate the square of 
the amplitude over all angles 6, which means 
averaging over n, without summation over a, 
we get 


(kFA)?3[bap— 25ap(€,-m) +(€a-n)*]. 
Averaging over the orientations of the axis 


system gives 


(kFA)*}(5ap+1) 
which means 
Ja = 4n(sin® 5/k®)i(1+ 8a). (105) 


In other words, if the incident particle is a 
“1 particle’ the number of scattered ‘1 parti- 
cles,” ‘‘2-particles,’’ and ‘‘3-particles” is }, }, }. 
This also holds if one takes for ‘‘2,”’ 
in any order, neutral, positive, and negative as can 
be shown by using the above-mentioned eigen- 
functions 


Uat = (1/V2)(ua +iua™), 
Ua = (1/v2) —it,”) 


as well as u,.. 

Finally we evaluate the integrals J, J, and F 
for the physically important case ka<1, xa<1. 
In this case we can first put F=1. Further we 
have, according to (97), 


3J= f —A)(cos kr/r) |K(x’)d VdV’ 
=4r f K*(x)dV 


+k? f f K(x)(cos kr/r)K(x’)dVdV’. 


In the same way 


f f K(s)| 


nite f K*(x)dV—«? f f 


r 


Therefore 


(x/)dVaV". (106) 
r 


So far this result is exact, but for ka<1, xa<1 
we can substitute 1 for cos kr and e~*’ and get, 
with w?=k?+ x? 


3(J—I)=w?/a for ka<1, xa<1, (106’) 


hence 
sin 6~tan 6~ (k*a/w*) 


and for (103), (104) we get finally on writing 
v/c=p/E for k/w 


dq= (=) +cos? 6)dQ 


q= (<) a*4r 


which was quoted in the second section. 


(107) 


8. THE CHARGED PSEUDOSCALAR THEORY 


Although the symmetrical theory is, from a 
mathematical standpoint, more elegant than 
the ‘‘charged pseudoscalar theory’’ which uses 
only positive and negative charged, but no 
neutral, mesons, it is also of interest to develop 
the corresponding results for the latter theory, 
because neutral mesons have never been observed 
and because in the case of strong coupling the 
charged theory also leads to forces between two 
like particles. 

With a complex pseudoscalar field g the 
Hamiltonian of this theory can be written 


ee f d V| (Ve) (Ve) 


|, (108) 
K 


where 
(109) 
In order to facilitate the comparison with the 


symmetrical theory we split + and ¢ into the 
real and imaginary parts which are normalized 
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according to 


1 1 


(110) 


1 1 
—im), g*= ties)» 


[= a a a 


(ar) f K(x)rao-Vead V}. (111) 


We can therefore take over the whole develop- 
ment of Section 3 if we simply drop the field 
components 73, ¢3 and let the summation index 
a take only the values a=1, 2. Of course there 
exists as an integral of motion only the 72 
component of the isotopic spin in this theory, 
which gives the electric charge minus 3. 

The determination of the eigenvalue of the 
interaction energy 

3 


g 
Him=— (rnent 
KV2 


here is somewhat simpler than in the case of the 
symmetrical theory. Computing the square of 
the interaction energy, we get, 

= (¢1°)? + 91°, 


The four eigenvalues of the interaction energy 
are therefore 


= + (449) 
K 
— $2")? (112) 

with all possible combinations of the + and 
— signs. ° 

The application of the formalism of Section 4 
leads here to the case where one of the eigen- 
values Q of (C) [Eq. (30)] or (C’) [Eq. (30’)] 
is zero because os: =0. The matrix C’ has 
therefore only two rows and two columns and is 


| 


10) 


1) 


»p- 
eld 
lex 
ore 


as 


he 
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given by 
C= (¢1°)? 
and Q,*, Q.? are the roots of the equation 


(G1°)? —x} | —x} — bo")? =0, 


which gives 


1 
(1°)? + (2°)? [ (1°)? + (62°)? 
— (1° 2")?) (113) 


The expression (36) for the interaction energy 
g 
Eine =—(4014Q2) (112’) 


is in agreement with (112). The orthogonal 
matrix B,. has here only two rows and two 
columns and can be assumed as 


cos @ sin 


—sin@ cos@ 


(114) 


If we write n™ for A;, the normal form (29) for 
gives us 

Qi cos sin én®, 

sin én’) +Q, cos én. 
The significance of n“, n® is therefore that n™ 
is orthogonal to the plane through $,° and 2° and 
the orientation of the n™, n® axes which we may 
call for the moment the é and 7 axes is defined by 
the condition 

0 0 0 
P1EPig t 0. (116) 
This can be immediately derived from (115) and 
means that the matrix C defined by (30) is 
diagonalized by going to the (é, 7, ¢) axes, the 
diagonal matrix elements being Q,, Qe, 0. 
In the Eqs. (45), (45’) we have to put 

ru =0, B3a=0, B,3=0 for a,r=1,2; P3=0; 


hence 


(115) 


(Lo!?+ T!2)? 
2°) 
2(Qi—Q2)? 
Lo")? 


2(Qi+Q2)? 


21 We denote corresponding equations in the symmetric 
and charged theories with the same number and add in 
the latter a c. 


The symmetrical! quantities gas, Pas have to be 
defined only for a, 8=1, 2 and according to (51) 


we get 
=n" cos sin 8, 


=n sin cos 8, (117) 
and from (52) 
qu=Q: cos? sin? 6, 
qi2=21= (Qi—Qz) sin (118) 
Q; sin? cos? 6. 
For the description of n™ and Lo” by Euler angles 
and their conjugate momenta, see Appendix. 
For the minimum of the potential energy we 


get from (62) the old value (63) for Q:=Q.=D; 
but because Q;=0 we now have 
D? 1 g? 
I 2«* 


In the approximation where que is disregarded, we 
now get 
Lo'* = D{ — 
= — Lo? = (a=1, 2). 
Therefore 
2D D 
hence 


0 1 
(x.°)? =) (Pap)? +-—— 


1 
3 2D? 


2, 


This is in accordance with (50c) if we put in the 
last two terms Q:=Q2=D, —Lo". In (II) 
we have therefore to substitute 2Lo?—(T»"*)* 
for 

In the same way we have to substitute in (75) 


1 1 
— iL aB 2e iL as f oes é iL af 


and get instead of (76): 
4n 


1 g? 
—-—I+H, (76c) 
2 x? 


| 

| 

| 

| 

: 

| 

| 

= | 

| 

he 

of 

| 

4 

n- 

is 
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where H/; is again given by (73), and instead of 
(80), (80a), (80b) 
_3 


2j(j+1) —n?—5/4], (80c) 


AE 3xa_ 
—n?—5/4] for xaX1, (80ac) 


pc? 


= n?—5/4] 


2 
for xa>1. (80bc) 


The restriction (74) —jSnSj is of course 
unchanged. 

In the computation of the scattering the 
following modifications of Section 7 take place. 
The functions p(x), gs(x) (a, 8B=1, 2 only in the 
following) fulfill, besides (87a, b), the condition 


es VKqd V=0, ex: VépadV=0, (87c) 


and the commutation relations (86) get a 
corresponding additional term: 


pa(X1), 9a(X2) ]= —X) 


1 x, OK 
2 Ax : OX1; OX2% 
1 OK 
— —. (86c 
With 
Pa=Ua, 


Ja = Ua =f Ualsi€sk 
k aK 


we now get 


es: Véuad VeusdV=0, 


ex: f V=0, (91c) 


21 ik, 8 


aK 
x f Bint esieax) — dV. (91c) 


The solution is 
Ua‘? (X) = dq exp [tk(n-x) ] 
+A ] 


eikr 


(n-e,)(e,-¥)} f (95c) 


where A is again given by (88), (89). Asymp- 
totically for large R we have 


]scatt.~ tk es) es) 


—(n-e,)(n,- 
>. | (x) |2scatt. =(kFA)?{[(n- n,) 


a=1,2 
+(n-e3)(n,-es) ]?—2(n-n,)(n-e,)(n,-e,) 
The process of averaging over the directions of 
the axis system @1, @2, @3 gives (for any value 


of p, and -for the solutions which correspond to 
an incident charged particle) 


dq = (sin? 6/k?)?(1+3 cos?@)dQ, (103c) 
q= (sin? 6/k*) 34x, (104c) 


where for sin? 6/k? we can again substitute 


if ka<1, xa<1. 
If the incident particle has a + charge, we 
find for the total cross section (integrating over 
all scattering angles) for scattering with the 
same (+) and the opposite (—) charge, re- 

spectively, 
= (sin? 6/k?)42(5/4), 


q-= (sin? 6/k?)4r(1/4), 


(105c) 


the sum of both agreeing with (104c). 


9. THE MAGNETIC MOMENT OF THE NUCLEON 


The meson field surrounding the heavy particle 
gives rise to a magnetic moment which will prove 
to be of equal magnitude, but of opposite sign 
for the proton and neutron, just as in the case 


ve 
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of the weak coupling theory.” The results of 
this section are equally valid for the symmetrical 
and charged pseudoscalar theories, since the 
neutral mesons give no contribution to the 
magnetic moment. 

The magnetic moment is related to the current 
density S by 


Mutsy f [x-Slv (119) 
or, in tensor notation, 


e being the absolute value of the elementary 
electric charge. The definition of the charge and 
current densities Sp and S becomes essentially 
ambiguous in a theory which uses finite distance 
operators, such as the source function K(x). 
Indeed, in such a theory it is not possible to 
fulfill the continuity equation 


dSo/dt+div S=0 (120) 


everywhere inside the source. It is fulfilled 
outside the source, as well as (on the average) 
over a volume V of dimensions of order a in 
which the source is included. 


[AS)/dt+div S]dV=0. (120a) 


Expressions for Sy and $ which fulfill this weaker 
condition and which agree with the definition (9) 
of the total charge are given by 


So= +3(1+73) K(x), (121) 
S=(g2 grad ¢: grad 


+5 (x), 


To prove that these satisfy the conservation 
conditions, we make use of the relations 


(-a+ gat" K)=0, 
or? K 


10 
2 at K 


grad gi—71 grad g2)K(x)d V. 


2H. Froéhlich, W. Heitler, and N. Kemmer, Proc. Roy. 
Soc. A166, 154 (1938). Their result, given in Eq. (63) of 
their paper, agrees formally (apart from a numerical 
factor) with our result, Eq. (124) below. 


These follow from the general rule éf/at=i(H, f] 
and the commutation relations (2), (6) applied 
to the Hamiltonian given by (1), (8) with 
f=ma(X), ba(X), 73. We obtain 


as 
+K(X)e- f dV grad grad o:)K(s). 


Therefore the continuity condition is fulfilled 
both for points outside the source, where 
K(x)=0, and on the average over the source as 
indicated by (120a). As a result of the lack of 
validity of the continuity equation inside the 
source, we may apply the assumed expression 
for the current density to the computation of the 
magnetic moment only if this moment is gen- 
erated in a volume large compared with the 
dimensions of the source; in other words, we are 
restricted to small source size, xa<1. 

Introducing the strong coupling approxima- 
tion, we apply the S transformation (65a, b) to 
the second part of the current (121b); then we 
insert the values (68) for 7,0,, namely 1,0, 
= —6,,. This gives 


S=(¢2 grad grad 


(2m) — gre) K(x). (122) 
K 


We insert for the field its value corresponding to 
the absence of free mesons and the lowest 
eigenvalue of the interaction energy, (16), (64). 


1 
Ga(x) )-grad (a=1, 2). 


If now we use 


g an 
[Eq. (63) ], Eq. (15), 


The current becomes 
OXm 


g? 3 1 
S —(ee,, —e 
2x? 4a 


ox 
(e,Pe, )—K (x) 
Ox) 


| 
| 
| 
np- | 
e3) 
| 
of 
lue | 
3c) 
| 
=e 
Pa(x) -grad X. 
gn 
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and the magnetic moment (119a) becomes 


( 0?X ox 
x: — XE 
Ox OXm 


ox 
— — (x) . (123) 
xi 


The first volume integral is transformed by 
partial integration with respect to 0/dx,, half 
the terms vanishing in the sum over / and m. 
In the second integral one substitutes for K(x) 
its expression (47)~!(—A+x«?)X taken from (13) 
and gets 


ox ox 
fx, 
Ox) 
=} 


445 avec (x dX 


Ox, Ox10x, Ox) Ox, 


f XAXd f XV 


dX 
Ox) Ox; OX10X; 


2 
= — f (grad X)%d V— bu f XV 


Using the spherical symmetry of X we have 


OX 6X 
———dV= f (grad X)dV. 
Ox; Ox 
We get 
e 
8a 2x? 


3 f xyar+e fear] 


As was shown in Section 4 the second integral 


can be neglected for xa<1 while 
1 
f (grad X)*dV =1/a. 


Changing to vector notation and_ using 
fe” xe® ]=e® we get the final result 


1 
M=—-~——e®, (124) 


We now compute the component of M in a 
fixed direction, say along the x3 axis. We 
introduce a matrix representation for the e,‘, 
where Li.°=Z;°, (Lo)?, and [cf. 
Eq. (73)] are diagonal with eigenvalues m, 
j(j+1), and n, respectively, with —j=mSj and 
—j=n=j. From the commutation relations 
(39), (55), (56) it follows that L = ZL" commutes 
with all LZ; and is therefore diagonal with respect 
to j, m, and L,. 

Moreover, e® commutes with L® and is 
therefore diagonal with respect to ». However, 
e;® is not diagonal with respect to 7 and we 
need the diagonal element (j, 2, m|e3“| 7, n, m). 

We have the following relations: 


(m | 

(m | 
To satisfy the commutation relations [analogous 
to (38) ] of e, with ‘ we set 

(j, m, n\es® | m,n) =(j,n|C| j, nm, 
(j, m, n|e: | 7, m—1, n) 
=(j,n|C|j, ]}, 
(j, m, n |e: —te, | 7, m+1, n) 
Inserting this in the equation 
—(L-e®) = —L® =T® 


—(j,n|C\|j, =n 
(j, m, n\es | j, =—— 
—(j,m,n\e ,m, n) =— 
iG+) 


Therefore, for a given isobar, we find that to the 
leading order 


_ eg? 1 nm 


x xa 12 iG+1) 


(125) 


| 
OX OX 
| +f 
Ox, OX; 
| 
| 
| 
| _ = 


ous 


)}'. 


the 
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For the proton n=m=j=}, hence 


(126) 
Ka 36 
For the neutron n=—}, hence M; simply 


changes sign, as was mentioned above. 

If, for a first orientation, we identify this 
magnetic moment with the total empirical 
magnetic moment of the neutron—1.93 proton 
magnetons—and assume a meson mass equal to 
1/10 proton mass (e/k=10 proton magnetons), 
we find g*/xa~6.95. Inserting the value of 
xa~0.1 determined from consideration of the 
scattering, we find g?~ 0.695, or g/xa~8.3, well 
within the range of the strong coupling 
approximation. 

However, the fact that the proton and neutron 
magnetic moments are not of equal magnitude 
cannot be explained in any simple way in this 
theory. Of course, (126) does not give the total 
magnetic moment, but only the part due to the 
meson cloud. To it must be added the moment 
due to the nucleon itself.2* If we assume that a 
“bare’’ proton or neutron would have a magnetic 
moment or one or zero proton magnetons, 
respectively, then we can compute the nucleon 
moment in the two lowest states by calculating 
the expectation value of 


eh 1 
M’=—— ~(1+73)e. (127) 
m,c 2 


In the weak coupling theory 73 is diagonal with 
the value +3 for the proton and —} for the 
neutron, giving rise to the same values for the 
moment which were assumed for the bare core. 
However, in the strong coupling theory we have, 
after performing the S transformation, 


e l 
M’= —— -e®), 
myc 2 


the expectation value of @ being zero. We 
compute the component of M’ along the x; axis 
as before and get 

1 mm e 


mye 


*% We are indebted to Dr. J. Schwinger for valuable 
discussion of this question. 


(128) 


This moment, while small compared to the 
empirical moments of proton and neutron, is 
again equal in magnitude and opposite in sign 
for the two. For the two states n= +4, m=j=} 


M,;'= +3 proton magneton. 


APPENDIX 
1. Introduction of Polar Angles 


It is possible to express by means of Euler angles the 
orthogonal matrices (or, in other words, the three or- 
thogonal-unit-vectors) introduced in Section 2. For 
instance, we can put 


=cos cos ¢ cos ¥—sin ¢ sin 
cos @sin gcosw+cos gsiny, —sin @cos y, 
Ax, = —cos 6 cos sin ¢ cos y, (1) 
—cos @sin gsiny+cos gcosy, sin @sin 
Ayx=sin @cos¢g, sin@sing, cos 


which means that the right side is the transposed matrix 
A. In the same way we put for the matrix 


Bia =cos @ cos cos ¥—sin sin ¥, 
cos @ sin ®@cos ¥+cos@sin ¥, —sin Ocos ¥, 
Boa = —cos cos ® sin ¥—sin cos ¥, (2) 
—cos @ sin @sin ¥+cos®@cos¥, sin Osin ¥, 
Bzsa=sin @cos*%, sin@sin%, cos9O 


and for the éax defined by (51) 


€i%=cos a cos b cos c—sin sin ¢, 
cos asin bcosc+cos bsinc, —sinacosc¢, 
€2%: = —cosacos b sin c—sin b cos ¢, (3) 
—cos asin bsinc+cosbcosc, sina sin ¢, 
éx=sinacosb, sinasin’b, cosa. 


Defining the operators pe, py, py simply by pe= — 10/08, 
pe=—id/dg, py=—id/dy, we find from the volume 
element sin 6d0dgdy that py, py are Hermitian, but 
pet = (sin sin 6; 
similar expressions are found for poet, pat. 
A direct evaluation of the angular momentum operators 
gives 
L,=L23= —sin gpe+(cos ¢/sin 6)(py —cos 
= — pot sin ¢+(py— py cos 8)(cos ¢/sin 4), 
L2=L3:=cos gppt+(sin ¢/sin @)(py—cos Opy) (4) 
= pet cos ¢+ (py — Py cos 8)(sin ¢/sin 4), 
L3=Li2= py, 
and also 
L, =L23= —sin bpa+ (cos 6/sin a)(p.—cos ap») 
= — pat sin b+ cos a) (cos b/sin a), 
L2=L3:=cos (sin b/sin a)(pe—cos apr) (5) 
= pat cos b+(p.— p» cos a)(sin b/sin a), 
L;=Li2= po. 
Analogously for the isotopic spin 


T1=T23= —sin + (cos @/sin @) (py —cos Ope) 
= — pot sin + (py — pe cos @) (cos /sin 9), 


| 
na 
a) 
[cf. 
| 
| is | 
m). 
| 
: | | 
| 
{ 
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T:=T31=cos (sin &/sin @)(py—cos Ope) (6) 
= pel cos 6+ (py — pe cos (sin &/sin 9), 

T;= Ti2= po. 

For the components L”* and 7°* we find 


L'=L%=sin ypet+(cos ¥/sin @)(cos Opy — py) 
= pet sin ¥+ (py cos py) (cos y/sin 4), 
L?=L*"=cos ppe—(sin ¥/sin @)(cos — py) (7) 
= pet cos (py cos py) (sin ¥/sin @), 
D=L"=p,, 
T!=T*=sin Vpe+ (cos V/sin @) (cos Opy — pe) 
= poet sin ¥+ (py cos O— pe) (cos W/sin 
T?=T*®=cos (sin ¥/sin @) (cos Opy — pe) (8) 
= poet cos ¥—(py cos O— pg) (sin W/sin @), 
T?=T"%=py. 
And, analogously for the components L° which we write 
L®, L®, L® or to distinguish them from 
L's 
L® =L®) =sin cha+(cos c/sin a) (cos ap.— pr) 
= pat sin c+(p. cos a— (cos ¢/sin a), 
L® =L®) =cos cha— (sin c/sin a) (cos ap.— (9) 
= pat cos cos.a— ps) (sin ¢/sin a), 
L® =L® = p,. 


We note that all components L", 77*, L% turn out to be 
Hermitian. For the total square of L and T we get 


= pot py + O08 


sin? 0 
6 2 
2 


= Pot po + (11) 


and also 
Lt = pat pat 
ante. +o Po)” 2, (12) 


sin? a 
The expressions (7), (8) can be substituted in the 
expression (50) for 


0 
(mat)? 
a, k 


which, however, gets rather complicated even in the 
approximation where one puts Q,+Q,=2D because of the 
terms containing Q,—Q,. The analogous expression (50c) 
for-the charged scalar theory is simpler, because here we 
have simply 7"°=T7,.=p@. In the approximation where 
we can put in the last two terms Qi+Q2=2D, Q,=D 
we get 


+ (a2)? = Pit Pot P2+ + 


sin? 


2. The S Transformation 


By using the expression (3) for the éax the S transforma- 
tion defined by Eq. (65a) is simply obtained by putting” 


S=exp (ies) exp (io:5) exp (ies3): 
exp i035) exp ) exp 
(-ie:5) exp exp (14) 
2 2 2 

which can also be written out as 

S=cos $(cos 5 cos $—sin 5 sin 5) 7 


+i sin 3(-sin cos 5 toos ; sin 


+i sin cos cos stein sin 


+i cos 5(sin cos 4 sin (15) 


The expression for S~! is obtained from it by changing i 
into —i. One easily finds 


SpaS = pa— sin cos 


SpoS = po— sin a cos 
+o2sinasinc+oe;cosa), (16) 
= pe— ho. 


Inserting this in the expressions (5) and (9) of the last 
section, one easily checks the Eqs. (69), (69') of the text, 
namely 

gala, 
5-1 = LB — 

See W. Pauli, Helv. Phys. Acta 12, 147 (1939), 

Section 3, where the special case c=0 is treated. 


| 

| 
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In studying the fine structure predictions of any pro- 
posed theory of nuclear forces it is important to know to 
what extent the results obtained are sensitive to the 
changes in the nuclear model employed in making the cal- 
culations. The present work is concerned with this. The 
splitting for He®, which should be very similar to that for 
Li®, is discussed on the alpha- and central field models. 
A deuteron group model for the alpha-particle is con- 
sidered in addition to the simple alpha-group. Li’ is treated 
from the standpoint of the theory of “holes” and also on a 
modified alpha-particle model. The contribution of the 
exchange terms to the splitting in He‘ is greater by a factor 
of the order of two or three on the alpha-model than on 
the central field model. The sign and magnitude of the 


total splitting are of about the same order on the two 
models. The deuteron group alpha-function gives essen- 
tially the same splitting as the ordinary alpha-function. So 
far as ordinary terms for Li’ are concerned, the spin-orbit 
interaction of a “‘hole’’ in a concentrated closed shell with 
a closed shell is the same, apart from a numerical factor, 
as the spin-orbit interaction of a single particle with a 
closed shell. The exchange terms for Li’ are approximately 
proportional to the degree of overlapping of the alpha- and 
triton groups. The exchange splitting has the same sign 
as the ordinary splitting. Energy splitting estimates on the 
alpha-model in which the exchange terms are neglected 
should give results reliable to an order of magnitude at 
least, as far as this model is concerned. 


I. INTRODUCTION 
ARIOUS experiments suggest that the 


lowest states of the nuclei He’, Li®, and Li’ 


are doublets. Calculations of Feenberg and 
Wigner! indicate that the ground states of these 
nuclei should be ?P. It is logical to interpret the 
doublets as consisting of the and *P 1/2 parts 
of a *P state. 

The evidence relating to He® comes from the 
experiments of Staub and Tatel? on the scatter- 
ing of neutrons by helium. The variation with 
neutron energy of the backward scattering cross 
section indicates a doublet structure with a 
splitting of about 0.3 Mev. The sign of the 
splitting is not fixed. 

Experiments by Heydenburg and Ramsey* on 
the scattering of protons by helium show a very 
broad maximum in the scattering of half-width 
over 1 Mev. There is no positive evidence of a 
doublet structure for 

In the case of Li’ the ground state is stable, 


* Submitted to the Graduate School of the University 
of Wisconsin in partial fulfillment of the requirements for 
the degree of Doctor of Philosophy. 
ae Now at Naval Ordnance Laboratory, Washington, 
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1077 (1940). 


and its spin is*~’? 3/2. Rumbaugh and Hafstad® 
found an excited state of this nucleus about 400 
kev above the normal level. 

Theoretical investigations of the fine structure 
of these nuclei can, as pointed out by Breit and 
Stehn,’ throw light on the form of the interaction 
energy between pairs of nuclear particles. Also, 
since the fine structure matrix elements involve 
the derivative of the potential and hence may 
magnify the influence of the shape of some poten- 
tials at small interparticle separations, studies of 
fine structure offer a possible way of investigat- 
ing the potential at small distances. It is difficult 
to obtain this knowledge from scattering experi- 
ments on elementary particles, since, according 
to Hoisington, Share, and Breit,’® scattering 
results depend mostly on the shape of the outer 
region of the potential well. 

In testing the fine structure predictions of any 
proposed theory of nuclear forces it is important 
to know to what extent the results obtained are 
sensitive to the changes in the nuclear model 
employed in making the calculations. The split- 
ting for He*®, which on present ideas should be 

*H. Schuler, Zeits. f. Physik 42, 487 (1927). 

5 A. Harvey and F. A. Jenkins, Phys. Rev. 35, 789 (1930). 

6 L. P. Granath, Phys. Rev. 42, 44 (1932). 
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8 L. H. Rumbaugh and L. R. Hafstad, Phys. Rev. 50, 
681 (1936). 

®G. Breit and J. R. Stehn, Phys. Rev. 53, 459 (1938). 


LL. E. Hoisington, S. S. Share, and G. Breit, Phys. 
Rev. 56, 884 (1939). 
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very similar to that for Li®, will be discussed on 
the alpha-particle and central field models. Li’? 
is treated on a modified alpha-model, the central 
field model having been discussed previously by 
Breit and Stehn.°® 

Dancoff" has given numerical estimates for the 
splitting in He® on a central field model, using 
the relativistic spin-orbit coupling and the 
tensor spin-orbit-spin coupling of meson theory. 
Inglis* has made estimates of the spin-orbit 
coupling on the alpha-model of light nuclei, 
including Li’. 

The theoretical basis for the alpha-particle 
model has been thoroughly discussed by Wheeler™ 
in his papers on the resonating group model. For 
He® the relatively large binding energy of an 
alpha-particle suggests that the most important 
grouping of particles will be He*+', while for 
Li’ it will be He‘+H’. The internal motion of the 
particles in each group is described by a wave 
function depending on the interparticle dis- 
tances; the relative motion of the groups is 
described by a function depending on the vector 
from the center of mass (c.m.) of one group to the 
c.m. of the other; the total wave function is the 
product of these functions. 

To simplify the treatment of Li’ the alpha- 
model is modified by introducing central field 
functions for the particles in the alpha-group and 
for the center of mass of the triton (H*) group. 
On the basis of this model it will be shown that 
the contribution to the splitting from the ex- 
change terms is approximately proportional to 
the degree of overlapping of the two groups, for 
small overlapping. This result should have 
application in approximate discussions of heavier 
nuclei. 

Breit has shown that the requirements of 
relativistic covariance of the wave equation for 
nuclear particles determines the spin-orbit inter- 
actions apart from certain adjustable parameters. 
From a consideration of the possible forms of 
wave equations for particles with spin he obtains 
the interaction energy 

H=> (1) 
— Vi Jur], (2) 
"S. M. Dancoff, Phys. Rev. 58, 326 (1940). 
2D. R. Inglis, Phys. Rev. 56, 1175 (1939). 


13 J. A. Wheeler, Phys. Rev. 52, 1083, 1106 (1939). 
4G. Breit, Phys. Rev. 51, 248 (1936). 
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for Wigner forces. The particles in the nucleus 
are referred to by subscripts k,/; ry: is the 
distance between particles k and /; the potential 
energy between two particles is —J(r;:); the 
momentum operator is written p,y=(h/i)V. >-’ 
refers to sums over pairs of particles, and > to 
sums over all values of k and / except k=1. The 
constant a is arbitrary. For particles interacting 
through the electromagnetic field the value of 
a is —1. The picture of each particle producing 
a scalar field defined in the reference system 
moving instantaneously with it corresponds to 
a=1. In the present work only a= 1 is considered. 


Note on Notation 


W4 is written for a complete antisymmetrized 
wave function, including spin; V is the unsym- 
metrized product wave function, including spin; 
y is the space part of ¥. 

R, Q, KR, p have different meanings on dif- 
ferent models and must be watched carefully. 

P;; is the permutation operator, and inter- 
changes 7 and j in the expression to its right. 

So(12) = (ai182—a281)/+/2, where a and B are 
the usual spin functions. @ is for spin up, 6 for 
spin down. Indices 1, 2, refer to the two particles. 


etc.; 


II. FINE STRUCTURE OF He‘ 


The resonating group wave function assumed 
here for the p32 state of He? is 


W= (1 — (1 +2+3+4)/4] 
Xu(rie, 713, 714, 723, 724, 134) S0(12)So(34)a(5). (3) 


Particles 1 and 2 are protons; 3, 4, and 5 are 
neutrons. « describes the internal motion of the 
alpha-particle group and is an s function: unless 
specifically mentioned to the contrary u will be 
assumed symmetric in all particles in the alpha- 
group. The function v describes the motion of the 
outer neutron relative to the center of mass of 
the alpha-group. v is a p function as is explained 
by its dependence on x, y, z in Eq. (4). 
Introducing 


o=5—(1+2+3+4)/4 


= (Xp, Vo Zp) 


q 
5 
4 
j 
i 
4 
‘ 
4 
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in Cartesian coordinates, one writes 
=(x,+iy,)R(p), (4) 


where R is spherically symmetric. 
The first-order perturbation energy for the p3;2 
state is 
(AE)3/2= (4, H'W4)/(¥4, (5) 


Neglecting the contribution from exchange terms, 
one finds 


ordinary (x°/4) (1234) R?(p) 
(6) 


a form which shows clearly the dependence on 
the radial component of the force on the outside 
particle. A more convenient form is obtained by 
integrating the original equation by parts. 

The total AE is three times the displacement of 
the state: 


(AE) ordinary = — (3X27 fu?(1234) Js:[2R?(p) 
+ (4/3)R(p)pR’(p) dr 
(2/3) fu?(1234) p?R*(p)dr], (7) 


where only the ordinary term has been included 
in the numerator as well as in the denominator. 


Evaluation of (AE) ordinary for Gaussian Functions 
Let 


J =B exp [—br?], (8a) 
u(1234) =exp (8b) 
R(p) =exp [—6p?*], (8c) 


where B, b, a, and 8 are arbitrary constants. 
The following general form of integral is needed: 


— —d — der Jd 
=n 
+6(asc+aia) (9) 


where 


This integral was evaluated with the help of the standard 
form 
S exp (—px*gx = (x/p)' exp [¢?/4p]. 
The integrand in (9) is separable into the form f(x) f(y)f(2), 
so that it is sufficient to evaluate /f(x)dx,dxedx3dx, and 
then take the cube of the result. 
The integral can also be expressed as (see Fig. 1) 


x®[all quadruplets of lines except those forming 
triangles and quadrilaterals (10) 


There are 125 terms in the square brackets. Other integrals 
used in this work are obtained by differentiation with 
respect to the constants in the integrand. 

One obtains 


(AE) ordinary = —6X*Bb[ (32a8) /(16ab+ (32a+3b)p) (11) 
Exchange Terms on the Alpha-Particle Model 
The exchange integral in (5) simplifies to 
JS v*(@)u(1234) (2B — 2Bys? w(E)u(1235)dr, (12) 
under the special assumption that w is of the form 
exp (13) 


The evaluation of these integrals with Gaussian wave 
functions and potentials is exceedingly lengthy, and the 
results difficult to interpret. It is more instructive to 
consider the values of the integrals in the case of a long 
neutron range. This has some bearing on the physical 
situation, because the forces between the alpha-particle 
and the incoming neutron in a scattering experiment are 
thought to be relatively weak. 

In the “long neutron range approximation” ™ it is 
assumed that a>>@. Particles 1, 2, 3 (in the unsymmetrized 
wave function) are assumed to be concentrated in a tight 
group. In what follows the “spread” of 1, 2, 3, will be 
neglected. The angular factors in the integral are con- 
siderably simplified under these conditions, and the integral 
over 1, 2, 3 may be replaced by.an integral over 1 alone. 
One obtains a long formula for (AE)exchange which the 
author will be glad to furnish to those interested. 


Numerical results are given in Tables I-IV. 
The values of the constants used are, in nuclear 


units: 

Formula Constant Value Range in 107" cm 
J=B exp (—bdr’) b 20 2.01 =10*3/4/b 
u=exp 6.1 3.64 = 10" 


0.81 10.0 =10"/4/8 
B 1.65 7.0 
5.06 4.0 


The value of 6 chosen is in fair accord with the 
interpretation of experiments on the scattering 
of protons by protons. A variational calculation 
to be described later is the basis for the value of a. 
In an expression of the form exp [—cr*], the 
range is defined as 1/1/c. 

From Table I it is seen that the ‘long neutron 
range”’ [I.n.r.] approximation is reliable for the 
range of values for 8 considered. Table II gives 
the contribution to AE due to the exchange 
terms. Table III gives the contribution to AE 
from the ordinary terms, but differs from Table I 
in that exchange terms are included in the 
normalization integral now. Table IV gives the 
ratio of exchange terms to ordinary terms, with 


R(p) =exp (—8p?*) 


16 Referred to as I.n.r. approximation. 
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the normalization integral including exchange 
terms in both cases. It is seen from this table 
that the exchange terms give a very appreciable 
addition to the splitting on the alpha-particle 
model. 


Central Field Model 


Certain results obtained by Breit and Stehn in 
their treatment of the fine structure of the 
nuclear ground level of Li’ can be used in the 
discussion of He® on the central field model. The 
wave function is 
= (1/+/3)[1 — P35 | 

X ]RiR2So(12). (15) 
Here §£=x+7y, Q describes the particle outside 


the alpha-shell, and R a particle in the alpha- 
shell. Suppose 


J=Aexp[—ar’], 
R=WNreexp[ —ur?/2]. 


Introduce the abbreviations: 
One has 
AE=[(AE)r— (16) 


where the expressions on the right side are 
defined and evaluated in the work of Breit and 
Stehn. From this 


AE= (17) 
For n=», it follows that d=g, and 
AE= (18) 


In this case the exchange term is 0.20 of the 
total AE. 
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As a preliminary to the comparison of the fine 
structure as calculated on different models, it is 
useful to have relations between the Gaussian 
function (1234) in 
the interparticle distances, and the product 
function 


u(1234) =exp exp[ 
Xexp[—uras’/2] exp (19) 


where »/2 is written for comparison with the 
central field function 
although in the present case fq: differs from 7, 
the former being taken relative to the center of 
mass of the alpha-particle group: 


Tai =1—(1+2+3+4)/4, 
Ta2=2—(1+2+3+4)/4, (20) 


If we let a=y/8, B=(8v+2y)/25, or 
u=8a, v=258/8—2a, (21) 


the product function in the five particles is 
found equal to the expression in 734, p. 

Central field functions are referred to a fixed 
point—therefore the above relations are not 
strictly applicable. They will be used, however, 
in estimating the uv and » to be used for compari- 
son with a given a, 8. This amounts to neglecting 
the center of mass motion on the central field 
model. 

For a=p/8; B=2y/5, or B=16a/5. With 
a=6.1, B=19.52 (range ~2X10-'* cm). This 
value of 8 is unlikely physically, but it is 
apparent from Table IV that on the alpha- 
particle model the exchange terms are relatively 
more important than on the central field model. 

Using (21) one has for (AE)ordinary on the 
central field model, 


(AE) ordinary = — OX (25a8 — 16a*)/(25a8 
— 16a? +6(258/8+6a)) (22) 


TABLE I. Table of values of (AZ)ordinary/B. Only ordinary 
terms included in normalization integral. 


8B Exact L.n.r. approx. 
0.81 —9.48 x 10-5 —9.50 x 10 
1.65 —4.46x 1074 —4.50 1074 
5.06 —3.33 x 1073 —3.41 X10 
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Geometrical Conditions for Vanishing of Ex- 
change Terms on Alpha-Particle Model 


It can be shown that on the alpha-group 
model of He® non-vanishing exchange integrals 
exist even when the alpha-group does not overlap 
the outside neutron. It is in fact necessary that 
rm», the radius of the sphere about the center of 
mass of the alpha-particle from which the outside 
neutron is excluded, be greater than 57,/3, where 
r, is the radius of the alpha-particle, in order that 
a non-vanishing exchange integral exist. The 
reason for this behavior is that the center of 
mass of the original particles in the alpha-group 
is displaced when the outside neutron is ex- 
changed with a neutron in the alpha-group. 


III. DEUTERON GROUP MODEL OF THE 
ALPHA-PARTICLE 


Up to now the alpha-particle wave function 
has been used in the form 


u(1234) = exp[ —a(ris?+ +7347) J. 


It is of interest to consider the effect on the level 
splitting of a change on the form of u. The 
alpha-group can be thought of as composed of 
two deuteron groups; u will be represented as the 
product of factors representing the internal 
motion of each deuteron with a factor describing 
the relative motion of the center of masses of the 
deuterons. The complete alpha-particle wave 
function is assumed to be 


W4 = fexp[ tres") 
J 
— So(12).S0(34). (23) 

The relationship of this wave function to a 
deuteron grouping is seen from the relation 
exp[ — +123") 

+12.) ]=exp[ — airy” 


(113? +124") ], (24) 


TABLE II. Table of values of (AE)exchange/B. 


where rp=1+3—2—4, or twice the vector be- 
tween the centers of mass of the two deuteron 
groups. This follows from 


rp? =n (25) 


Two terms in the spatial part of ¥4 are necessary 
to make the space part symmetric in 1, 2 and 3, 4. 

The ranges a; and ay will be determined by 
minimizing the energy integral 


E=(W4, H¥4)/(¥4, (26) 
H=T+V, (27) 


Ore? rat Oras 


4 4 3? 

+ © (cos (ten) )} (28) 
l<n 


Ji; ; 
J=Aexp[—br*], (29) 


where P™” is the Majorana exchange operator, 
and P¥ the Heisenberg exchange operator. The 
effect of g is to decrease the attraction between 
like particles relative to that between unlike 
particles. 

The minimization of E is done by first mini- 
mizing the symmetric function with ar=an. 
With the value of a; found in this manner, ar is 
varied in the general expression for E, and the 
minimum determined graphically. It is too labo- 
rious to minimize E simultaneously with respect 
to both a; and ay. The results are given in 


TABLE III. Table of values of (AEZ)ordinary/B. Exchange 
terms included in normalization integral. 


B l.n.r. approx. 
0.81 —9.48x 107% 
1.65 —4.43 x 
5.06 —2.94x 


TABLE IV. Exchange terms included in 
normalization integral. 


exchange/‘-) ordinary} 


06 —2.48X 10-3 


8 l.n.r. approx. B l.n.r. approx. 
0.81 —4.50x 10 0.81 0.475 
1.65 —2.47X 10-4 1.65 0.557 
5. 5.06 0.843 
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Table V, in which we use as constants A = 85.0, 
b=20.0, g=0.22. Nuclear units are used. For 
a1=ay, the minimum occurs at ay=6.10, where 
E= —52.54 mc*. With this value of az, a graph of 
E against ay indicates a minimum at ay=6.5 
+0.05. The effect of g is to increase the internal 
binding of the deuteron subgroups. 

An alternate form of deuteron grouping was tried in 
which the spin terms are built up of deuteron spin func- 
tions: 


—exp 
(30) 


$1 = (1/ V3) (Si(13)S_1(24) — So(13)S0(24) 
+5S_1(13)5i(24)), (31) 

— So(14)So(23) 
+S_1(14)Si(23). (32) 


For a, =ay1, this function becomes identical with the one 
previously considered. The binding energy is slightly less 
with the present function. It can be shown that 4 
differs from 4 only by the addition of a small term with 
like particle spins coupled parallel. Since Heisenberg 
forces favor antiparallel coupling, ¥:4 is not as good a 
function as ¥4, The fine structure was not calculated for 
v4. 


Estimates of (AE) ordinary on the Different Models 


Constants used : b= 20, a=a;=6.10, ar = 6.60, 
B=75. Equation (22) used for central model. 
B=10 

(AE) aipha = — ] = —0.605 
(AE) centrat = — ] = —0.516 
B=19.52 (when 
(AE) alpha > (AE) central 
= — 6x°Bb(0.224] = —1.1 me’. 
Concentrated Alpha-Particle (a— = ) 
(AE) sipha= — 28) 
= for B=10 
= —6x?Bb[0.36] for B=20. 


In this limit the exchange terms vanish, and the 
above values give the total splitting. 


Central Model, letting ~ 


It is seen from these results that so far as the 
ordinary terms are concerned the ‘‘deuteron- 
grouping” alpha-function gives essentially the 
same AE as the simple alpha-function 


u=exp [—a(rie? +--+ +734) J. 


TABLE V. a1 =6.10. Energy in me’. 


ay, = 5.0 6.1 6.2 6.6 7.0 


Kinetic energy 103.1 109.80 110.32 112.631 114.961 
Potentialenergy -—155.07 —162.34 —162.86 -—165.244 —167.296 
Total energy —52.0 —52.54 -—52.54 —52.61 —52.34 


The results are not very different from those 
calculated on the central model with corre- 
sponding ranges. In comparing with the central 
model no correction has been made for the center 
of mass motion. 


IV. SPLITTING DUE TO INTERACTION OF A CON- 
CENTRATED s GROUP WITH AN 
ALPHA-GROUP 


It will be shown that the spin-orbit interaction of an s 
group, in an $1/2 state of internal motion, with an alpha- 
group is the same apart from a numerical factor as the 
interaction of a single particle with an alpha-group, in the 
limit of a very concentrated space function for the outside 
group. This calculation is of interest in connection with the 
theory of holes. The result will be given for the case of a 
Psi2 state of the complete system. 

Let N be the total number of particles, of which N—4 
are in a group in an S$1/2 state of internal motion, and 4 are 
in an alpha-group. The exchange terms vanish in the limit 
of a concentrated function for the outside group; thus 
ordinary terms only will be considered. The result will be 
given for any odd integer N in order to exhibit the nu- 
merical constant mentioned above; in practice only N=5 
and 7 are possible because of the exclusion principle. 

The wave function is 


X So(34).So(56) So( N—2, N—1)a(N); 
v(@) = (x, +typ)R(\e|), 
where 
(33) 
One finds 
(W, H’V) = (1234) N) 
X [2R*(p) +4/3pR'(p)R(p) dri...w. (34) 
To see the resemblance to the corresponding result for He’, 
one introduces a delta-function for w: 


6(N—[N—1])]!, (35) 
where 6(N—5) denotes 6(xy —x5)5(yw — ys)6(sy—25). Then 
the integral becomes 


S In u?(1234)[2R2(p) +4/3pR'(p)R(p) 
where @=N—(1+2+3+4)/4. Finally 


(AE) = — fu2(1234) Jy iL2R2(p) 
+4/3pR'(p)R(p) 
(LN —4 #2/3 f-u?(1234) p?R*(p)drizaw 
= [N—4] *(AE) ordinary, alpha-model, (36) 
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The factor [N—4]~* can be understood from a classical 
model, according to which 
— Fradial/2c*. (37) 


For a given angular momentum 7, is [V—4]" that for 
He’, and fradial is that for so that [V—4]? 
enters the result. 

The preceding results are also applicable to a model in 
which central field product functions are used for the 
alpha-particle, and also for the center of mass of the 
outside group. In this case one finds 


(AE) =[N—4]-*(AE) ordinary, central model, Hes. (38) 
V. FINE STRUCTURE OF Li’ 


The nucleus Li’ is pictured on the alpha- 
particle model as built up from an alpha-particle 
group and a triton (H*) group. The unsym- 
metrized wave function for the p3,2 state is taken 
to be of the form 


WV =u(1234)w(rs6, 767)0(@)So(12) 

X S(34)S0(56)a(7), (39) 
so that the antisymmetric wave function is 
w4=(1—Pi;— P2;)(1—P35— Ps6— Pas 

(40) 


y will be used to denote the space part of ¥. 
Particles 1, 2, 7 are protons; 3, 4, 5, 6 are: 
neutrons; u is a central field product type 
function, 


u(1234) =R(1)R(2) R(3) R(4), (41) 
v(o) = (x, +iy,)Q(p), (42) 


o= 3(5+6+7). 


where - 


This is a central field representation of v(). w is 
an s function, and will be assumed symmetric 
in 5, 6, 7. 


(W4, H’W4) =18(, (43) 
(WY, (44) 


where J) contains the terms classified as ordinary, 
I, the terms with a single permutation operator, 
I, the terms with two permutation operators, and 
]; those with three permutation operators. 

The ordinary terms can be evaluated from the 
equation for central field functions analogous to 


(34), with N=7. Setting 


J=B exp [—br*], 
R=exp [—ar’], (45) 
Q=exp [—8p*], 


w=exp ], 
one has 


(AE) ordinary = — (2X?/3) (2aBy)/(2aBy 
(46) 


Consider various limiting cases of this result: 
For y—* (concentrated alpha-group, 
concentrated triton) 


(AE) ordinary = — }*?. (47) 
For a— (concentrated alpha-group) 


(AE) ordinary =< (2x?/3)Bb 
(48) 


For (concentrated triton) 


(AE) ordinary == 
X[(2a8)/(2a8+6(8+a)) 
= 1/ 9(AE) central mode! He5, ordinary terms. (49) 


Discussion of Exchange Terms for Li’, Alpha- 
Triton Model 


It will be shown below that the exchange inte- 
grals in the limit of small overlapping of the 
group wave functions have a value proportional 
to the degree of overlapping. It is of great 
assistance to be able to eliminate double and 
triple exchange integrals from the calculations. 
This may be done by the following reasoning in 
the case of a small overlap. On a central field 
model of Li? one has the following sequence, 
where v is a central field function for a particle in 
the triton: 


Ordinary: 
Single exchange: 
Double exchange: 
S RER Sve d tise S 
XS 2". 
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Triple exchange: 
S Ryv3Rsvsd 735 
X S 727 = 


If x is small, then x? and x* may be neglected. On 
going over now to the case of a group model for 
the triton, it will be seen that the corresponding 
double and triple exchange terms are smaller 
than x? and x’, respectively. Suppose that the 
contribution x to the single exchange integral 
comes from a small region of space. It may not be 
possible for more than one particle in the triton 
to occupy this region simultaneously, because the 
center of mass of the triton is restricted in its 
motion to remain, say, outside of a certain 
spherical region. In this example the double and 
triple exchange terms would vanish. These terms 
may become as large as x? and x* only if the 
triton is completely included in the sphere of the 
alpha-particle. 

The general plan of this section is to work out 
the single exchange terms in a form depending on 
the overlapping, and certain vector and scalar 
products. These products will first be estimated 
roughly in the limit of a very small overlap, and 
the values obtained in this manner will be the 
basis for a classification of the terms as large or 
small. The small terms will then be neglected. 


The terms which are not neglected will be esti- — 


mated on a more refined model. Henceforth one 


assumes 
(Vv, (50) 


Division of Terms into Large and Small Parts 


The division proceeds by calculating the aver- 
ages of various vector and scalar products ap- 
pearing in J;, such as [rs }? and (r;-o), in the 
limiting case of small overlap for very simple 
group wave functions. 

The result of the calculation of the averages is 
that 


I, = — (7/108) CS PR sw (567) 
Jdri---7, (51) 


where several small terms contributing about 1 
percent of the result have been neglected. This is 
free from any direct assumption about the wave 
functions. It need be supposed only that the 


“small’”’ terms are always negligible and that 
J'/r=C. 


Measure of Overlap 


The term P3sy) in the normalization inte- 
gral will be used as a measure of the overlapping. 


It is supposed that (y, y) =1. 
(v, Paw) = S sw(567)w(367) 
X QQ 5) (X_ dt 
= (2/3) S sw(567) 
Xw(367)0,0,p%dr, (52) 


for small overlapping, in which case the whole 
contribution to the integral comes from 9=n. 
Define 


Q = (8d?/3) S 
Xw(367)0,0.d7, (53) 


which is the value of the above integral for triton 


radius=alpha-radius=d in the limit of small | 


overlap. It is, however, approximately valid as a 
measure of overlap for finite overlap—the error 
is that of replacing (9-n) in the integral by the 


constant 
On this model, 


I,= —7x7CQ/16, (54) 
and to the same approximation 
¥4)=(Y, ¥) —3(Y, Posy) 
(55) 
when we take (¥,¥)=1. Thus the complete 
splitting is 
(AE) = — (x?/2) MMs)? 
]w*(S67)dr 
+7€2/16]/T SR: 
XQ,w(567))*dr—32]. (56) 


It is instructive to evaluate (AE) wholly on the 
model used for the exchange term: 


(AE) = —X*[—Cly, 
+7€2/16]/L(y, —32] 
=Cx2(1+2.62). (57) 


The exchange term has the same sign as the 


hat 


'te- 
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ordinary term, in agreement with the central 
model result. The doublet is inverted, since C is 
negative. 

The overlapping can be calculated directly, 
with Gaussian functions and a product type 
central field representation of the alpha-group 
as in (46): 


2=[12a8/(3a*+ 68y+6ay+ 1408/3) (58) 


where several very small terms have been omitted. 

The range of the wave function 9% should be 
corrected for the center of mass motion. This can 
be estimated from the result for the mean square 
displacement of the center of mass of the alpha- 
particle: 


X 
[16S 
=3a/16, (59) 
where ¢ is the coordinate of the c.m. of the alpha- 
particle group. Suppose now a@ on an independent 
particle model is given, and one wishes to find the 


corrected range. The apparent range is 1/\/a=f.; 
then the corrected range is 


ta— = (3/16) (1/ax)! 
=(1—v3/4)ra. 


The overlap integral (y, Pssy) = can be ap- 
proximately related to the periods of motion of 


the particles in the alpha-group and the center of , 


mass of the triton, considering the triton group 
to be very concentrated. By utilizing the uncer- 
tainty principle, we have 


Q=range particles in alpha-group/ 
range c.m. of triton 
=momentum c.m. of triton/ 
momentum particles in alpha-group 
=v period particles in alpha-group/ 
V period c.m. of triton. (60) 


If the alpha-group is considered to be very 
concentrated, 


period c.m. of alpha-group/ 
/ period particles in triton. (61) 


Similar estimates can be made by using the 
W. K. B. approximation as a solution of the 
radial wave equations for the various particles 
and groups, but no exact connection can be made 
between the overlap and the periods of motion, 
because the spatial characteristics of wave func- 
tions such as the presence of ‘“‘tails” affect the 
overlap integrals. 

The splitting for Li’ on a true alpha-particle 
model in which only relative coordinates appear 
has not been discussed because of the great 
amount of labor involved in treating such a 
model. An investigation of the geometrical con- 
ditions for the vanishing of the exchange terms 
has been carried out, however, in a manner 
similar to that employed on the same question 
for He®. It turns out that for non-overlapping 
groups the single and double exchange integrals 
I, and J, vanish, but the triple exchange integral 
I; does not vanish until the distance between the 
c.m.’s of the two groups is at least 7/3 the radius 
of the alpha-group. 


VI. CONCLUSIONS 


1. The contribution of the exchange terms 
relative to ordinary to the splitting in He® is. 
greater by a factor of the order of two or three on 
the alpha-model than on the central field model. 
The sign and magnitude of the total splitting are 
of about the same order on the two models. 

2. The ‘‘deuteron-grouping’’ alpha-function 
gives for He’ a splitting within a few percent of 
that for the ordinary alpha-function. 

3. The exchange terms for Li’ are approxi- 
mately proportional to the degree of overlapping 
of the alpha- and triton groups. The exchange 
splitting has the same sign as the ordinary 
splitting. 

4. Energy splitting estimates on the alpha- 
model in which the exchange terms are neglected 
should give results reliable to an order of magni- 
tude at least so far as this model is concerned. On 
reasonable assumptions the value of the exchange 
terms is of the order of one-half of that of the 
ordinary terms. 

I wish to acknowledge my great indebtedness 
to Professor G. Breit for many discussions and 
suggestions on this work. 
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Electromagnetic Properties of Nuclei in the Meson Theory 
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(Received May 22, 1942) 


The mathematical treatment of the interaction of a meson field with an electromagnetic 
field is simplified by using Pauli’s method. A discussion of calculation of exchange multipole 
moments of nuclei by the method of canonical transformation and the results for dipole mo- 


ments in the pseudoscalar theory are given. 


1. INTRODUCTION 


HE theory of the interaction of a meson 
field and an electromagnetic field has been 
developed by several authors.! More details have 
been worked out by Bhabha. The latter has 
given complete expressions for the electro- 
magnetic interaction of a meson field in the 
presence of nuclear particles and also equations 
of continuity for the meson current and the 
current of nuclear particles separately. Both the 
mathematical treatment and the final result can 
be simplified if we adopt a method used by 
Pauli.? Pauli’s treatment is based on the general 
properties of the Lagrangian function instead of 
the explicit form of Proca’s equations which are 
derivable from the Lagrangian function as 
variational equations. The method is therefore 
quite general and is applicable not only to the 
vector meson theory, but also to the scalar, 
pseudovector, and pseudoscalar theories as well. 
Pauli treated only the electromagnetic inter- 
action of a free meson. We give below a general- 


ization of his method to include nuclear particles. - 


Owing to the exchange nature of the meson 
theory, it becomes now necessary, in treating 
the interaction of a nucleus with an external 
electromagnetic field, to take into account not 
only the direct influence of the electromagnetic 
field on the nuclear particles, but also the 
interaction of the external field with the exchange 
meson field. Investigations in this direction have 
recently been made by several authors.’ It has 


1 Yukawa, Sakata, Kobayasi, and Taketani, Proc. Phys. 
Math. Soc. on 20, 319, 720 (1938); Fréhlich, Heitler, 
and oe, = roc. Roy. Soc. A166, 154 (1938); Kemmer, 
Proc. Roy. Soc. A166, 27 (1938) : Bhabha, Proc. Roy. 
Soc. A166, 501 (1938). 

Post lg of the Solvay Conference 

a Phys. Rev. 52, 787 (1937); 
I. Schiff, Phys. Rev. 53, 651 (1938); 
Frohlich, Heitler, and Hahn, Proc. Roy. Soc. Al74, 85 


been shown that the electromagnetic interaction 
of a nuclear system can be expressed in terms of 
interactions between the external field and the 
electric and magnetic multipole moments of the 
nuclear system. The electric and magnetic dipole 
moments and the electric quadripole moment 
have been calculated on the vector theory by 
means of the ordinary perturbation method in 
quantum theory. It might appear objectionable 
to use a perturbation method for this purpose, 
because perturbations of as high as the third 
order are involved in such calculations. It can 
be shown, however, that all the previous results 
can be obtained, in a simpler manner, as the 
first approximation in the canonical transforma- 
tion carried out by Mller and Rosenfeld.4 We 
give below a discussion of the calculation of the 
exchange multipole moments by the method of 
canonical transformation and also give results 
in the pseudoscalar theory. 


2. GENERAL THEORY OF THE INTERACTION OF 
MESON FIELD WITH ELECTROMAGNETIC FIELD 


We give below a mathematical treatment of 
the interaction of a meson field with an electro- 
magnetic field in the vector meson theory on 
the basis of Pauli’s method. Pauli’s method may 
be easily extended to other forms of meson 
theory. We give at the end of this section the 
results for the pseudoscalar theory, which can 
be similarly derived. 

The vector meson field is described by a four 
vector U,, with Us=e€Uo, € being the square root 
of —1 introduced by Proca.®’ Let the complex 
conjugates of the field variables U, be denoted 


Hong S. T. Ma, Proc. Camb. Phil. Soc. 36, 351, 438 
1 

‘Moller and Rosenfeld, K. Danske Vidensk. Selskab. 
17, 8 (1940). 

5 Proca, J. de phys. et rad. 7, 347 (1936). 
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by U,*. We introduce the operators 
n,=V,—1%,, 2,*=V,+i®,, (1) 


where V,=0/dx,, Po, #3) = (e/hc)(A ly A», 
and )=@,/e=(e/hc)Ao, A and Ag being the 
electromagnetic potentials from which the anti- 
symmetric tensor f,, representing the field 
strengths can be derived. We also adopt the 
notations 


From these we define the new field variables 


(2) 


Gp = Us, — Ups, Us» = (3) 
It follows from the relations 
(4) 
hc hc 
that 
te 
Us» Upr» = es 
(5) 


* * 


te 
hc 


Furthermore, for two complex scalar quantities 
A and B, we have 


(6) 


Let the wave functions of the nuclear particles 
be denoted by ¥ which consists of the relativistic 
wave functions of the proton and the neutron. 
On introducing the notations 


Wi=eW*y', y’=—iBa’, at=e, (7) 
where a’, a, a, 8 are Dirac’s matrices, we obtain 
the following four vectors and six vectors: 

>t Vn = Ir, = re. 

Vy 
The complete Lagrangian function for the 
system of nuclear particles, the meson field, 
their mutual interaction, and their interaction 


with the external electromagnetic field may be 
written as follows: 


(9) 
L?= pt (chy’x,+M Vp, 
Ly=- (chy?V,+ Myc?) Vy, 


(8) 


(10) 


MESON THEORY 119 


(11) 


LPN=— 
(12) 


LY? = —g 
kK 


It follows from the variational principle 


if [ f (13) 
AL 


aL aL 
aU, au,’ “aus 


that 


OL OL 


(15) 


AL OAL 
Vy ’ ° ( 1 6) 
OVy, 


— 
OVy 


From these equations we obtain 


hv 

v.( » OL & 


(17) 


}) = + 


The Lagrangian function defined above, be- 
sides being relativistically invariant, is invariant 
under a gauge transformation. A gauge transfor- 
mation for the wave functions of the protons, 
neutrons, and mesons is 
Vp—WVpe*, 
VyVy, (18) 
U,-U,e*, U,*-U,*e-*, 

where a is an arbitrary function of x1, x2, Xs, 


and x4. The corresponding gauge transformation 
for the electromagnetic fields is 


Via. (19) 
From (18) and (19) it follows that 
(20) 
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and therefore and also that 


(22) 


* * 
U, U, U, py U, pane 
(21) Since L is invariant under a gauge transforma- 
tion, we have 


Vv Vn», aL/aa=0. 


Vp, > 


Putting a=0 in this equation we obtain 


With the help of this equation we obtain 


The energy _ momentum tensor may be defined as 


aU* P.N 
By (17), 
aU* dU,* oVv,* ay, ov* P.N 


Substituting (24) into this expression we obtain the result 


where 
1e OL OL OL ) (28) 


may be interpreted as the four vector expressing the current and charge densities of the nuclear 
system including the nuclear particles and the meson field. From (28) we have 


. (29) 


The equation of continuity for the total current 


follows from (29) and (23). 
The current s, can be separated into two parts: the proton and the meson current defined, re- 


spectively, by 


) 
—— 


I 
V,s,=0 (30) 
te 


2) 


3) 


4) 


5) 


6) 


7) 


8) 


ar 


9) 


1) 


2) 
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The two parts do not satisfy (30) separately because of the exchange of charge. We have 


he ov" OVp OV py, 
On account of the relations 
OV p* ov" OV OVp 
we obtain the result ; 
= = (33) 


This equation expresses the rate of flow of the exchange current into or out of the nuclear particles. 
It has been given by Bhabha! but it is here given in a much simpler form. 
The explicit expressions for the variational equations and the current densities are 


(34) 
ie 
where 


Let us now write the four and six vectors in terms of ordinary vectors and scalars as follows: 


U,=(U, V), G»=(G,F), ete.), 
gil’?=(M, giVp* Vy), (37) 


=(S, T) = “Sy 
k kK 
Equations (34) and (36) become then 
10 
Cc 


(38) 
G=(vV —i®) XU+S*, 


1a 
c at 
together with the complex conjugates of these equations. The meson current and charge densities 
are, in the present notations, 
ie 
j=c(S1, Se, $3) = V*F— VF*), (39) 
1e 


In order to link up the present notations with that of reference 3, we divide U and F into longi- 
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tudinal and transverse parts and put 


1 K 


The longitudinal part of U and the transverse part of F may then be expressed in terms of @ and ¢ 
with the help of (38) with &o, ®, M, N, S, and T omitted. We have 


The Hamiltonian is equal to the 44 component of the energy-momentum tensor 
Typ = — + 7p) V (42) 


multiplied by —1. Since 


and 
= + U,*(x2U, — J 


], 
we may write 


T te) «U,*U,—gil’U, the complex , 
conjugate. (43) 


Hence 
H= V 


kK 


3 


i, j=1 


+5 F— N*|?+ |F—T* |?—M-U+3 > S4G;;+the complex conjugate. (44) 


i, j=1 


This result can be separated into a term independent of the electromagnetic field and a term de- 
scribing the electromagnetic interaction, namely, 


Hy= 
1 1 
xU)— 
K 


+the complex conjugate+terms containing go”, (45) 
and 
1 
H,=-—-j-A. (46) 


In the second equation terms proportional to the square of the A’s have been omitted. In order 
to obtain the interaction with an external electrostatic field Ao we have further to introduce into 


the Hamiltonian the zero term 
1 
(47) 


where p is the charge density given by (40) and E is the electric field strength.® 
6 Heisenberg and W. Pauli, Zeits. f. Physik 59, 168 (1930). 
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Equation (44) can also be reduced to the form 


= chy: (V —i@rp)+(Mpre+ Myrty)c? . 
(48) 


if we add the relativistically invariant term 


1 22" 
=§$*-S—T*-T 
2 «* 

to the Lagrangian function (9). 

The treatment of other forms of meson theory is similar to the above. We state briefly the results 
in the pseudoscalar theory. Let y be a pseudoscalar describing the meson field and y,=2, be the 
pseudovector derived from y. For the nuclear particles we have the following pseudovector and 
pseudoscalar : 

TA pty yy Vy, Te Vp, 
(r24*, [321*) = re. res. re"). 


pty =I", = — ie Wyt Wp =I". 


(49) 


While the Lagrangian function for the nuclear particles is still given by (10), for the meson field 
and its interactions we have 


LP%= —fil’y —fir'*y* (50) 
The variational equations for the pseudoscalar meson field are 
(51) 
or 
TpPp= — Ky fil”, (52) 
where f f 
= =-y, (53) 


The present scheme is also gauge invariant and therefore (23) still holds. From the energy-momentum 
tensor 


aL 
T= 


—Lé,, (54) 
P.N 


aw, ot Y, 
it follows that the four vector of the current and charge density is given by 


$,=5,?+s,™, 
where 


(55) 
he 


The proof of the equations of continuity (30) and (33) is the same as in the vector theory. 
Let us rewrite the above four-dimensional pseudovectors and pseudoscalars in terms of three- 
dimensional pseudovectors and pseudoscalars as follows: 


= (Tr, (56) 


fil”=R, r’, =P, (57) 


- | 
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Then (52) and (53) become 
1a 
(- < +i) (- < —i®)-T—R*, =—(V—i®)y+P*, (58) 
c c 


together with the complex conjugates of these equations, and (55) takes the form 


= (59) 
p™ = (ie/hc)(y*e— (60) 
The Hamiltonian is given by 
= —1Vt {chy —i@rp)+(Mpret+ Myry)e?} V+ yt | ay o—fil’y 
3 2 
K K 


k=1 


in which the term describing the electromagnetic interaction is of the form (46), j” being now given 
by (59). Equation (61) may be rewritten in the form 


H= —iWvt {chy-(¥ —i@rp)+(Mpret+ Myr )c?} V+ — Ry 
(62) 


which is the expression given in reference 4. 


3. CALCULATION OF EXCHANGE MULTIPOLE MOMENTS 


As in the case of the nuclear forces the multipole moments of a nuclear system due to exchange 
processes may be calculated either by using the ordinary method of quantum theory or by adopting 
the method of canonical transformation in reference 4. Calculations of the multipole moments by 
the former method have been given in reference 3. We discuss below the method of canonical trans- 
formation and also calculate the multipole moments in the pseudoscalar theory. 

Let us take as an example the exchange magnetic moment of the deuteron in the vector theory. 


The static parts of the field variables are given by 
= 


U%(r') = f XS*(r) o(r—r’)dr, 


(63) 
F°(r’) = — f VN*(r)o(r—r’)dr, 
| together with the complex conjugates of these equations, in which 
1 exp (—«|r—r’}) (64) 


in 


58) 
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62) 
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The exchange magnetic moment is given by 
Take any point O as the origin. Let the vectors OP and ON from O to the proton and the neutron, 


respectively, be denoted by R, and Ro. We write R,—R.=R, R:=0, r—R,=n, 
Then the first term in the integrand of (65) gives a contribution 


f f f 


With the help of the relation 
1 
f (66) 
82x 


we obtain 


1 1 
167x 167 R 


The first two terms in the integrand of (65) give therefore a contribution 


f f dR. f f (: XR ) Nit (67) 


16xch 2ch 


It is known that (see reference 3) the total exchange current between the proton and the neutron 
is equal to (ie/h)J, J being the nuclear potential between the proton and the neutron due to exchange 
processes. Equation (67) is therefore just one part of the magnetic moment due to the exchange 
current with respect to the point O. If we choose the origin at the midpoint of the two particles 
then 9=0 and (67) vanishes. Henceforth we shall omit all terms containing @ in the final result. 

The contribution of the third and fourth terms in the integrand of (65) can be similarly calculated 
with the help of the general formula for partial integration 


fov-cexwar= f(xw)-vudr+ [uw-(w xv)dr. (68) 
The final result for the exchange magnetic moment with respect to the midpoint of P and N is 
tegs" 
167, 


in agreement with the result obtained in reference 3. 
In the pseudoscalar theory the exchange magnetic dipole moment is given by 


le 
u= rxX 


P.)(¥1- Py*) — (x2Py* +1 X01 X P,*) Ps) o(r2)drdRidRo. 


On evaluating the integral as above, we find 


nf ‘RR = (1+) fe (70) 
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The result is antisymmetric in the spins of the two nuclear particles. A deuteron in any stationary 
state has therefore no extra non-additive magnetic moment due to meson exchange in the pseudo- 


scalar theory as well as the vector theory. 


In contrast to the result in the vector theory, in the pseudoscalar theory a proton-neutron pair 
has no exchange electric dipole moment that reacts with an electrostatic field. There is also no 


electric quadripole moment that reacts with an electrostatic field. 
In conclusion the writers would like to thank Professor W. Heitler for a comment at the early 


stage of this work. 
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It has been suggested that the interaction between a pair of nuclear particles contains, in 
addition to a central force term, a tensor term, and quantitative estimates of the relative 
magnitudes of these two terms have been given. In this paper a rough quantitative estimate 
of the relative contributions of these two terms to the binding energy of a nucleus is made on 
the basis of a highly simplified nuclear model. It is found that a tensor interaction between two 
nuclear particles of the form suggested by the neutral meson field theory of nuclear forces with 
constants determined to fit the experimental data on the two-particle neutron-proton system 
would lead to the existence of moderately light nuclei having binding energies, spins, and 
deviations from spherical shape much greater than those observed. 


1. INTRODUCTION 


HE discovery of the electric quadrupole 

moment of the deuteron! suggests the 
existence of a term lacking spherical symmetry 
in the interaction between a proton and a 
neutron. The current assumption of equality of 
forces between all pairs of nuclear particles then 
implies the existence of a similar term in the 
interaction between two like particles. General 
considerations of invariance? indicate that such 
a term should contain the operator 


— 01-02. (1) 


Field theories of nuclear forces have been 
proposed? which lend support to this view. They 


* The greater part of this work was completed during 
the spring term of 1940 while the author was in the 
Department of Physics at Princeton University, Princeton, 
New Jersey, on a Research Fellowship granted by the 
Royal Society of Canada. 

c M. B. Kellogg, I. I. Rabi, N. F. Ramsey, Jr., and 
— Phys. Rev. 55, 318 (1939) and 57, 677 

1 


2E. Wigner, Phys. Rev. 51, 106 (1937). 
3H. A. Bethe, Phys. Rev. 57, 260 (1940) (gives a list 
of references); R. E. Marshak, Phys. Rev. 57, 1101 (1940). 


predict the existence of terms containing Sj. as 
a factor, but little reliance can be placed on the 
detailed form of the interaction given by these 
theories because of divergence difficulties from 
which they suffer. 

The influence of the S,. tensor term on the 
two-particle proton-neutron system has been 
investigated by Rarita and Schwinger* and by 
Bethe.* The former develop a purely phenomeno- 
logical theory adopting simplified rectangular 
well potentials with constants chosen to fit the 
binding energy and quadrupole moment of the 
deuteron, and the scattering of slow neutrons 
in hydrogen. The range of the forces is taken 
equal to that deduced from proton-proton 
scattering. Bethe employs the form of the 
interaction predicted by the single-meson theory 
of nuclear forces cutting off at small distances to 
avoid divergence. The value of the “cut-off” 
distance and the strength of the interaction is 
determined to fit the deuteron binding energy 


*W. Rarita and J. Schwinger, Phys. Rev. 59, 436 (1941) 
and Phys. Rev. 59, 556 (1941). 
5H. A. Bethe, Phys. Rev. 57, 390 (1940). 


lary 
ido- 


pair 
no 


arly 


62 


TENSOR FORCES AND HEAVY NUCLEI. 127 


and the scattering cross section of slow neutrons 
in hydrogen. A discussion of H* and He‘ based 
on the use of tensor forces is given by Gerjuoy 
and Schwinger.® 

It seems of interest to investigate the effect of 
the tensor interaction on the binding energy, 
spin, and shape of heavy nuclei. Wigner and 
Eisenbud? pointed out that forces of the general 
type of those used by Rarita and Schwinger and 
by Bethe will not have the property of saturation. 
Bethe* also mentioned the possibility of the 
tensor part of the interaction becoming pre- 
dominant in heavy nuclei with the consequences 
not only of non-saturation, but also of large 
nuclear spin and markedly non-spherical shape 
of the nucleus. The object of the present paper 
is twofold : (a) On the basis of a highly simplified 
nuclear model, we obtain a rough estimate of 
the relative contributions of the spherically 
symmetrical and of the tensor interaction to the 
nuclear binding energy. (b) The application of 
these results to the “neutral” theory type 
interactions using the numerical constants as 
determined by Bethe® shows that the difficulties 
mentioned above will become pronounced well 
within the range of existing nuclei. The same 
holds for the values of the constants as de- 
termined by Schwinger and Rarita‘ if the 
spherically symmetrical part of the interaction 
is due to an ordinary potential. The case that 
it is an exchange interaction is not considered 
in the present paper and may well lead to a 
different conclusion. This makes it very probable 
that in the Schwinger and Rarita theory, the 
Jo(r) of (2) has to be replaced by an exchange 
interaction. 


2. THE FORM OF THE NUCLEAR INTERACTION 


The general form of the interaction between 
two nuclear particles predicted by the neutral 
meson theory is given by 


(2) 


The symmetrical theory introduces an additional 
factor t,-t2, while the charged meson theory 
gives no interaction between like particles, and 


(1943) Gerjuoy and J. Schwinger, Phys. Rev. 61, 138 
1942). 
7 E. Wigner and L. Eisenbud, Phys. Rev. 56, 214 (1939). 


Sin 


=~ Gr 


Fic. 1. The factor giving the dependence of the nuclear 
interaction on the distance (in units of 1/x=2.1810-" 
cm) between the particles. F(r) and G(r) are Bethe'’s 
results on the neutral straight cut-off theory for the 
spherically symmetric and the tensor parts of the inter- 
action. S(r) is the step function used to replace both in 
the present calculations. 


introduces a factor (1+ :-*2) in the interaction 
between a proton and a neutron. 

Bethe® arbitrarily assumed the spin inde- 
pendent term Jo to be entirely absent, and took 
the other two terms to be of the form predicted 
by the meson theory of nuclear forces: 


Jo(r) =0, ’ 
3 


where 1/xk=h/yc=2.18X10-'* cm on the as- 
sumption that ».=177 m, and f? is a constant. 
To avoid divergence these interactions must be 
cut off at some value of r=a. Using straight 
cut-off (J(r)=J(a) for r<a) and the interaction 
(2) of the neutral theory Bethe obtained 


f?=0.0770hc, «a=0.405. (4) 
J, and Jz may then be rewritten as 
Ji(r)=AF(r), J2(r)=BG(r), (S) 


where F(r) and G(r) are unity for r=a, and for 
r>a 


F(r)=Ji(r)/Ji(a), G(r) =Ja(r)/J2(@). (6) 
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TABLE I. Values of the function W;(k, 6). 


k 5=0 0.5 0.6 0.7 0.8 0.9 1.0 1.5 2.0 3.0 
0.1 1.000 0.968 0.801 0.624 0.342 0.207 
0.2 1.000 1.000 1.000 0.946 0.805 0.483 0.300 
0.4 1.000 1.000 1.000 1.000 0.999 0.946 0.641 0.404 0.172 
0.6 1.000 1.000 1.000 1.000 1.000 0.988 0.723 0.448 0.183 
0.8 1.000 1.000 1.000 1.000 1.000 1.000 0.999 0.759 0.465 0.187 
1.0 1.000 1.000 1.000 1.000 1.000 1.000 1.000 0.768 0.469 0.188 
1.5 1.000 1.000 1.000 1.000 0.997 0.985 0.741 0.457 0.185 
2.0 1.000 1.000 1.000 1.000 0.975 0.943 0.692 ‘| 0.436 0.180 
3.0 1.000 0.998 0.975 0.893 0.843 0.598 0.394 0.170 
4.0 1.000 0.994 0.923 0.810 0.755 0.524 0.357 0.161 
8.0 1.000 0.879 0.717 0.588 0.537 0.361 0.259 0.134 
20.0 1.000 0.606 0.437 0.346 0.311 0.208 0.152 0.091 
TABLE II. Values of the function W2(k, 4). 
k 6=0 0.5 0.6 0.7 0.8 0.9 1.0 1.5 2.0 3.0 
0.1 —0.939 —0.908 —0.742 —0.566 -—0.290 —0.162 
0.2 —0.833 —0.833 —0.833 —0.779 —0.641 —0.336 —0.177 
0.4 —0.595 —0.595 —0.595 —0.595 —0.594 —0.543 -—0.283 —0.121 —0.029 
0.6 —0.370 -—0.370 -—0.370 -—0.370 —0.370 —0.359 -—0.175 -—0.069 —0.016 
0.8 —0.172 -—0.172 -0.172 -0.172 -0.172 -0.171 -0.076 -—0.029 —0.006 
1.0 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 
1.5 0.338 0.338 0.338 0.338 0.333 0.310 0.124 0.046 0.010 
2.0 0.582 0.582 0.582 0.582 0.535 0.483 0.194 0.073 0.016 
3.0 0.906 0.906 0.858 0.708 0.621 0.263 0.103 0.023 
4.0 1.111 1.101 0.962 0.754 0.654 0.290 0.119 0.028 
8.0 1.496 1.252 0.937 0.693 0.595 0.290 0.140 0.036 
20.0 1.778 0.953 0.662 0.484 0.419 0.223 0.128 0.043 


A and B are constants given by 
A=J,(a)=7.7 Mev, B=J2(a) = —102 Mev. (7) 


F(r) and G(r) are plotted in Fig. 1 together with 
the step function S(r) of range ro defined by 


S(r)=1 
=0 


Bethe also found that the factor 1-2 introduced 
by the symetmric theory alters the values of a 
and f? and leads to several serious difficulties. 
Rarita and Schwinger* assumed rectangular 
well potentials of equal range, but different 
depths: 
Jo(r) = — (1—4g) VoS(r) = —13.39S(r), 
Ji(r) = — 3g VoS(r) =— 0.50S(r), 
J2(r) = —yVoS(r) = —10.75S(r). 
S(r) is the step function defined by (8), and the 
constants g, y, Vo, and 79 were evaluated to be: 


g=0.0715, y=0.775, Vo=13.89 Mev, (10) 
= 2.80 X cm. 


(8) 


(9) 


The known properties of the ground state of 


the deuteron and low energy collision phenomena 
are insufficient to discriminate on the basis of a 
purely phenomenological theory between the 
interaction (2) and those obtained from it 
by the inclusion of the factors —4:-t2 or 
—43(1+ 1-2), which both have the value unity 
for states antisymmetrical in the charge. 


3. THE SIMPLIFIED NUCLEAR MODEL 


The following model is used in the present 
calculations: 

(a) The interaction V between two particles 
is assumed to be of the form (2) corresponding 
to the neutral meson theory. The particles are 
taken to have all their spins parallel. V then 
becomes an ordinary potential depending on the 
distance r between the two particles, and the 
angle @ which r makes with the direction of 
their spins: 


V=[Jo(r) +Ji(r) ]+(3 cos? @—1)Je(r). (11) 


(b) The functions J(r) are replaced by rec- 
tangular potential wells of appropriate range fo. 
The ranges need not all be the same, but will 
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be taken the same for convenience. Then: 
V=AS(r)+BS(r)(3 cos? @—1). (12) 


(c) The nucleus is taken to have the shape of 
an ellipsoid of revolution with its axis of sym- 
metry parallel to the spins of the particles. 
The shape of the ellipsoid is specified by the 
ratio k of its length along the axis of symmetry 
to the diameter of its greatest circular cross 
section (k>1 corresponds to prolate, and k<1 
to oblate ellipsoids). The volume 2 of the 
ellipsoid is specified by 6 which is defined by 
Q= (2/6) (5ro)*: i.e., 6 is the ratio of the diameter 
of the sphere having the same volume as the 
ellipsoid to the range ro of the step functions 
used. 

(d) The density p of the particles is assumed 
to be uniform over the volume of the nucleus. 

(e) The Coulomb energy, and the exchange 
integrals of the purely nuclear interaction are 
neglected, the total energy being calculated as 
the sum of the direct nuclear interaction and 
the kinetic energy of all the particles. 

If the Coulomb and the exchange integrals 
are not neglected, and if the surface effect® is 
properly taken into account then calculations 


based on the above crude model must necessarily 


give an overestimate of the total energy and 
hence an underestimate of the binding energy, 
so that any conclusions as to lack of saturation 
for such a model would a fortiori hold for a 
better one. The additional approximation made 
in neglecting the Coulomb, the exchange, and 
the surface energy makes this argument some- 
what less conclusive, but the results obtained 
below are still suggestive enough to make such 
a model of interest. 


4. EVALUATION OF THE POTENTIAL ENERGY 


On the above model the potential energy is 
given by 


P.E.=4 f f p(x) V(x—x’)o(x')dxdx 


=}N*[AW(k, 6) + BW2(k, 5) ]= W(k, 5) 


Here JN is the number of particles in the nucleus 


8H. A. Bethe and R. F. Bacher, Rev. Mod. Phys. 8, 82 


(1936); E. Feenberg, Phys. Rev. 60, 204 (1941). 


and 
Wilk, 3) =— f f S(\x—x' | )dxdz’, 

(14) 
Wilk, 3) =— f f S(|x—x'|)(3 cos? @—1)dxdz’, 


both integrations being extended twice over the 
volume of the nucleus. The integrals in (14) 
may be carried out analytically, and the results 
are stated in the appendix. From these results, 
Tables I and II, giving numerical values of 
Wi(k, 6) and W2(k, 6) as functions of k and 6 
for several selected values of these parameters, 
were constructed. 


5. KINETIC ENERGY 


In the statistical model the wave function of 
each individual particle is taken to be that of a 
free particle in a box of a size and shape corre- 
sponding to the particular nuclear model under 
investigation. The number of particles in each 
level is restricted by the Pauli exclusion principle. 
Usually two particles are allowed in each level 
corresponding to the two possible spin orienta- 
tions. In our case, where the spins are all taken 


Fic. 2. The total kinetic energy (in units of Ko= 
mh? /2M%) of N particles obeying the exclusion principle 
confined to a box of volume &. I—Approximate calcula- 
tion. Il—Cubical box, periodic boundary conditions. 
II1I—Cubical box, zero boundary conditions. [V—Rec- 
tangular box, ratio of sides 1 : 1 : 3, zero boundary condi- 
tions. V—Spherical box, zero boundary conditions. 


\ 

3.0 
).172 
).183 
).187 
).188 
1.180 
1.170 
1.134 
1.029 
.006 
.010 
.016 
.023 
.028 
.036 
.043 
ena 
fa 
the 

it 

or 
| 
ity 

L 

| 
11) 
ec- 
To. 
ill 


130 


TABLE III. Values of the function 


Walk, 8) =4[7.7Wilk, 8) —102W2(k, 8)J 
in Mev. 


k | 8=0.50 0.70 0.90 1.0 1.5 20 30 
1.0 3.8 3.8 3.8 3.8 3.0 1.8 0.7 
1.5 | —13.4 —13.4 —13.1 —12.0 — 3.5 —06 0.2 
2.0 | —25.8 —25.8 —23.5 —21.0 — 7.2 —2.1 —0.1 
3.0 | —42.4 —40.0 —32.7 —28.4 —11.1 —3.8 —0.5 
4.0} —52.3 —45.5 —35.3 —30.4 —12.8 —4.7 —0.8 
8.0 | —60.5 —45.0 —33.1 —28.2 —13.4 —6.2 —1.3 
20.0 | —46.3 —32.1 —23.3 —20.2 —10.6 —5.9 —1.8 


parallel, there are still two particles in each 
level—a proton and a neutron. 

The familiar formula for the kinetic energy of 
a large number N of particles obeying the 
exclusion principle and restricted to a volume Q, 
obtained by assigning two particles to each cell 
of size h® in phase space, is 


3/3\ 2/3 wh? 
K.E.--(-) (15) 
2MQ?/3 
where 
(16) 


This is seen to depend on the mass M of each 
particle, and on the volume 2 of the box, but 
not on its shape. However, a more detailed 
calculation for a small number of particles shows 
that the result does depend on the shape of the 
volume 2, and also on the type of boundary 
conditions used. Results are given below for a 
few simple cases. 

The simplest case is that of a rectangular box 
of sides a, koa, k3a and volume 2=k2k3;a*. Two 
types of boundary conditions are possible: 
periodic or zero boundary conditions. The wave 
functions are complex exponentials in the former 
case, and sines in the latter. The corresponding 
expressions for the kinetic energy are: 


(K.E.) exp = 4(kek3) iKy 


X +037/ks?]; (17) 
(K.E.) sin = (Rok3) 
Di? J. (18) 


Ko is defined by (16). The summation in (17) 
extends over zero, positive, and negative integral 
values of 1, m2, 3, while in (18) it extends 
only over the positive integral values. 
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TABLE IV. Values of No! = —107/5*Wa(k, 8). 


k 6=0.50 0.70 0.90 1.0 1.5 2.0 3.0 

1.0 

1.5 | 32.0 163 101 89 13.7 44.7 

2.0 | 16.6 8.4 5.6 5.1 6.6 13.0 91.7 

3.0 | 10.1 5.5 40 3.8 4.3 7.1 224 

4.0 8.2 4.8 3.7 35 3.7 5.7 15.2 

8.0) 7.1 4.8 40 3.8 3.5 4.3 9.1 
20.0 | 9.3 6.8 $7 $3 4.5 4.5 6.5 


TABLE V. Average energy per particle in Mev (from 
Bethe’s interaction) for a nucleus containing 45 particles. 


k 6=0.9 1.0 1.5 2.0 
3.0 202 76 101 169 
4.0 82 —16 25 126 
8.0 183 82 —2 61 
20.0 621 445 125 72 


TABLE VI. Average energy per particle in Mev (from 
Bethe’s interaction) for a nucleus containing 216 particles. 


k 6=0.7 0.9 1.0 
3.0 — 778 — 2298 — 2283 
4.0 — 1965 — 2871 — 2723 
8.0 — 1864 — 2388 — 1349 


Another simple case is that of a spherical box 
with zero boundary conditions for which 


K.E. = (4/32?) (19) 


where Ko is given by (16), and Z,,; (¢=1, 2, 3, 
--) are zeros of Bessel’s functions of order 
n=1/2, 3/2, 5/2, ---, the sum being taken over 
a sufficient number of the lower lying zeros. 
Figure 2 gives the ratio K.E./Ko as a function 
of the total number of particles for five different 
cases : 


I. Calculated from (15). 

II. Calculated from (17) with ke=k;=1. 
III. Calculated from (18) with k2=k3;=1. 
IV. Calculated from (18) with ko=1, k3=3. 

V. Calculated from (19). 


It is of interest to note that although for a given 
volume of the box the K.E. does not depend 
strongly on its shape (cubical, oblong, and 
spherical, respectively, in III, IV, V), it does 
depend markedly on the type of boundary 
conditions used (periodic in II, zero in III). 
For a very large number of particles the ratio of 
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the K.E. obtained from any of the formulae 
above to that given by (15) must necessarily 
approach unity. But for small N it is seen that 
zero boundary conditions, which lead to a greater 
concentration of particles towards the center of 
the box, give a considerably higher energy than 
the periodic boundary conditions, which for all 
values of N give results very close to those 
obtained from (15). 

Since in our simple model of the nucleus the 
density of particles is assumed to be constant 
over the nucleus, it is consistent to use the 


_ periodic rather than the zero boundary conditions 


on the wave functions in calculating the kinetic 
energy, the former giving a more uniform particle 
density. Since the dependence of the kinetic 
energy on the shape of the nucleus is slight, for 
simplicity expression (15) is used for all values 
of k and 6. Rewritten in terms of 6, (15) becomes: 


K.E. = (20) 
TABLE VII. Values of the function 


Wars(k, 6) = 3[—13.89W,(k, 8) —10.75We(k, 


in Mev. 


6=0.5 0.7 0.9 1.0 1.5 2.0 3.0 


—6.94 —6.94 —6.94 -694 —533 -3.26 —1.30 
—8.76 —8.76 —8.72 -851 -5.81 -—3.42 -—1.34 
—10.07 —10.07 —-965 -9.15 -585 -342 -—1.34 
—11.81 -11.39 -1001 -9.19 -5.57 -3.29 -—1.31 
—12.83  —11.58 —9.71 -8.76 -5.20 -3.12 -—1.27 
—12.84  —10.02 —7.81 -693 -3.94 -2.55 —1.12 
—9.33 —6.60 —-5.01 -—442 -2.64 -1.74 —0.86 


| 


TABLE VIII. Values of No! = —23.3/5?Wers(k, 5). 


k | 6=0.5 0.7 0.9 1.0 1.5 2.0 3.0 
1 13.41 684 4.14 3.35 1.92 1.79 1.98 
1.5| 1063 5.42 3.30 2.74 1.77 1.70 1.93 
2.0} 9.24 4.72 2.98 2.54 1.75 1.70 1.93 
7.88 4.17 287 2.53 1.84 1.77 1.98 
4.0) 7.26 4.10 2.96 2.66 1.97 1.87 2.04 
8.0; 7.25 4.74 3.68 3.36 2.60 2.28 2.31 
20.0; 9.98 7.20 5.74 5.27 3.88 3.34 2.99 


TABLE IX. Average energy per particle in Mev (from 
Rarita and Schwinger’s interaction) for a nucleus con- 
taining 9 particles. 


k 6=1.0 1.5 2.0 
1 38.2 —3.6 —4.1 
1.5 24.1 —7.9 —5.6 
2.0 18.4 —8.3 —5.6 
3.0 18.0 —5.7 —4.4 


where 


A? 8.4 
) Mev, (21) 
2n Mr,? ~ Cary)? 


with 1/x=2.18X10-* cm. 


6. CALCULATION OF TOTAL ENERGY FROM 
BETHE’S CONSTANTS 


Our approximation to the total energy of a 
heavy nucleus is obtained by adding (13) and 
(20): 

E=W(k, (22) 


In order to apply this to Bethe’s® neutral theory 
interaction with straight cut-off the functions 
F(r) and G(r) of Fig. 1 are replaced by the step 
function S(r) whose range ro is so chosen that 
the area under G(r) is the same as that under 
S(r). F(r) could be replaced by a step function 
of a different range, but A is numerically so 
much smaller than B in (7) as to make the 
contribution of W,(k, 6) to the total energy 
relatively unimportant, so that for convenience 
F(r) and G(r) are replaced by the same step 
function S(r) of range determined by 


The value of Kx is then obtained from (21): 
Kz=107 Mev, (24) 


and W2(k, 6) calculated from 


W2(k, 6) =3(7.7Wi(kR, 5) 
—102W2(k, 6)] Mev (25) 


is given for selected values of k and 6 in Table 
III. The value of N for which E=0 is given by 


Kp 107 
=— (26) 
4) 6°W2(k, 4) 


For values of N greater than No, E rapidly 
becomes more and more negative, very soon 
numerically exceeding the observed value of 


TABLE X. Average energy per particle in Mev (from 
Rarita and Schwinger's interaction) for a nucleus con- 
taining 216 particles. 


k 6=0.9 1.0 1.5 
2.0 — 1048 — 1137 — 891 
3.0 —1126 — 1146 — 830 
4.0 — 1062 — 1053 — 750 
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roughly —8N Mev. Table IV gives the values 
of N,' from (26) for selected values of k and 6. 
In this table the lowest value of No= (3.5)'=43 
is seen to occur for k=4, 6=1. Table V gives 
the values of the average energy per particle 
E/N in Mev when N=45, and Table VI gives 
the average energy per particle when N=216. 
If the range ro of the step function S(r) replacing 
F(r) and G(r) in the above calculations is 
increased so as to satisfy xro=0.7 instead of 
kro =0.6, the lowest value of No is further 
reduced to No=17. 

The dimensions of the ellipsoid specified by 
k=4, 6=1, xro=0.6 are: major axis=3.3X10-" 
cm, minor axis=0.83X10-" cm. The ellipsoid 
is of the order of magnitude of nuclear dimen- 
sions, but is excessively elongated. 


7. CALCULATION OF TOTAL ENERGY FROM RARITA 
AND SCHWINGER’S CONSTANTS 
Rarita and Schwinger‘ use the step function 
S(r) of range = 2.80 X10-" cm, i.e., 


kro = 1.285. (27) 


The corresponding value of Krs from (21) is 
Krs=23.3 Mev, (28) 


and Wrs(k, 5) calculated from 


Wrs(k, 6) 5) 
—10.75W2(k, 5) | Mev (29) 


is given for selected values of k and 6 in Table 
VII. The value of N for which E=0 is given by 


Krs 23.3 
6°Wrs(k, 4) 5) 


Table VIII gives the values of N¢! from (30) for 
selected values of k and 6. In this table the 
lowest value of No=(1.7)'=5 is seen to occur 
for k=2, 5=2. Table IX gives the values of 
the average energy per particle E/N in Mev 
when N=9, and Table X gives the average 
energy per particle when N= 216. 

The dimensions of the ellipsoid specified by 
k=2, 6=1.5, xro=1.285 are: major axis=6.7 
cm, minor axis= 3.4 cm. 


No'= (30) 


8. DISCUSSION 


Although the model used in the above calcu- 
lations has little meaning when applied to nuclei 
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containing as few particles as the values of No 
obtained above, nevertheless the inference may 
be drawn from these calculations that inter- 
actions of the type suggested by the neutral 
meson field theory of nuclear forces, although 
useful for describing the lighter nuclei, will allow 
stable nuclei having large spins, excessive binding 
energies, and deviating appreciably from spheri- 
cal shape to exist well within the range of the 
atomic weights of the periodic table, contrary to 
experimental evidence. 

The small values of No given by (26) are due 
to the large negative values of Wx(k, 5) obtained 
by using Bethe’s interaction. The feature of the 
interaction responsible for this is the excessive 
depth of the attractive tensor term as compared 
with the magnitude of the repulsive spherically 
symmetric term. This disparity in the relative 
strengths of the two parts of the interaction 
between a proton and a neutron is, however, 
required to give the observed binding energy of 
the deuteron, if the wave function for the ground 
state of the deuteron is not to deviate appreciably 
from spherical symmetry. 

The values of Wers(k, 6) obtained by using 
Rarita and Schwinger’s interaction are numeri- 
cally much smaller than those of Wa(k, 6) 
because in their interaction the tensor and the 
spherically symmetric terms are both attractive, 
and consequently the observed binding energy 
of the deuteron can be obtained with relatively 
shallow potential wells. In spite of this the 
values of No given by (30) are even smaller 
than those given by (26) because Kes in (30) is 
much smaller than Kz in (26). This is because 
the width of the potential well used by Rarita 
and Schwinger is twice the width of the square 
well introduced above to replace Bethe’s po- 
tential. The greater is the width of the square 
well used, the more important is the potential 
energy relative to the kinetic energy, since the 
particles can occupy a larger volume, and hence 
have less kinetic energy, for the same value of 
the potential energy. 
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APPENDIX 
The analytical expressions for the functions W,(k, 4) 
and W,(k, 6) defined by (14) are given below. 
k=1 (Prolate Ellipsoids) 
W,(k, 6)=1 
B(k)5-*+ C(k)5-* 
— for 
B(k)d-*+- C(k)5-* for 


for 


where 
C(k) = 4[3+2k-*+ 3k?(k?— (k) J, 
D(x) = 34(x-*— 
+194 1)¥?— — cos™ x, 


and 

W2(k, 6) = E(k), for 
= F(k)5-*—G(k)5 *+H(k, 5) for 
= for 

where 


E(k) = 1) cos“! (k), 
F(k) = 4) (k2— cos“! 
G(k) = 14 [3k2(k® — 2) 1)-8 cos“! (k-2) 
+ 
H(k, 8) =3¢k2(k?— 3(k2—4) 
X +154 — 1)-1/2( — 1/2 
— (k2 +44) — 1) — 1) 3/2 
+3¢ 2) 1) — 1) 5/2, 


k=1 (Oblate Ellipsoids)® 


Wi(k, 6) =1 for 
for 


where 


C’(k) = 4[3+2k*+3k*(1 cosh (k-)], 
D' (a) =34(1 (1 — 3/2 
+154(1—x-*)"2— (x4 — cosh™ x 


and 

W2(k, 5) = —E’(k) for 6k", 
= — E’(k)+H'(k, 6) for 
= — for 


where 


E’(k) 1) — 3k2(1 — cosh™ (k™), 
F'(k) = %k*3[ cosh 
G'(k) = [3k?(k? cosh™ 
+ (3k? —4+44k-*) (1 —k*)), 
H'(k, 5) = 3gh?(1 (hk? — 2) 
—3(k?—4)k-*/3§-4— 8] cosh! 
— 154 (1 — 
+1¢ (k®+44) (1 — 3/2 
+36 (k®—2)(1 — 8/2, 


The three separate ranges of validity of all the above 
results correspond, respectively, to the range 7o of the step 
function S(r) being greater than the maximum dimension 
of the ellipsoid, intermediate between the maximum and 
minimum dimensions, and less than the minimum dimen- 
sion. 

As 6—« the values of W,(k, 5), Wo(k, 6), and 
W2(k, 6)/W,(k, 6) all approach zero for all values of k. 

As k->~ for any finite value of 6 the values of W,(k, 8) 
and W2,(k, 6) both approach zero, while their ratio 
W2(k, 6)/W,(k, 6) approaches 2. 

As k-0 for any finite value of 6 the values of W,(k, 8) 
and W2(k, 6) both approach zero, while their ratio 
W2(k, 5)/Wi(k, 6) approaches —1. 

W,(k, 5) is always positive, while W2(k, 5) is positive for 
prolate ellipsoids (k>1), negative for oblate ellipsoids 
(k<1), and zero for a sphere (k=1). 

The values of W,(k, 6) and W2(k, 6) calculated from the 
above formulae for selected values of k anc 6 are given in 
Tables I and II. 


* The formulae for the oblate case were also obtained independently by W. Rarita and J. Schwinger. I am indebted 
to them for yp ane ag he me their results and the following numerical values of W:(k, 5) and W2(k, 6) some of which 
i 


supplement those already listed in Tables I and II: 


é k 1 0.9 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 

1 W,i(k, 5) 1.000 1.000 0.999 0.996 0.988 0.973 0.945 0.896 0.805 0.624 
W2(k, 6) 0.000 —0.082 —0.172 —0.266 —0.360 —0.454 —0.542 —0.612 —0.642 —0.566 

0 W1(k, 4) 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 
W2(k, 6) 0.000 —0.082 —0.172 —0.268 —0.370 —0.480 —0.594 —0.714 —0.832 —0.940 
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A Note on Exchange Tensor Forces in Heavy Nuclei 


G. M. Vo_korr* 
Department of Physics, University of British Columbia, Vancouver, British Columbia, Canada 
(Received June 9, 1942) 


Short range spherically symmetrical exchange forces between elementary nuclear particles 
were first introduced into nuclear theory to account for the observed saturation of nuclear 
forces. More recently additional tensor forces between elementary particles have been proposed. 
The object of this note is to point out that if tensor forces, whether of the ordinary or the 
exchange type, are assumed to be the dominant forces between nuclear particles, they will not 
lead to saturation, but will permit the existence of nuclei with excessive values of binding 
energy, spin, and isotopic number. The presence of a sufficiently strong spherically symmetrical 
exchange term in the nuclear interaction thus seems to be required to counteract this effect of 


the tensor forces. 


HE Hartree method is used to obtain an 

estimate of the nuclear energy.! The nu- 
cleus is described by the completely antisym- 
metrical wave function 


Wilf) 


vilév) ++: 


The wave functions y,(£;) describing the indi- 
vidual particles are taken to be products of three 
factors, one depending only on the space co- 
ordinates r; of the particle, the second only on its 
spin, and the third only on its isotopic spin, i.e., 
specifying whether a given particle is a proton or 
a neutron. Let the spin wave function be a or 8 
depending on whether the spin component in a 
given direction Z is +} or —}3, and the isotopic 
spin wave function be p or ” according to whether 
a given particle is a proton or a neutron. Then a 
given y;(é;) will consist of some spatial wave func- 
tion ¢,(r;) multiplied by one of the four possible 
combinations: a(s;)p(7;), B(s;)p(7;), a(s;)n(r;), 
B(s;)n(7;). Let the interaction between a pair of 
particles be denoted by V(é;, &;), and the kinetic 
energy of each particle by 7;. The total energy of 


* This note was written while the author was visiting at 


' the Department of Physics, University of California, 


Berkeley, California. It is a pleasure for him to record his 
gratitude for the hospitality extended to him. 

1Cf. H. A. Bethe and R. F. Bacher, Rev. Mod. Phys. 8, 
82 (1936). 


the nucleus on this model is then given by:! 


N NN 
Ti+} 


i=1 i=Lk=1 


f f V1, (2) 


The integrations in (2) imply summations over 
the spin and isotopic spin as well as integrations 
over the space coordinates. After the explicit ex- 
pression for V(1, 2) has been substituted in (2) 
and the summations over the spin and isotopic 
spin coordinates have been carried out two types 
of space integrals may appear in the result: the 
“direct” integrals 


Jaw f f $:*(1)bx*(2) (3) 
and the “exchange” integrals 
Kae f f (4) 


Using plane waves for ¢,(r;) the following prop- 
erties of these integrals may be demonstrated.! 


The sum 


N N 
f f p(1) (5) 
im1 k=l 
is, for small R, proportional to N?, while the sum 
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is proportional to N. Here 


p(t) |? 


i=1 


and 


N 


are the total particle density, and Dirac’s ‘‘mixed 
density,’’ respectively. The former may be taken 
constant over the nucleus and of the order NR 
where R is the nuclear radius. The latter is of 
order NR~ when r; and r- coincide, and falls off 
rapidly as soon as |rfi:—r2| becomes greater 
than RN-}. 

In order to bring out the characteristic fea- 
tures of spherically symmetrical and tensor forces 
responsible for the difference in their behavior 
with respect to saturation we substitute for 
V(é:, 2) in (2) the following four expressions 
in turn: 


Ordinary Exchange 
Spherically symmetrical J (r) — (r) (8) 
Tensor J(r)Si2 — (1) Sie. 


Here J(r) is a spherically symmetrical potential. 
Siz is defined by 


Si2=3 


(9) 


which reduces in the case of the two spins point- 
ing in the same direction to 


Si2=3 cos? @—1. (10) 


The exchange operator — }#,-*2 has the value +1 
for a state antisymmetrical in the charge (e.g., 
the ground state of the deuteron), and —}4 for a 
state symmetrical in the charge (e.g., two pro- 
tons, or two neutrons). Also by analogy with 
similar results' for the case of the ordinary spin 
operator 


t1°%2p(1)p(2) = p(1)p(2), 
t2p(1)m(2) = — p(1)m(2)+2n(1)p(2). 
With V(é1, &)=J(r) the direct integrals Ji 
all have the same sign as J(r) irrespective of the 


values of the spin and the isotopic spin of the 
individual particles, so that (5) is directly ap- 


(11) 


plicable resulting in lack of saturation for an 
attractive J(r). The introduction of the addi- 
tional exchange operator — }¢,-t2 helps to avoid 
this difficulty by making the direct integrals 
between like particles contribute positive amounts 
to the total energy counteracting the negative 
contributions of the direct integrals between un- 
like particles. In virtue of (11) the contribution 
of all the direct integrals from (2) will now be 
given not by (5) but by 


—pn(2) Jdridre, (12) 


where p, and p, are now separately the proton 
and the neutron densities. If the individual wave 
functions in (1) are chosen in such a way that 
there is always both a proton and a neutron in 
each space-spin state, then p,—pn, and hence 
(12), will be identically zero. For any other dis- 
tribution of protons and neutrons among the 
space-spin states (e.g., the extreme case of a 
nucleus consisting of only protons or only neu- 
trons) such that p,—p, is always of one sign, the 
total contribution of the direct integrals given by 
(12) is always repulsive for an attractive J(r). 
When (12) vanishes the binding energy is due 
entirely to the contributions of the exchange in- 
tegrals Kx, so that according to (6) saturation 
will result. Saturation in this case therefore de- 
pends on the fact that the direct integral of the 
interaction —4¢,-t2J(r) between two like par- 
ticles is necessarily repulsive for an attractive J(r). 
The feature of the tensor forces which produces 
lack of saturation is the fact that by a proper 
choice of individual particle wave functions in 
(1) it is always possible to make all the direct 
integrals J give an attractive contribution. 
Take all the nuclear particles with their spins 
parallel, so that Siz. can be replaced by (10). The 
ordinary tensor interaction is then given by 


V(r12) = J(r)(3 cos? @—1). (13) 


The contribution of all the direct integrals will 
be given by (5), and can always be made attrac- 
tive by choosing a particle distribution concen- 
trated along the polar axis (prolate) if J(r) is 


negative, or in the equatorial plane (oblate) if 
J(r) is positive. To show that the exchange 
tensor force —}*1-t2V(fi2) also leads to lack of 
saturation choose the individual particles not 
only with all their spins parallel, but also all of 
the same kind, so that —4}¢,-t2=—43 for all 
pairs, which is equivalent to reversing the sign 
of J(r) in (13). With spherically symmetrical 
forces such a reversal of sign necessarily leads to 
a repulsive total contribution proportional to N?. 
With tensor forces the reversal of sign of J(r) can 
be compensated for by a suitable change (prolate 
to oblate) in the distribution of particles. There- 
fore if tensor forces alone are assumed to act 
between elementary nuclear particles they would 
allow the existence of heavy nuclei composed of 
only protons or only neutrons, deviating appre- 
ciably from spherical shape and having excessive 
values of the spin and the binding energy. 

The presence of a sufficiently strong spherically 
symmetrical exchange term of the right sign 
will counteract this tendency of the tensor forces 
to produce non-saturation. If the spherically- 
symmetrical exchange term is assumed to be the 
dominant one, while the tensor force is considered 
as a perturbation, then it will favor nuclei com- 
posed of proton-neutron pairs in identical space- 
spin states, as was seen above. The exclusion 
principle will further favor the presence of two 
protons of opposite spin, and two neutrons of 
opposite spin in each space state. For such a 
combination of individual particle wave functions 
the direct integrals J, of the tensor part of the 
interaction will all cancel out pair by pair, as in 
the spherically symmetrical case, because of the 
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relations 


T1, 72 81, 82 


X t2a(1)p(1)a(2)p(2) =3 cos? @—1, 
14 
a(1)p(1)8(2)p(2) 


TL, 81,82 


X t2a(1)p(1)8(2)p(2) = —3 cos? 6+1, 


and similar relations for neutrons, and likewise 
for the ordinary tensor force with the +;-+2 oper- 
ator left out. Since the contribution of the direct 
integrals in this case vanishes saturation again 
becomes possible. 

In a recent paper? the influence on heavy nuclei 
of non-exchange interactions combining spheri- 
cally symmetrical and tensor terms proposed by 
Bethe and by Rarita and Schwinger was investi- 
gated. It was concluded that such interactions 
would lead to difficulties even for relatively light 
nuclei. On the basis of the above discussion it 
would seem that the inclusion of the exchange 
operator —}¢,-t2 in Rarita and Schwinger’s po- 
tential will overcome non-saturation and related 
difficulties, since the spherically symmetric term 
in their interaction has a depth of 13.39 Mev, 
while the depth of the tensor term is only 10.75 
Mev. On the other hand in Bethe’s interaction 
the tensor term having a depth of 102 Mev far 
exceeds the 7.7 Mev spherically symmetrical 
term, so that the inclusion of an exchange opera- 
tor will not prevent non-saturation. 

It is a pleasure to thank Professor J. R. Oppen- 
heimer for his interest in this work, and for 
valuable discussions. 


2G. M. Volkoff, Phys. Rev. 62, 126 (1942). 
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The Absence of the Mg X-Ray Satellite Intensity Anomaly* 


F. R. Hirsu, Jr. 
Pasadena, California 


(Received January 19, 1942) 


The absence of the M@ x-ray satellite intensity maximum is shown, and explained in terms 
of the relative probability of Auger transitions, the relative probabilities being determined by 


means of diagram line widths. 


INTRODUCTION 


N a former study of the Ma; and M§ satellites,' 
the writer especially noted the similarity of 
the ‘“semi-Moseley”’ graphs for the maxima of 
these two groups: they are practically identical. 
Moreover, in a subsequent paper? the writer 
pointed out how similarly shaped the Ma; and 
MB lines of Bi(83) were (Fig. 7, reference 2). 
In that same paper,’ the present writer proposed 
that the Ma; satellites were due to the transition 
MyMw.v—NyuMiv,v effected by the radiation- 
less transition My,— My with ejection of an Niy,yv 
electron; that the M@ satellites were due to the 
transition Niv.v effected by the 
radiationless transition Myy—>Mjy with ejection 
of an Njy,y electron. Thus the satellites of both 
the Ma; and M8 groups could be accounted for 
by the Auger effect. Moreover, the very similar 
nature of the two, single-electron transitions 
between doubly-ionized states, just mentioned, 
practically ensures the similarity of the semi- 
Moseley graph groups for Ma; and M8 as well 
as that of the over-all line shapes. In Munier, 
Bearden, and Shaw’s article,* the Ma and MB 
groups for W(74), (p. 541), are very similar, 
except for fine detail. 

With this situation in mind, the present 
writer investigated the plates from which the 
measurements of reference 1 were made. It was 
simple to show photographically that the Ma; 
lines revealed a maximum of intensity‘ for their 
satellites near the proper atomic number. In the 
case of the M8 satellites, the energy curve for 


*The absence of the Mf satellite intensity anomaly 
was briefly commented on in a Letter to the Editor of 
The Physical Review [Phys. Rev. 59, 766 (1941) ]. 

1F.R. Hirsh, Jr., Phys. Rev. 38, 914 (1931). 

?F. R. Hirsh, Jr., Phys. Rev. 50, 191 (1936). 

3J. H. Munier, J. A. Bearden, and C. H. Shaw, Phys. 
Rev. 58, 537 (1940). 

4F. R. Hirsh, Jr., Phys. Rev. 57, 662 (1940). 


the radiationless transition intersects 
and passes above the ionization-energy curve for 
the Niy,y shell (necessarily for Z+1 as an inner 
M electron is missing), at Z=84 (see Fig. 1). 
Thus we might expect the Auger effect to occur 
and doubly ionize atoms below Z = 84, producing 
an Auger intensity maximum (anomaly), as in 
the case of the Ma, satellites. 


EXPERIMENTAL 


The plates for the M@ lines! (all lines except 
those for Pt and Ir were taken with a calcite 
crystal ; the two named were taken with a quartz 
crystal) were photometered by using a Koch 
microphotometer; the constants were secured 
and density plotted as a function of wave-length 
for each x-ray line. The range of the present 
wave-lengths, from 3.7 to 6.0A, has been studied 
by the writer® who has shown density to be 
proportional to exposure (intensity constant), 
over a considerable range of densities. The 
present writer is not alone in this result. Having 


Atomic Number 


Fic. 1. Coster-Kronig diagram for the satellites of the 
MB line: energy released by the radiationless transitions 
Mjy and plotted as a function of atomic 
number (full line); ionization energies for the various NV 
—_— plotted for Z+1 as an M electron is missing (dashed 
ines). 


5 F. R. Hirsh, Jr., J. Opt. Soc. Am. 25, 229 (1935); 28, 
463 (1938). 
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Th,(90) 


Bi,(83) 


«| Pb(82) . TL (8) 


Fic. 2. Density (intensity) plots for MB vig | lines of the elements from U(92) to Ir(77). Vertical lines indicate visual 


measurements O 


secured the true over-all intensity shapes for each 
line (see Fig. 2), a suitable and ‘“‘reasonable”’ 
background for each satellite group was drawn in 
with a spline and the satellites were delineated. 
Next, the microphotometer magnification was 
secured, the A’s (satellite separations) were com- 
puted easily from the Ad values of reference 1, 
and the spectrograph radius was known to be 
36.728 cm. These A values are seen to be quite 
accurately related to the satellites of Fig. 2, with 
one exception: the satellite MB‘' is seen from the 
density (intensity) plots to be double (two com- 
ponents) for Pb(82), T1(81), and Au(79), showing 
without a doubt that the semt-Moseley graph group 
for the MB satellites is the exact counterpart of the 
semt- Moseley graph group for the Ma, satellites. 
In each case the M£*' satellite is seen to be 
the most intense component, and easily the best 
defined. For this reason the intensity of the Mp" 
group, with respect to the parent line, was taken 
as a measure of x-ray line energy produced by 


maxima on the original plates (reference 1). 


double ionization as compared with the x-ray 
line energy produced by single ionization. These 
data are presented in Table I. While there should 
be an increase in the relative energy (ratio S/P) 
starting just below Z = 84, coming to a maximum 
and falling off with decreasing atomic number, 
it is quite evident that the intensity maximum is 
not present for the MB satellites. 

The absence of the intensity maximum for the 
MB satellites is made very obvious in an illustra- 
tion (made from these same plates taken by the 


TABLE I. Relative energy (by areas) for the Ms" satellite. 


Parent line area 
(P) 


MB* satellite area 
Element (S) Ratio S/P 
U(92) 0.01 X10 sq. in. 0.50 < 10 sq. in. 0.020 
Th(90) 0.01 : 0.44 0.023 
Bi(83) 0.01 0.49 0.020 
Pb(82) 0.03 3.22 0.025 
T1(81) 0.045 1.97 0.023 
Au(79) 0.07 2.76 0.025 
Pt(78) 0.02 0.78 0.026 
Ir(77) 0.028 0.86 0.032 
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present author') on page 305 of H. E. White’s 
book, Introduction to Atomic Spectra. Here the 
satellites of Ma; show an unmistakable maximum 
of intensity at Pb(82), while it is absent for the 
MB satellites. 


DISCUSSION OF EXPERIMENTAL RESULTS 


One fact that is evident in Fig. 1, is that the 
radiationless transition My— My can eject an N; 
electron at Z=84, and below that atomic num- 
ber. Thus there are two Auger transitions possible 
at and below Z=84: one which will produce the 
satellite group, (Min- and one which will 
not, (My—Myy). If it could be shown that the 
radiationless transition was more prob- 
able than the radiationless transition My,;— 
the first transition would rob the second to its 
effective exclusion. This seems to be actually the 
fact, and will now be discussed. 

A very interesting fact, which is not necessarily 
essential to this argument, but which tells us 
what intra-level transfers are effective as radia- 
tionless transitions, is that all the known radia- 
tionless transitions, which cause the known x-ray 


satellite intensity anomalies, are permitted by the 


x-ray quantum number selection rules: Aj= +1, 0; 
Ak= +1; An is arbitrary. For (La; and 
LB: satellites), Aj= —1, and Ak= —1; for 
(LA; satellites), Aj =0 and Ak = —1;for Min—My 
(Ma; satellites), Aj= —1 and Ak = —1. Similarly 
for the two radiationless transitions under dis- 
cussion: for My—Myy, 4j=0 and Ak= —1; for 
My, Aj= —1 and Ak= 

Theoreticians might immediately raise a pro- 
test as to the significance of these last observa- 
tions. However, it will be remembered that in 
1922 D. Coster® was searching for radiations 
corresponding to the transitién Li—Lm (per- 
mitted by selection rules). He failed in this first 
search ; he later reported’ that radiation had been 
observed very faintly for a few atomic numbers. 
We now know why this search failed at first: 
I;—Ln is much more probable as a radiationless 
transition (see reference 7). Such electron jumps 
are intra-level transfers as compared with the 
usual inter-level transfers which give rise to the 
ordinary x-ray diagram lines. Thus it would seem 
evident that these intra-level transfers, while prob- 


®D. Coster, Phil. Mag. 43, 1089 (1922). 
7D. Coster and R. DeL. Kronig, Physica 2, 13 (1935). 
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_ Fic. 3. Energy level diagram showing relation of radia- 
tionless transitions to radiation transitions in interpreting 
Auger probabilities in terms of line widths. 


able to some extent as ordinary permitted transitions 
yielding radiation, are much more probable as 
radiationless transitions. The significance of the 
radiationless transitions being permitted by selec- 
tion rules is that the selection rules point out where 
the radiationless transitions are sure to occur. 
With this last fact well in mind, let us proceed 
to the solution of the problem caused by the 
absence of the M8 x-ray satellite intensity 
anomaly. Figure 3 shows the situation as regards 
the radiationless transitions: all are permitted 
by selection rules; two result in ionization of the 
Mwy shell (Min—My ionizes the My shell), but 
the Min— Mwy transition alone could cause en- 
hancement of the M@ satellite group, as has been 
previously discussed in this paper. The facts 
are as follows: My;—My and My—-My both 
broaden the Mi level by decreasing its mean 
ionized life; Mu—My broadens the My level 
similarly by decreasing its mean ionized life. If 
the My level is broader than the My level, we 
would immediately see that My—Myry has an 
extremely high probability of occurrence. Ram- 
berg and Richtmyer® ascribe practically the entire 
breadth of such levels to Auger transitions. How 
may we measure the breadth of the Mn and 
My levels? Fortunately, there are two radiation 
transitions starting from these levels and ending 
in a common level: LB, and L; (see Fig. 3). 
Any difference in the width of the Lf, line as 


8 E. Ramberg and F. K. Richtmyer, Phys. Rev. 51, 913 
(1937). 
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TABLE II. Corrected line breadths (ev). 

Yb(70) Lu(71) Ta(73) w(74) Re(75) Os(76) Ir(77) Pt(78) Au(79) TH(81) 
LBs 13.9 _— 14.1 15.1 16.4 17.1 17.7 18.4 19.6 20.4 
LB; 12.9 12.8 11.7 12.3 13.6 14.7 14.9 16.1 17.2 18.4 
Diff. 1.0 —_— 2.4 2.8 2.8 2.4 2.8 2.3 2.4 2.0 
compared with the L@; line, must be due to the From (5), for the Jj level: 
greater width of the My level compared with AE h . 
the Min level. The breadth of a diagram line iaaliiiat teases (6) 
has been shown to be equal to the sum of the For the Mi level: 
breadths of the two energy levels between which 

~h{ yonn-Mw) +yOnu—-My)}. (7) 


the electron transition takes place.® 
One way of expressing the Bohr-Heisenberg 
indeterminacy principle is: 


AE,- AT, =h, (1) 


where AE; is the uncertainty in measuring the 
energy of a certain atomic phenomenon, and 
AT; is the uncertainty in measuring the time of 
its occurrence. Equation (1) is dimensionally 
equivalent to the more commonly quoted form. 
Let AE; be considered to be the energy breadth 
of a given electron energy level, and AT; the 
mean ionized life: we then see that the inde- 
terminacy principle applies to our problem. Let 
I',=1/AT), where IL; is the reciprocal of the mean 
ionized life. Then 


AE, ~h-T). 


(A) (R) 


(2) 


But 
(3) 
which says that for a particular level I’; is the 
sum of all the y’s contributed by the Auger 
transitions, plus the sum of all the y’s con- 
tributed by ordinary radiation transitions.* 

Now Ramberg and Richtmyer® have shown 
that in the M and N levels where the Auger 
transitions are the chief contributors to the 
breadth of levels: 


(4) 
(A) 
i.e., 
(A) (R) 
Then from (2) and (4) we have 
AE, ~h > (5) 


(A) 


°V. Weiskopf and E. Wigner, Zeits. f. Physik 63, 54 
(1930) ; 65, 18 (1930). 


For the L@, line Weiskopf and Wigner have 
shown :° 


AEMy+AExy. (8) 
For the L8; line: 
= AEM + (9) 
Subtracting (9) from (8) we have: 
AEts,— = AEMy — AEMu. (10) 


Substituting in (10) from (6) and (7) we then 
have: 


AE1g,—AE1g;~h { 


(11) 


But AExg,—AEzs; is always + (positive) be- 
tween Z=70 and Z=81 as is shown in Table II, 
data from the thesis of Dr. John N. Cooper,” 
written at Cornell University in 1940. I am in- 
debted to Dr. Cooper for an admirable presenta- 
tion of the subject of line breadths. 

This, then, affirms the idea that Wy— My is 
more probable, according to line widths, than 
Mm—-My and Miy—My combined; this fact 
prevents the occurrence of the M@ satellite in- 
tensity anomaly. 
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AUGUST 1 AND 15, 1942 PHYSICAL REVIEW VOLUME 62 


The Absorption Spectra of CdF and SnCl 


CuarLes A. FOWLER, JR. 
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a (Received May 21, 1942) 


Cadmium fluoride, vaporized in an electric furnace, shows no absorption in the region AA1950— 
7000 at temperatures below 1350°C, but above this temperature a band system of CdF appears 
at \A2716-2924. Strong continuous absorption is observed at higher temperatures. The spec- 
(6) trum of SnCl is photographed in absorption for the first time. The spectrum consists of the two 

band systems observed in emission, a new system in the farther ultraviolet, and two continua. 


Its similarity to the absorption spectrum of the homologous molecule SnF is discussed. 


(7) 
iia I. CdF their intensities decreasing from band to band 
EARSE and Gaydon have shown! that the in the direction of shorter wave-lengths. Although 
only spectrum previously attributed to CdF the first two bands of longest wave-length are the. 
(8) was actually due to the presence of CaF as an strongest of the system, they have very diffuse 
impurity. They have listed six other bands heads. The third has a fairly definite head while 
between 42819 and 43005 as possibly due to CdF those of the fourth, fifth, and sixth are of increas- 
(9) inasmuch as they could not be traced to any i" Sharpness. Each of the third, fourth, and 
known impurity, but have hesitated to attribute fifth bands has a weak satellite associated with . 
them definitely to this molecule because of their the strong head. A narrow continuum with its 
10) weakness and general behavior. Since those in- ™@ximum at \2824 and a group of three closely 
vestigations were carried out in arc emission, it SPaced heads apparently not belonging to the | 
len was thought that an attempt to observe the principal system constitute the remaining fea- 
spectrum in absorption might give less ambigu- tures of the spectrum. Measured wave-lengths of | 
ous results. the heads and the corresponding wave numbers 1 
The absorbing vapor was produced by heating in vacuum are given in Table I together with an 
[1) a small amount (~ 10g) of CdF2 in a carbon-tube assignment of probable values for v’—v’’. The 
furnace previously described.2. The spectral individual quantum numbers cannot be assigned 
m4 region 1950-7000 was surveyed under low With certainty because only a single progression | 
= dispersion and regions showing band structure 8 Observed, the higher members of each sequence i 
ry : i . The values of 
were photographed in the first order of the being covered by the first member . 
call 30,000 lines-per-inch, 21-foot concave grating v’—v"’ in Table I are the most probable when one 
fa with average dispersion of 1.3A/mm. The ex- considers the great likelihood of a v’ progression ' 
: perimental technique has been fully described in absorption, despite the fact that the trend of 
is - 
in a previous article.’ 
sale At the lower furnace temperatures no bands - 
| were observed but at about 1350°C some narrow a. 
n- 
regions of absorption around \2800 appeared. ,_ ee 
With the vapor at 1600° the bands were quite | 
strong (Fig. 1), but at higher temperatures con- a | 
tinuous absorption, coming in from shorter 
wave-lengths, began to cover up the region. 
High dispersion spectrograms revealed six bands 
. shaded to the red and spaced almost equally, 
he IR. W. B. Pearse and A. G. Gaydon, Proc. Roy. Soc. Fic. 1. Low dispersion absorption spectra of (a) CdF, 
50, 711 (1938). (6) SnCl with the vapor density high, (c) SnCl with the 
2F. A. Jenkins and G. D. Rochester, Phys. Rev. 52, 


1135 (1937). 
3C. A. Fowler, Jr., Phys. Rev. 59, 645 (1941). 


vapor density low. High dispersion spectrograms of (d) 
part of the CdF band system and gh de ower energy 
component of the C*II, X*II system of SnCl. 
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TABLE I. Band heads of CdF. 
0 2923.5A 34196 
509 
1 2880.6 34705 
514 
2 2838.5 35219 
522 
3 2797.1 35741 
529 
4 2756.3 36270 
535 
5 2716.2 36805 
Satellite 
heads 
2 2837.2 35236 
519 
3 2796.0 35755 
526 
4 2755.5 36281 
Unclassified 
heads 
— 2788.2 35855 
22 
2786.5 35877 
20 
— 2784.9 35897 


the first differences indicates otherwise. Such a 
trend of the differences is possible, however, if 
the rotational and vibrational constants have 
certain, nearly equal values in the two electronic 
states involved in the transition. 

A visual inspection shows that the over-all 
extent of each band is not large and that the 
rotational structure is unresolved except near the 
tail of the band, indicating a very small value for 
B’—B”", so that the band heads are formed at 
high J values and at considerable distances from 
the band origins. Furthermore, the distance 
between each band head and its origin will 
decrease rapidly with increasing v’. Therefore, 
although the first differences of the origin wave 
numbers will decrease with v’, the first differences 
of the band heads may increase. Two features of 
the CdF bands support this assumption. The low 
(long wave-length) members of the progression 
whose heads should be formed at very high 
values of J show an almost headless nature, while 
the heads of the higher members are increasingly 
sharp. Furthermore, the third, fourth, and fifth 
bands possess satellite heads of increasing sharp- 
ness, which would not be expected in the lower 
members owing to the high J values at such 
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heads. The diffuseness and lack of resolution is 
increased by the approximate superposition of 
the many bands of each sequence, superimposed 
because of the near equality of w,’ and w,’’. It is 
evident from the foregoing that all of the v’—v’’ 
might well have values less by one than those 
given in Table I. 

The vibrational constants cannot be accurately 
determined because only one progression of heads 
is measurable. Moreover, the first differences are 
poor representations of the AG(v+}) when B’ 
and B” are nearly equal. The data indicate that 
we’ and w,”’ are almost equal and somewhat 
greater than 535 cm. A comparison of the 
normal state vibration frequencies of CdCl 
(330.5 CdBr (230.0 and CdI (178.5 
cm~') with those of the various alkaline earth 
halides‘ reveals strikingly similar trends, and 
predicts for CdF a vibration frequency of about 
540 cm. Such a comparison also predicts for 
ZnF a normal state vibration frequency of about 
630 cm~, in close agreement with the 620 cm~! 
reported by Rochester and Olsson.°® 

No other bands were observed between \1950 
and A7000. It is to be noted that although they 
lie in the same general region, there is no cor- 
relation in wave-lengths or first differences 
between these bands and those of Pearse and 


Gaydon. 


II. SnCl 


The spectrum of SnCl has been observed only 
in emission where it shows two band systems, a 
2A, system® at AA3486-3910 and a 
system’ in the region \A2830—3405. The *II level, 
with doublet separation of 2360 cm, is common 
to both systems and is assumed to be the normal 
state. The occurrence of a *II lower state in the 
absorption spectrum of the SnF molecule® 
implies the same type of normal state for all of 
the halides of Si, Ge, Sn, and Pb. 

This investigation was undertaken to test the 
possibility of observing any shorter wave-length 
band systems analogous to those of SnF by 
observing the spectrum in absorption. 


4G. Herzberg, Molecular Spectra and Molecular Structure, 
I (Prentice- Hall, 1939), pp. 483-494. 
’G. D. Rochester and E. Olsson, Zeits. f. Physik 114, 
495 
W. F. C. Ferguson, Phys. Rev. 32, 607 (1929). 
TW. Jevons, Proc. Roy. Soc. A110, 365 (1926). 
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With furnace temperatures ranging from below 
the boiling point of SnCl. (605°C) to 1000°C a 
weak, general absorption due to SnCl. was ob- 
served, but only by superheating the vapor to 
higher temperatures did the bands and continua 
characteristic of SnCl become evident. At 1200°C 
the higher energy component of the B?2<~—X?II 
system is of considerable intensity, but the 
lower energy component is hidden by a strong, 
sharply defined continuum covering the region 
\A3000-—3500. A second continuum extends from 
2500 to the ultraviolet limit of observation at 
1950. At these temperatures the salt was quickly 
vaporized, and the vapor density slowly de- 
creased by diffusion from the carbon tube. The 
rate of this diffusion can be governed by regu- 
lating the temperature, and a series of low 
dispersion spectrograms showed that the con- 
tinua decrease in extent as time goes on. At one 
stage the more refrangible component of the 
A°?AX?Il system (A3488) could be observed 
under the edge of the weakening continuum. In 
the same manner the short wave-length con- 
tinuum in retreating to the ultraviolet revealed a 
new band system at \A2250—2450 consisting of a 
number of heads shaded to the ultraviolet 
(Fig. 1). 

High dispersion plates of this system were 
taken with a furnace temperature of 1300°C and 
were exposed only during a predetermined 15- 
minute interval for which the vapor density was 
suitable for the bands to appear. The heads of 


.the system fall readily into two groups, three 


faint bands at AA2436, 2416, 2393 and six 
stronger ones between \2249 and 2326. These 
constitute the two components of a doublet 
system, the heads of the higher energy com- 
ponent being represented by the equation 


v=43656.1+ 392.8’ — 1.40u" — 352.8u'’+1.00u'”. 


TABLE IT. Wave-lengths and intensities of the band heads 
of the CX system of SnCl. 


2,0 2248.86 8 0,1 2307.42 5 
1,0 2268.62 10 2436.38 0 

2393.08 1 0,2 2326.16 2 
0,0 2288.88 9 1,3 2323.80 3 
2415.67 1 


The observed wave-lengths and the estimated 
intensities are given in Table IT. 

A difference of 2292.2 cm~' exists between the 
two 0,0 heads. Inasmuch as the *II normal state 
has a separation of 2360.1 cm~', the new upper 
state must also be double, i.e., *II or *A, since a 
*~ upper level would produce two sub-bands 
whose separation equals the lower state interval. 
In both and 2A, *II transitions only two 
sub-bands are expected because of the selection 
rules, and the wave number difference must be 
the difference (or the sum) of the doublet separa- 
tions of the two levels involved. Hence the new 
upper state is a doublet with a separation of 67.9 
cm~!. That the heads are apparently single 
implies that the upper state is II, but this is not 
conclusive because of the low intensity of the 
bands. 

The intensity ratio of the two sub-bands 
depends upon the temperature of the vapor and 
the splitting of the normal state if a Boltzmann 
distribution is assumed, and is given by e~#/*7 =$§ 
for the temperature of 1300°C. This agrees well 
with a visual estimate of the 0,0 bands. Com- 
parison of the intensities of the two sub-bands 
is very difficult, however, because of a rapid fall 
of background intensity toward shorter wave- 
lengths. This is an inherent characteristic of the 
He continuum used as the source, and is here 
accentuated by a partially overlying absorption 
continuum. 

The analogy between the electronic states of 
SnCl and SnF is quite obvious. Both possess *II 
normal states. The A*A and B*> states of SnCl 
correspond to the B?A and A?’ states, respec- 
tively, of SnF as pointed out by Jenkins and 
Rochester,? while the new state of SnCl is 
probably “II corresponding to the C?II state of 
the fluoride. Failure to observe the D state can 
be attributed to the strong overlying con- 
tinuum. The two continua, because of their close 
association with the band systems, must be 
attributed to SnCl. They indicate the existence 
of the repulsive states E and F, and these, like 
the stable excited states, have somewhat smaller 
energies than their SnF analogues. 

In conclusion the writer wishes to thank 
Professor F. A. Jenkins for his assistance and 
criticism in these investigations. 
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The theory of the action of a recently developed 8-ray spectrograph is discussed briefly and an 
account is given of several improvements. In particular a baffle system of greater resolving 
power and a G-M counting system of greater stability are described in detail. The 8-ray spectra 
of rubidium®, strontium®’, ekatantalum, and protactinium were measured. The use of thin 
collodion for the backing of sources and windows of counter permitted accurate measurements 
down to 50 kev. The upper energy end points were found to be: rubidium*®—1.60+.03 Mev; ~ 
strontium*®*—1.32+0.03 Mev; ekatantalum—230 kev. The Fermi and K-U plots of Rb® and 


are given. 


NEW solenoidal type 8-ray spectrograph 

has been described recently by Witcher and 
co-workers! which has as its greatest advantage a 
good combination of resolving power and effi- 
ciency. The relatively high value of the efficiency 
(1 percent of the total number of particles 
emitted can be detected) has been found ex- 
tremely useful for the study of many artificially 
radioactive substances of small intensities. One of 
the major advantages of this instrument is that 
radioactive sources covering fairly large circular 
areas may be utilized without decreasing the 
resolution seriously. This permits the use of ex- 
tremely thin sources and minimizes dispersion 
due to scattering within the source material. 

As originally constructed the resolution of the 
instrument (half-width of line divided by H7p=5 
percent) was satisfactory for continuous spectra, 
but the complete resolution of many line spectra 
such as ekatantalum required a baffle system of 
greater resolving power. Further work has there- 
fore been done on this instrument to increase its 
usefulness by designing an alternate baffle system 
which could be substituted when stronger sources 
were available and when it was desirable to study 
a line spectrum in more detail. The new baffle 
system, besides reducing the spread of the beam, 
introduces some improvements in the geometry 
of the spectrograph, i.e., it enables the source, 
slits, and counter to be more accurately lined up 
and centered. 

This type of spectrograph requires a special 


* Publication assisted by the Ernest Kempton Adams 
Fund for Physical Research of Columbia University. 

1 Witcher, O’Conor, Haggstrom, and Dunning, Phys. 
Rev. 55, 1135A (1939); Witcher, Phys. Rev. 60, 32 (1941). 


form of Geiger-Mueller counter which, in general, 
does not have the characteristics which make the 
ordinary counter a satisfactory detector for 
6-rays. Further work has therefore been done to 
improve the counter characteristics and to make 
it useful for work at low energies. 


THE BAFFLE SYSTEM AND THE RESOLUTION OF 
THE SPECTROGRAPH 


In this spectrograph a homogeneous magnetic 
field is used to focus the electrons. The field is 
produced by a solenoid and that portion of the 
solenoidal field is utilized which is homogeneous 
to within 1 percent. The path of an electron with 
an initial momentum making an angle a with the 
magnetic field is a helix of radius (mcv/elH) sin a 
and pitch (27mcv/elH1) cos a, whose axis is parallel 
to the magnetic field. An electron emerging from 
a point source on the solenoidal axis will return to 
the axis after completing one revolution, the dis- 
tance between the source and the point where it 
returns being equal to the pitch of the helix. If z 
denotes the displacement of the electron in the 
direction of the axis (the source taken as origin) 
and R the displacement of the electron from the 
axis in a direction perpendicular to it, then 
R=(2mcv/el1)-sin a-sin cos a. If we 
write D=2mcv/eH, then 


R=D-sin a-sin (z/D cos a). 


For fixed z; and D, (at a plane distant 2, from 
the source and for a fixed ratio of electron velocity 
to field strength), R will vary with the angle of 
emission a, and the mathematical condition that 
will give a point of focus at R=R,, is that 
dR/da=0. The symmetry of the field about the 
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CORRECTION COILS 
COUNTER FILLING SYSTEM 


MAIN WINDINGS 


-+—TO PUMP | 
= 

/—aIR LOCK 
— SHUTTER CONTROL 

SOURCE HOLDER COUNTER SHIELD 

BAFFLE SYSTEM ORS RR 


COUNTER LEADS 


Fic. 1. Solenoidal 8-ray spectrograph. 


axis then leads to the focusing of all electrons 
with the same velocity and an angle of emergence 
a in a ring of radius Rj, in the plane z=2z;. A 
point focus is obtained if R, vanishes concurrently 
with dR/da. This requires that a be equal to 0, 
and since in this case there is no separation of the 
different velocities the focal ring must be of finite 
radius. 

A diaphragm situated at a distance z; from the 
source and having an annular slit of radius R, 
will permit the focused electrons to pass through 
the instrument. A system consisting of a source 
and one such diaphragm would allow other 
electrons besides those which are focused at the 
slit to pass through, and an additional diaphragm 
is necessary to determine completely the helix 
within the limits set by the finite widths of the 
slits, if it can be assumed that no electron makes 
more than one revolution. By having a rod of 
material situated along the axis for most of the 
length of the spectrograph, the possibility of an 
electron making more than one revolution is 
ruled out. 


FIRST BAFFLE SYSTEM 


While only two slits are necessary, the baffle 
system originally designed contained a large 
number of slits which outlined the paths taken by 
the electrons of maximum and minimum velocity 
and angle of emergence that were allowed to pass 
through the spectrograph. An inner series of disks 
which are supported on a rod coincident with the 
axis outlines the curve 7.16 sin (27/90). The brass 
wall of the solenoid at its central point together 
with the inner diameters of two rings on 
either side of the center lies on the sine curve 


12.5 sin (zr/90). This baffle system, though more 
complicated than necessary and especially diffi- 
cult to study exactly when computations of the 
average efficiency of extended sources are being 
made, has been designed to reduce scattering. 
There has been no evidence of any error due to 
scattering when this baffle system has been used. 
It was constructed so that the electrons were 
focused 86 cm distant from the source ; the inner 
radius of the annular slit at that plane was 1 cm, 
and the outer radius was 1.5 cm. The value of D 
for which the efficiency of the spectrograph is a 
maximum is 94.7/x7. The emergent angles of the 
electrons which constitute the beam vary be- 
tween 12° 36’ and 23° 30’. The solid angle at 
maximum efficiency is 1 percent and the half- 
width of the resolution curve (the plot of effi- 
ciency against D) is 5 percent. The choice of this 
geometry which leads to a good combination of 
resolution and efficiency under the practical re- 
strictions of the size and power of the available 
solenoid has been discussed previously by 
Witcher.' He has also discussed the theoretical 
and experimental resolution curves for an ex- 
tended circular source 2 cm in diameter. 


SECOND BAFFLE SYSTEM 


In the new baffle system the number of defining 
slits has been reduced to three. Two additional 
diaphragms whose slits are so large that they do 
not intercept the beam have been introduced to 
eliminate scattering. The annular ring slits are 
carefully centered with respect to the axis by 
means of bars connecting the outside ring to the 
inner disk. The 8-ray source is in the form of a 
circular area 0.5 cm in diameter centered on the 
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axis of the instrument. The Geiger-Mueller 
counter has an open area greater than that of the 
last slit and is situated 1 cm behind the dia- 
phragm where the electrons come to a focus.The 
maximum solid angle is obtained for D=30 and 
for this value of D the angles a range from 8.9° to 
18.7°. The half-value of the resolution curve 
divided by the Hp value at the peak is 2.5 percent. 
Figure 1 shows the spectrograph with the new 
baffle system in place. 


5 CM 


Fic. 2. G-M counter for use with spectrograph. 


SOURCE 


While the geometry of the instrument de- 
termines the efficiency for any point on an ex- 
tended source, the intensity for any given source 
is a function of the diameter of the source. The 
resolution becomes greater as the diameter of the 
source is decreased. Sources up to 2 cm in diame- 
ter may be used. 

The source, which is made by dropping the 
radioactive material in the center of a film of 
collodion of the order of 0.7-u thickness, supported 
by an aluminum ring, is placed in position by 
fitting the ring into a diaphragm which has a 
central hole. Scattering from the collodion film 
support should be negligible for electrons above 
30 kev.? 

To minimize the small residual scattering that 
might take place by multiple scattering from the 
brass wall of the instrument, a sheet of fiber 
board 1/64” thick is rolled to cover the wall of 
the solenoid between the source and the final slit. 


GEIGER-MUELLER COUNTER 


Inasmuch as the solenoidal spectrograph fo- 
cuses the electrons in an annular ring, the elec- 
trons must enter the counter at one end and the 
cylindrical electrode must therefore be at least 3 
cm in diameter. To keep the background of the 


2 A. Flammersfeld, Zeits. f. Physik 112, 727 (1939). 


counter low, the length of the cylinder was made 
relatively small (4 cm) and the central wire was 
covered by glass tubing leaving an active region 
of only 2-cm length exposed. The effect of this 
restricted geometry is that the pulses generated 
are small. Previously, the counter window was 
mounted on a metal grill and it was found neces- 
sary to apply a potential intermediate between 
that on the central wire and the cylindrical 
electrode to prevent spontaneous discharges 
through the counter to the grill. The counter 
characteristics have been improved by substi- 
tuting a plastic grill with 60 percent open area 
(3” holes) which results in increased stability and 
a wider plateau. 

0.5-mil thick Cellophane windows are satis- 
factory for use down to 1900 Hp and approximate 
corrections can be made down to 800 Hp. For 
accurate spectra below this energy, thinner collo- 
dion windows were used to close the counter. The 
film is made by dropping a half-and-half solution 
of collodion and amyl acetate on distilled water. 
The size of the holes in the grill determines the 
strength, i.e., the thickness necessary to with- 
stand the pressure in the counter. In order to 
obtain a film that was tight over the compara- 
tively large area of 10 cm?, a film made of 12 
superimposed layers of original film was neces- 
sary. Since each film is approximately 0.08. 
thick the total thickness is about 1p. 

The collodion film when laid over the open 
area of the grill extends only about 2 mm beyond 
the open area and the pressure in the counter is 
sufficient to seal the thin film to the grill. A thin 
film of stopcock grease provides a satisfactory 
seal between the ground edge of the counter and 
the ground plastic grill, the two being held 
together by means of rubber bands that pass over 
hooks on counter and grill. A counter gas mixture 
of 9 cm of argon and 1 cm of alcohol has been 
found satisfactory, a five-gallon ballast flask 
being connected to the system to stabilize the 
pressure over a long period against possible 
porosity in the thin films. The operating charac- 
teristics of the counter are now very satisfactory 
with the threshold at 800 volts and a plateau of 
more than 150 volts. The counter is shown in 
Fig. 2. 

Such a counter closed with a thin window of 
approximately 1-u thickness of collodion will give 
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accurate spectra for electrons ranging down to 50 
kev without appreciable corrections.*? Since the 
collodion films will vary somewhat in thickness, 
it is not possible to extrapolate accurately to zero 
thickness by placing one, two, and three addi- 
tional thicknesses in the path of the beam by 
means of a shutter. It is therefore necessary to 
determine the region in which the spectra are not 
distorted by absorption in another way. This is 
done by placing one additional film in the path 
of the beam, varying the energy of the beam, and 
determining the point at which the intensity of 
the beam is not reduced by the film. For all 
energies above this point, the spectra are 
accurate. 


COUNTER STANDARDIZATION 


As a check on the sensitivity of the counter, a 
source of uranium oxide electroplated on copper 
is placed in one of the shutter windows. The use 
of this 8-ray source has been found preferable to a 
radium y-ray standard outside the instrument, 
especially when it is necessary to compare the 
operation of the counter over long periods of 
time. It is thus possible to measure the decay of a 
spectrum with accuracy. The counter is com- 
pletely shielded from the uranium when the 
shutter is turned through 180°, and the back- 
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NUMBER OF ELECTRONS DETECTED 


793 1586 2379 
Mp 
Fic. 3. Uncorrected spectrum of s,ekatantalum*, 


ground of the counter, which is about 20 per 


minute, is not affected by the presence of the 
uranium in the spectrograph. 
THE SPECTRA OF EKATANTALUM AND 
PROTACTINIUM 
A. V. Grosse, E. T. Booth, and J. R. Dunning*® 
have reported on the preparation of a 8-emitting 


3A. V.Grosse, E. T. Booth, and J. R. Dunning, Phys. Rev. 
59, 322 (1941); G. T. Seaborg, J. W. Gofman, and J. W. 
Kennedy, Phys. Rev. 59, 321 (1941). 


isotope of the naturally radioactive protactinium. 
Thorium in the form of a solution freed from the 
naturally radioactive protactinium is irradiated 
by slow neutrons from the Columbia cyclotron. 
The irradiated thorium shows a strong 8-activity 
having a half-life of 25 minutes which could be 
explained by the formation of Th?** by resonance 
capture in an (m, y) reaction which disintegrates 
by 6-emission into an isotope of protactinium— 
giekatantalum**, 


NUMBER OF ELECTRONS DETECTED 


793 1586 2379 


Fic. 4. Uncorrected spectrum of »:protactinium™'. 


The chemical separation of ekatantalum from 
the irradiated thorium gives a product which is 
again 8-active, tests on an ionization chamber 
giving the half-life as 27.4+0.3 days. The decay 
of this substance by the emission of an electron 
would give rise to a new nucleus, presumably an 
unstable isotope of uranium, 92U**. To establish 
definitely that the ekatantalum emitted nuclear 
electrons, the 8-ray spectrum was examined. 

The source after preparation by Dr. Grosse 
was placed over a 2-cm diameter area on a thin 
collodion film about 0.74% thick. The layer of 
material was very thin so as to produce no effect 
from back scattering above the energy at which 
the accuracy of the instrument is limited by ab- 
sorption in the Cellophane counter window. The 
comparison of the spectra obtained from two 
sources whose thickness and backing were differ- 
ent, established that errors due to back scattering 
were not apparent in the region of the spectrum 
reported on. The counter window was of 0.5-mil 
Cellophane. 

The spectrum of ekatantalum shown in Fig. 5 
has been obtained from the experimental results 
plotted in Fig. 3 by correcting for decay, division 
of ordinates by // to allow for the varying resolu- 
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tion of the instrument, and by correcting for 
absorption in the Cellophane window. This last is 
done by placing one, two, and three additional 
layers of Cellophane in front of the counter 
window and extrapolating to zero thickness. The 
correction curve as a function of the energy is 
shown in Fig. 6 and gives the factor by which the 
actual intensity must be multiplied to yield a 
corrected spectrum. These corrections have been 
applied down to Hp=800 although it is difficult 


192 KV 


Fic. 5. Corrected 8-ray spectrum of »,ekatantalum**. 


to establish just where the approximation be- 
comes to poor to apply. 

The spectrum of ekatantalum contains at least 
four y-ray conversion lines at 63 kev, 77.2 kev, 
192 kev, and 293 kev. Two of these lines, those at 
192 kev and 293 kev, appear to be complexes of 
unresolved lines. The spectrum above 230 kev 
seems to consist of conversion lines alone; the 
lines below 230 kev appear to be superimposed on 
a continuous background. The spectrum has been 
followed with an accuracy of 2 percent on the 
peaks and the points have been taken too close 
together to admit of any undiscovered line of any 
intensity in the portion of the continuous back- 
ground between 1000 and 1500 //p. If we rule out 
the possibility of this portion of the curve having 
been formed by a large number of lines all of 
small intensity, it can be concluded that the 
ekatantalum emits nuclear electrons, and that a 
uranium isotope 92U?* is formed. The maximum 
energy end point is masked by the line at 192 kev. 
By subtracting this line from the continuous 
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spectrum, the end point of the continuous spec- 
trum appears to lie in the neighborhood of 
230 kev. 

A source of naturally radioactive protactinium 
91Pa**! prepared by Dr. Grosse was examined in 
order to compare it with its isotope and thereby 
establish independently the purity of the eka- 
tantalum. Although Figs. 3 and 4 (plot of number 
of electrons detected against current or corre- 
sponding Hp) show a similarity in the positions 
of some of the most intense conversion lines, the 
relative intensities and positions of the lines differ 
markedly. The spectrum of protactinium ob- 
tained checks as far as the resolved lines are con- 
cerned with the spectrum obtained by L. 
Meitner.* 

To check further the purity of the sources and 
the modes of decay, both the ekatantalum and 
protactinium were remeasured after 27 days. The 
curve 2 of Fig. 5 shows that all parts of the 
ekatantalum curve decay at the same rate. The 
half-life of the ekatantalum measured on the 
spectrograph by taking an average of the decay 
at various points along the curve was 26.0+0.5 
days which checks well with the observed decay 
period for the bulk material. Since there was no 
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Fic. 6. Correction factor for absorption in 0.5-mil 
Cellophane window. 


observable change in the shape of the ekatanta- 
lum curve, it was not possible to make any 
discovery as to the mode of decay of the product 
uranium nucleus. The natural protactinium 
showed no decay at the second measurement as is 
to be expected if the sample is pure. The spectra 
shown in Figs. 3 and 4 are the experimental plots 
of the electron intensity against the current or 
corresponding Ip. 


4L. Meitner, Zeits. f. Physik 50, 15 (1928). 
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THE SPECTRUM OF RUBIDIUM* 


Helmholz, Pecher, and Stout® have succeeded 
in producing a radioactive product with a period 
of 19.5 days by the bombardment of Sr with 
16-Mev deuterons by the (d, a) reaction. They 
definitely established this product to be rubidium. 
By comparing their results with those of Snell® 
(who found an 18-min. and 18-day period after 
the bombardment of Rb®* and Rb*’ with slow 
neutrons giving rise to the isotopes Rb** and 
Rb**), they concluded that the radioactive Rb 
was Rb**, since this is the only Rb of the two 
that could be formed by the Sr (d, a) reaction, 
Sr®8 being the heaviest Sr. 

Three microcuries of Rb** were prepared in the 
form of Rb Cl for examination in the spectro- 
graph. The solution was dropped onto a collodion 
film about 0.7 thick resulting in an invisible 
layer of Rb Cl when dry. The spectrum was first 
studied with a counter window of 0.4-mil thick 
Pliofilm, and a careful study made of the upper 
energy region. For the low energy region a 
counter window of approximately 1-y thick collo- 
dion was employed. The curves from these two 
measurements were reduced to the same intensity 


INTENSITY 
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21435 42670 
4p 


64305 
Fic. 7. 8-ray spectrum of rubidium™, 


at a point above the energy at which absorption 
occurs in the Pliofilm window. 

Figure 7 shows the corrected spectrum of Rb**. 
It is seen that the upper end point of the experi- 
mental spectrum lies at Hp=7288. After cor- 


5A. C. Helmholz, C. Pecher, P. Stout, Phys. Rev. 59, 
902 (1941). 
® A. H. Snell, Phys. Rev. 52, 1007 (1937). 


recting this upper energy end point for the finite 
spread of the resolution curve the value 1.60 Mev 
+.03 is obtained as the true end point. The value 
obtained by Helmholz, Pecher, and Stout from 
absorption measurements is 1.56 Mev calculated 
from Feather’s rule. 


THE SPECTRUM OF STRONTIUM*® 


Sr** is produced by bombarding Sr** with 
deuterons in the (d, p) reaction. A source of 2.78 
microcuries of this Sr** was prepared for study. 


tt \ 


2248 4496 $6744 
Hp 136 Mev 


Fic. 8. 8-ray spectrum of strontium®’. 


The solution was prepared in the lactate form 
and after placing a drop of the concentrated 
solution on the center of the 0.7-u thick collodion 
film the source after drying was invisible. There 
is little probability that an error in the spectrum 
has occurred due to back scattering from the 
source. 

As in the case of rubidium the spectrum was 
first examined with a 0.4-mil Pliofilm window 
closing the counter and a second run was taken 
at the low energy end with a collodion film of 
approximately 1-y thickness on the counter. The 
source strength was such as to give a counting 
rate of only 144 counts per minute at the peak. 
The counting period was such as to give an 
accuracy of +2 percent on most of the points of 
the curve except at the upper and lower energy 
end of the spectrum. Figure 8 shows the spectrum 
thus obtained. The upper energy end point when 
corrected for the finite spread of the resolution 
curve is 1.32 Mev +0.03. The half-life of this 
substance is 55 days. 
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Fic. 9. Fermi and K-U plots for rubidium®. 


THE SHAPE OF THE SPECTRA 


Since both rubidium* and strontium*® lie on 
the second forbidden Sargent curve, it is not to 
be expected that the Fermi plot for these ele- 
ments be straight lines. To facilitate comparison 
between these and other doubly forbidden spectra 
such as P®, the Fermi and K-U plots of Rb** and 
Sr®® are shown, in Figs. 9 and 10. N represents the 
intensity and F(n, Z) is the Fermi function for 
low Z 


F(n, Z) ; 


Fic. 10. Fermi and K-U plots for strontium®’. 


It is seen that toward the upper energy end 
points of both of these spectra the Fermi plot is 
very nearly a straight line but at low energies it is 
convex toward the axis and lies above the ex- 
tended straight line. 

I wish to express my sincere thanks to Pro- 
fessor J. R. Dunning and Dr. E. T. Booth for 
their constant suggestions and advice during the 
course of this work, and to my co-workers, Dr. 
C. M. Witcher and J. O’Conor, S. J., for their co- 
operation and assistance in the research. I am 
indebted to Dr. A. V. Grosse and Dr. J. M. 
Kenny for their aid in procuring radioactive 
sources. 
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The infra-red bands of SiH,, measured originally by Steward and Nielsen, have been 
remeasured at considerably higher resolving power. From a rotational analysis of the bands 
the equilibrium value of B and values of B for the vibration states v2, v3, v4, ve+vs, and vs+% 
have been determined, B being h/8x*Jc. These have, respectively, the values By = 3.00 cm™, 
B” =2.96 cm™, B’(v3) =2.95 cm™, B,"(v2) = 3.06 cm™, B_'(v2) = 3.02 cm™, B’(v4) = 2.89 
=3.05 cm, =2.97 cm, B_'(vs+v4) =2.94 cm, By,'(v2+ =3.05 
cm™!, and B_'(v2+¥3) =2.98 cm™!. The centers of the bands have been determined to be the 
following: v2=974.6, 2190.6 cm™, »y=914.2 cm™, (v2+ 3) = 3148.8, = 3100.2. The 
value of {4= 4—{£;=0.454 is in good agreement with the value predicted by Jahn for valence 
force fields. From the band centers the valence force constants are calculated. 


I. INTRODUCTION 


HE infra-red spectrum of SiH, was investi- 
gated originally by Steward and Nielsen! 

who observed three principal absorption regions 
near 900 cm~'!, 2200 cm-, and 3100 cm~, all of 
which were partially resolved into rotational 
structure. The first two of these were bands 
exhibiting P, Q, and R branches and were 
identified with the fundamental vibration fre- 
quencies v4 and ps3, respectively. The third region 
was found to consist of two overlapping bands 
with centers at 3095 cm and 3153 cm~ which 
may be identified as v3+»,and 73, respectively. 
These measurements revealed, moreover, that 
the rotational structure, especially in the 11-z 
region, was much more complex than was con- 
sistent with the simple theory of tetrahedrally 
symmetric molecules prevalent at that time. This 
structure has since been studied by Murphy? who 
has given a qualitative description of it on the 
basis of a theory proposed first by Childs and 
Jahn’ and by Jahn‘ and subsequently verified by 
Shaffer, Nielsen, and Thomas,® that because the 
frequency », has nearly the same value as the 
twofold degenerate frequency v2 a resonance 


* Now at Mt. Holyoke College, South Hadley, Massa- 
chusetts. 
** Now at University of Toledo, Toledo, Ohio. 


2G. M. Murphy, J. Chem. Phys. 8, 71 (1940). 

3 W. H. T. Childs and H. A. Jahn, Proc. Roy. Soc. 
A169, 451 (1939). 

*H. A. Jahn, Proc. Roy. Soc. A168, 469 (1938). 

5 W. H. Shaffer, H. H. Nielsen, and L. H. Thomas, 
Phys. Rev. 56, 1052 (1939). 


interaction of the Coriolis type is set up. This 
interaction removes the degeneracy in the rota- 
tional quantum number K so that a rotation level 
is split into (J+1) component levels. The result 
is that a multiplet line should occur where on 
the basis of the simple theory a single line was 
predicted. 

The data of Steward and Nielsen were not 
sufficiently good, however, to yield very accurate 
information concerning the structure of the 
molecule. Since then the experimental techniques 
have been so much improved that almost com- 
plete resolution of the bands might now be hoped 
for. In an effort to obtain data which might yield 
reliable information concerning the intermolecu- 
lar distances and, perhaps, the force constants we 
have undertaken to remeasure these bands. 


II. EXPERIMENTAL 


The spectrometer used to remeasure this spec- 
trum was the same used in similar investigations 
carried out in this laboratory. For these measure- 
ments two echellette gratings were used, one 
ruled with 900 lines per inch for the region 9.04 to 
12u and another ruled with 4800 lines per inch for 
the regions 3.2u and 

Two absorption cells were required; one 6 cm 
long was most useful for the regions near 11.0u 
and 4.6u and another 20 cm long was best suited 
for the work at 3.24. These were of glass tubing 
closed at the ends by polished plates of rocksalt. 

In all cases it was possible to operate the 
spectrometer with slits adjusted to subtend an 
angle equivalent to 0.25 cm~. Deflections of the 
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Fic. 1. 4.6-u region of SiH,. 


galvanometer were recorded with the absorption 
cell in the light beam and with the cell out of the 
light beam so that the data could be recorded in 
terms of percent absorption. Records of the data 
were made at intervals on the spectrometer circle 
equivalent to about 0.2 cm7!. 

The SiH, gas used was taken from the same 
container from which the samples used by 
Steward and Nielsen were obtained. The gas had 
been prepared for them by Professor Warren C. 
Johnson of the University of Chicago. 


The Absorption Band at 4.6u 


The measurements made on the absorption 
region near 4.6u, plotted in percent absorption 
against a scale of frequencies, are shown in Fig. 1. 
The upper curve (A) represents the pattern ob- 
tained by using the 6-cm absorption cell filled 
with SiH, gas to a pressure of 12 cm Hg and the 
lower curve (B) was obtained with the same cell 
filled with SiH, to a pressure of only 6 cm Hg. 
Each curve represents the characteristics which 
were repeated in three independent runs over this 


region. Curve A may be said to reproduce the 
data of Steward and Nielsen with high fidelity 
when one considers that the slit widths and the 
corresponding intervals on the circle at which 
data were recorded were in this case roughly one- 
sixth of those of the earlier attempt. The second 
curve is much better resolved, shows much more 
detail and much sharper lines than the upper 
curve. The actual heights of the principal lines 
are not altered appreciably, but the relative 
intensities of the satellite lines to the principal 
lines are quite different. This suggests that in 
curve A the spectrum may be regarded as much 
“over-exposed” and that curve B presents a far 
more accurate picture of the 4.64 band than does 
A. In the first column of Table I are set down the 
frequency positions of the principal lines and the 
more important satellite lines. 


The Region 9.0u to 12.0u 


Figure 2 represents the details which were re- 
peated in three independent sets of data taken on 
the absorption by SiH, gas between 9.0u and 
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12.0u. The region from 800 cm~ to 950 cm is by 
far the more intense and for the measurements on 
this region a cell 6 cm long filled with SiH, gas to 
a pressure of 5 cm Hg was found suitable. In 
order to make measurements on the absorption 
by SiH, in the region from 950 cm™ on, it was 
necessary to increase the pressure of the gas to 
12 cm Hg. The absorption pattern shown in Fig. 2 
reproduces quite faithfully the earlier measure- 
ments by Steward and Nielsen except that much 
more detail is revealed. In the first column of 
Table II are set down the frequency positions of 
the more prominent peaks measured throughout 
this region. 


The Absorption Bands at 3.2u 


In Fig. 3 are shown the reproducible details of 
three independent sets of measurements made 
over the absorption region near 3.2u in the 
spectrum of SiH,. This region consists of two 
overlapping bands with centers near 3.174 and 
3.23u. The bands are considerably better resolved 
and much more detail is discernible in the present 
set of data, but in general details the pattern 
obtained by S and N for this region is well re- 
produced. To obtain these measurements it was 
necessary to use a cell 20 cm long containing gas 
at a pressure of 40 cm Hg. In Table III are re- 
corded the frequency positions of the more impor- 
tant lines measured in this region of absorption. 


III. DISCUSSION OF THE RESULTS 


The vibration and rotation of the tetrahedrally 
symmetric X Y,type molecule have been studied in 
considerable detail by several authors. Such mole- 
cules may oscillate in four different normal modes 
with the frequencies 1, w2, w3, and ws. The first 
is frequently referred to as the breathing fre- 
quency and corresponds to a symmetrical periodic 
motion of the Y particles to and from the X 
particle. It is non-degenerate and induces zero 
electric moment. The vibration w2 is perhaps best 
visualized as one where the four Y particles move 
approximately in small ellipses about their posi- 
tions of equilibrium. This frequency is twofold 
degenerate and like w; induces zero electric mo- 
ment. The frequency ws is essentially a vibration 
of the four Y particles along the XY bonds in 
such a manner that two Y particles are ap- 


proaching the X particle while the other two are 
receding. This can take place in three independent 
ways so that the frequency ws; is threefold de- 
generate. The frequency w, is essentially a de- 
formation frequency where one XYX angle 
decreases in magnitude while the opposite XYX 
band angle increases in size. Also this type of 


TABLE I. Frequency positions and identifications of 
absorption lines near 4.4. 


. Frequency Identifica- | Frequency Identifica- 
positions tions tions tions 
2091.2 2210.1 R”(3) 

12. R(3) 
| 
R(4 
2104.8 | 9218.4 R(4) 
2107.4 | 2221.0 R'(5) 
2109.0 P(i4) | 2222.5 R'(5) 
2109.9 P(i4) | 22241 R(5) 
2111.3 | 2226.8 R"(6) 
115. RO 
2115.8 } P(13) 2232.2 R"(7) 
2118.4 P’(12) | 2233.4 R'(7) 
2120.0 P’(i2) | 2234.8 R(7) 
2121.1 P(i2) | 2237.1 R"(8) 
2124.1 P'(11) | 2238.5 R’(8) 
2126.1 P'(i1) | 2240.3 R(8) 
2127.2 P(it) | 2242.9 R'(9) 
2130.1 P’(10) | 2243.9 R’(9) 
2131.3 P’(i0) | 2245.2 R(9) 
2132.8 P(io) | 2246.2 
2135.5 P'(9) | 2248.0 
2137.0 P’(9) 2250.2 R(10) 
2138.6 P(9) 2251.2 
2141.4 P'(8) | 2252.6 
2142.9 P’(8) 2253.9 
2144.5 P(8) 2255.1 R(11) 
2147.2 P'(7) 2256.8 
2148.5 P’(7) 2258.2 
2150.2 P(7) 2260.0 R(12) 
2153.1 P'(6) | 2262.8 
2154.9 P’(6) 2264.5 
2156.1 P(6) 2265.8 R(13) 
2158.8 P’'(5) | 2268.1 
2160.1 P’(5) 2269.6 
2161.9 P(5) 2270.8 R(14) 
2165.3 P’(4) 2272.9 
2166.5 P’(4) 2274.8 
2168.0 P(4) 2276.0 R(15) 
2170.6 P’(3) 2277.9 
2172.2 P’(3) 2279.2 
2173.5 P(3) 2281.2 R(16) 
2179.2 P(2) 2282.5 
2185.1 P(1) 2284.5 
2193.2 R’(0) | 2285.7 R(17) 
2194.6 R’(0) 2287.2 
2196.2 R(0) 2289.5 
2199.6 R'(1) | 2290.6 R(18) 
2200.8 R'(1) 2292.0 
2201.3 R(i) 2294.0 
2204.9 R'(2) | 2295.8 R(19) 
2206.0 R’(2) 2300.3 R(20) 
2207.2 R(2) 2305.5 R(21 
2310.2 R(22 
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TABLE II. Frequency positions (cm~) and identifications of absorption lines from 9.0u to 12.0u. 


i F n Lin i i 
Now Ident. | Wine  Freqency Ident. | Frequency Ident. 
8138 38 872.5 927.1 103 983.3 

2 15. 927.5 104 985.9 R(1) 
3 816.2 | 40 } P10) | 75 105 «988.5 (1) 
4 817.9 | 41 876.1 928.6 106 991.7. R-(2) 
5 820.6 | 42 878.9 Pus) 930.0 107 995.0 
6 821.7 | 43 879.2 ‘ 76 930.5}  R,(6) 108 997.7, 
7 823.2 | 44 880.7 931.1 109 1000.0 

8 824.6 | 45 882.6 P,(8) 931.4 110 1001.0 

9 825.5 | 46 884.2 77 932.2} —_R,(7) 111 1002.6 (3) 
10 826.7 | 47 886.0 PA) 932.7 112 1004.1 
827.9 | 48 886.4 933.8 113 1005.2 

12 831.5 | 49 887.1 78 934.2} R,(8) 114 1007.1 

13 834.2 | 50 888.8 934.8 115 1008.9 R,(4) 
14 835.1 | 51 889.2 P.(6) | 79 9365 | 116 «1010.1 RAS 
15 836.9 | 52 891.8 80 938.2 117 1015.7 
16 838.7 | 53 892.5 81 939.1 P_(6) 118 10175 R (6 
17 842.8 | 54 893.5 pus) | 82 9413 119 1021.0 

18 844.4 | 55 895.3 ‘ 83 944.0  P,(6)P_(5)| 120 1022.5  R_(7) R,(6) 
19 846.4 | 56 897.5 | 84 945.1 121 1026.5 

20 847.6 | 57 898.7 ‘ 85 946.6 122 10284 ~—«R, (7) R_(8) 
21 349.2 | 58 900.4 86 949.2 (5) 123 1031.6 

22 850.0 | 59 901.4 | 87 950.3 124 1032.8 

23 851.3 | 60 902.0 ‘ 88 9514 125 1035.0 (8) R_(9) 
24 854.6 | 61 905.1 P.(2) | 89 954.1 P,(4) 126 1039.3 

25 856.2 | 62 907.1 | 9 9865 127 1041.5 

26 857.6 | 63 909.0 ‘ 91 959.1 P,(3) 128 1045.0 

27 859.6 | 64 912.2 92 9618  P-(2) 129 1048.4 

28 862.7 | 65 913.4 93 964.1  P,(2) 130 1054.1 

29 864.7 | 66 914.7 | 94 9668 Pi) 131 1054.7 

30 865.7 | 67 915.5 95 969.5 P,(1) 132 1061.9 

31 866.5 | 68 916.8 96 972.7 133 1065.5 

32 867.6 | 69 917.7 Ri) | 97 973:8 134 «1068.5 

33 868.0 | 70 919.3 98 976.2 135 1071.7 

34 868.4 | 71 920.3 R.(2) | 99 978.2 136 1074.9 

35 369.0 | 72 926.6 100 979.7 R(0) 137 

36 369.9 | 73 923.1 R.(3) | 101 9808 138 1088.6 

37 871.6 | 74 925.4 R.(4) | 102 982.0 R,(0) 139 1091.6 

140 1094.5 


motion may take place in three independent ways 
so that w, is threefold degenerate. The frequencies 
ws and w, will evidently give rise to a variation in 
the electric moment. In zero order of approxima- 
tion the vibration energy of the XY, type of 
molecule will be 


E,/he=¥ (Vitgi/2)wi, (1) 


i=1 
where g; takes the value of the order of the 
degeneracy of wi. 

When all the nuclei are in their positions of 
equilibrium the tetrahedrally symmetric XY, 
type molecule will have three identical moments 
of inertia. For such a molecule, regarding it as a 
rigid rotator, one obtains as the rotational energy 


Er/he=BoJ(J +1), (2) 
where ], being the equilibrium 


value of the moment of inertia of the molecule. In 
higher orders of approximation the vibration- 
rotation energy is found to be 


T=E/hc=Ty+B(V, J, K)J(J+1) 
+C(V, J, 
(3) 


In (3) Ty is the oscillational term value, 
B(V, J, K)=h/8rI(V, J, K)c, I(V, J, K) being 
the effective moment of inertia for a given 
vibration-rotation state, C(V, J, K) isa constant 
for any one vibration-rotation state and Dy, 
D;,.x, and Dx are the centrifugal distortion 
coefficients. The term ¢B(0) is the energy contri- 
bution of the Coriolis interaction originating with 
the degenerate states Viwi+ Vsws+ Viaws. When 
V3= Vs=0 the coefficient is always zero, but when 
either V3; or V4 is equal to one, ¢ takes the three 
values €,=2Jf3,4, — 23,4, €-= —2(J+1)f 3.4 


where ¢; and £4 are constants depending upon the 
nature of the normal frequencies w3 and w,. It is 
readily shown that {3+{4=3. 

The constants B(V, K, L) and C(V, J, K) de- 
pend in a complicated manner upon J and K, but 
the dependence is always of the order (J XK) to 
the zero power. It is, therefore, probably legiti- 
mate to neglect entirely their dependence on J 
and K and consider them as constants for any 
one vibration state. Moreover, these constants 
contain terms which have as denominators 
(w2—w3) and (w2—w,). These originate with the 
Coriolis forces between w2 and the frequencies ws 
and w,. In the case of such molecules as CH,, 
SiH,, and GeH, the potential function is such 
that w2 and w, are nearly the same. This paves the 
way for a resonance interaction which manifests 
itself in that as the molecule makes a transition 
from the normal state to the states we and w, the 
B(V, J, K) and C(V, J, K) values increase and 
decrease, respectively, above and below their 
values when the molecule is in the normal vibra- 
tion state. Except for those states involving V2we 
and V4w, the constant C(V, J, K) is so small that 
it may completely be neglected. 

The selection rules for the rotation quantum 
numbers for molecules of this type are AJ = +1, 0. 
In the case of the vibration levels Viwit+ Vow2+ws 
and Viwit Vewe+ws the threefold degeneracies 
are removed by the Coriolis perturbation and the 
levels split up into three components. In these 
cases the selection rule for J becomes AJ=1, 0, 
or —1 depending upon whether the transition is 
from the normal state to the +, 0, or — com- 

. ponent levels in ws or ws. For combination states 
between w3; and w, the situation is considerably 
more complicated and we shall consider only one 
case, namely, w3-++w4, in this paper. We defer this 
discussion to a somewhat later point. 

With the aid of the above selection rules and 
the energy relations (3) one may derive the 
following combination relations if the term con- 
taining C(V) and the centrifugal stretching terms 
are neglected : 


R(J—1)—P(J+1) =2(2J+1)(B”—¢:Bo), (4a) 
R(J) — P(J) =2(2J+1)(B’—§:Bo). (4b) 
These relations may be taken to hold quite 


generally for all bands except for combination 
bands between w; and w,. In the above ¢; is the 
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amplitude of the internal angular momentum 
associated with the frequency w;; B’’ and B’ are 
h/8n*c times the reciprocals of the moments of 


inertia of the molecule in the normal state and in 


the excited state, respectively. For the frequency 
w; which is non-degenerate and for w2 which is 
twofold degenerate it is readily verified that the 
¢’s are identically zero. Neither of the quantities 
¢3 nor ¢4 can in general be equal to zero, but must 
satisfy the relation 3. The relations (4a) 
and (4b) differ from the usual combination rela- 
tions only by the presence of the term ¢;Bo. They 
do not, except in cases where ¢;=0, enable one to 
evaluate the B” or B’ themselves unless the value 
of ¢; can be determined by some other method. 
If, however, a band occurs in the spectrum for 
which ¢ is zero the value B’ can be determined. 


TABLE III. Frequency positions and identifications of 
absorption lines near 


Frequency Identifica- Frequency Identifica- 
position tions position tions 
2958.5 P_(15) 3114.0 P_'(6) 
2968.5 P_(14) 3115.3 Ro(2) 
2978.4 . P_(13) 3115.5 R_(6) 
2990.5 P_(12) 3116.8 P.,'(6) 
2999.3 P_(11) 3118.6 R,(1) 
3006.1 P_(10) 3119.4 R_(7) 
3015.9 P_(9) 3119.6 P’(5) 
3026.7 P_(8) 3120.4 R_(8) Ro(3) 
3033.4 P(11) 3121.0 P,'(5) 
3035.7 P_(7) 3124.2 R_(9) 
3041.0 P,(10) 3124.5 P_'(4) 
3045.5 P_(6) Po(9) 3127.0 P.,'(4) 
3049.0 P.(17) 3127.5 Ro(4) 
3050.6 P (8) 3128.0 R,(2) 
3053.8 P..(15) P_(5) 3130.7 P_’(3) 
3056.7 P(7) P.(14) 3133.0 P,'(3) 
3058.6 P.(13) 3135.0 Ro(5) 
3061.4 P,(6) P(12) 3136.0 P_*(2) 
3063.3 P_(4) 3137.0 R,(3) 
3065.4 P (11) 3138.7 P,‘(2) 
3068.2 P(5) P(10) 3140.9 Ro(6) 
3070.2 P (9) 3142.0 P."(i) 
3073.1 P_(3) Po(4) P.(8) 3144.4 P,'(1) 
3076.0 3145.5 Ro(7) 
3079.8 P.,(6) Po(3) 3153.0 R_'(0) 
3081.6 (?) P_(2) 3155.7 R,'(0) 
3083.7 P.,(5) 3159.0 R_'(1) 
3086.6  Po(2) P.(4) 3161.8 R,’(1) 
3088.6 P (3) (2) 3164.4 R_’(2) 
3091.6 P_(1) P,(2) 3167.8 R.'(2) 
3094.3 P.(1) 3170.5 R_'(3) 
3100.9 R_(0)° 3173.7 R,'(3) 
3102.7 Ro(0) R_(1) 3175.6 R_"(4) 
3105.2 R_(2) 3179.0 R,'(4) 
3108.0 P_'(7) 3182.0 R_’(5) 
3108.4 R_(3) 3185.0 R.'(5) 
3109.4 R,(0) Ro(1) 3187.5 R_'(6) 
3110.6 R_(4) 3192.5 R,'(6) 
3111.4 P,'(7) 3193.0 R_'(7) 
3112.7 R_(5) 3198.7 R,'(7) 
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Since B”’ is characteristic of the normal state its 
value must be the same when determined from 
other bands where ¢;#0. This principle may be 
employed to determine the value of ¢; for other 
bands when the quantities (B’”’ —¢;Bo) have been 
obtained from experiment. 

The regions which have been measured in the 
spectrum of SiH, have been interpreted in the 
following manner. The region near 4.64 (2190 
cm) which is a region of intense absorption is 
identified as the band v3 and the very intense 
band at 11.0u (914 cm~") is taken to be v4. The 
absorption band near 10u (976 cm~") is much less 
intense than the band v4 which it partly overlaps. 
This band is identified with the vibration v2. This 
frequency should, according to what has already 
been said, be optically inactive in the infra-red. 
There exists, however, a resonance Coriolis inter- 
action between these two nearly degenerate 
oscillations and the rotation of the molecule 
which is large enough so that the wave functions 
of the two vibration states become sufficiently 
“mixed” that the frequency v2 may borrow 
enough optical activity from v4 to become ob- 
servable in the spectrum. There remains to con- 
sider the identity of the two overlapping bands 
which make up the region of absorption near 3.2. 
Of these the band of higher frequency (3150 cm~) 
is probably v2+v;3 while the other (3097 cm~") is 
undoubtedly »v3+ v4. The reasons for this choice 
we shall have occasion to consider subsequently 
with more care. Each of the bands referred to 
above shall be considered in turn and an interpre- 
tation of them follows hereafter. 


The Frequency v; 


This band consists of a set of quite sharp and 
well-defined lines lying on each side of the central 
line or Q branch. Each intense line is accompanied 
by two less intense peaks also observed by 
Steward and Nielsen. These satellite peaks we 
believe to be the lines of the P and R branches 
of the bands »; from the spectra of the molecules 
Si®®H, and Si®®H,. From Fig. 1, curve B, it is 
evident that they are much less intense than the 
principal lines, a fact which is consistent with 
the relative abundance of the isotopes of silicon. 
The identification of the lines in this band 
is given in Table I. With the aid of Eggs. 
(4a) and (4b) we obtain from these sets of lines 


the values for (B’’—¢£3By) and (B’—¢3Bo) tabu- 
lated in Table IV. 


The Frequency v, 


Because of the close proximity of this band to 
the frequency v2 it is not justifiable to neglect 
C(V) in the expression (3) for the vibration- 
rotation energy. The effect of this term is to 
remove the degeneracy in the quantum number 
K and thus split each level up into J+1 com- 
ponents. This means, of course, that the combi- 
nation relations (4a) and (4b) are not correct as 
they stand for this case, but should contain other 
terms also. The splitting of the rotation levels due 
to the term K?C(V) in the energy is not extremely 
large, however, particularly for small values of J 
so that the effect upon the spectrum will be to 
break up what otherwise would have been a single 
line into a group of components frequently too 
close together to be resolved experimentally. We 
have attempted to identify the lines in the P and 
R branches in »4, taking a line to mean the center 
of what in some cases becomes a fairly broad 
band of unresolved component lines. This has 
been fairly satisfactory, although it will subse- 
quently be highly desirable to make a detailed 
study of the structure of this band from the 
standpoint of the K splitting. This we hope to do 
at a later date. In this manner we have accounted 
for about nine lines on either side of the center. 
For lines corresponding to higher J values the 
structure becomes too complex to proceed in this 
manner.® On the high frequency side the struc- 
ture of v4 is further obscured by the overlapping 
with ve. Our identification of the lines is given in 
the third column of Table II and from the combi- 
nation relations we have arrived at values for 
—§4Bo) and (B’ (v4) —§4Bo). These are found 
in Table V. 


The Frequency v2 


As we have already seen, a Coriolis resonance 
interaction exists between the two frequencies w, 


®It is suggested that the rather complex region near 
875 cm is due to the Q branch of the first upper sta: 
band. It will be seen to resemble the Q branch of the 
fundamental to a considerable degree. No attempt has 
been made to identify any of the lines in the P and R 
branches of this band since they must be weak and can 
scarcely be distinguished from the complex fine structure 
lines. 
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TABLE IV. Rotational constants for SiH, molecules 
derived from the band 73. 


TABLE V. Rotational constants of SiH, for the states 
V2, vot v3, and V4. 


Si@9H, Si®29H,4 SiGH, 
(B” — ¢3Bo) 2.82 cm™ 2.83 2.84 cm™! 
(B’— ¢3Bo) 2.81 2.82 2.82 


vot v2 (va+vs)* vs 


2.96cm™ 2.96cm™! 2.86cm-! 2.82cm™ 1.60 cm-! 
(B’-fiBo) 3.06cm™ 3.02 cm™ 2.91 2.84cm™! 1.53 


and we. Here as in the preceding case we may 
expect the rotation lines to be split or broadened 
because here also the quantity C(v2) is large 
enough to be of some importance. It is readily 
shown, moreover, that for the frequency w2 the 
internal angular momentum is equal to zero; i.e., 
¢2=0. We should, therefore, expect the interval 
between two rotation lines to be about Av=2B. 
We take, as before, a rotation line to mean the 
center of a group of unresolved component lines. 
From the data on the bands v3 and v4 we can 
estimate this interval to be about 6 cm~. Two 
sets of lines, apparently independent of each 
other and of about the same intensities, can be 
found in the band v2. When the combination rela- 
tions are applied to the lines in the two sets the 
same value of B” is obtained, indicating that both 
sets of lines originate in the same lower state. The 
values of the B’ are slightly different. The expla- 
nation for this, we believe, is found in the papers 
by Shaffer, Nielsen, and Thomas. Because of the 
tetrahedral symmetry of the molecule a term is 
permitted in the anharmonic part of the potential 
energy which removes the degeneracy in the fre- 
quency w: so that it will split into w2+ and w2-, 
say. From their discussion it is also clear that the 
effective moments of inertia associated with these 
component states will be slightly different. The 
principal lines in vz have been accounted for on 
this basis up to about the ninth line on either side 
of the center. Beyond this the K splitting ob- 
scures the identification. Our identification is 
given in the third column of Table II and the 
values obtained for (B’’—{2Bo) and (B’—£2Bo) 
(where, however, {2=0) are given in Table V. 


The Frequency v.+v; 


The band whose center lies near 3160 cm™ 
must almost certainly be v1: or v2+v3. The 
frequency v2 is known from Raman data to be 
2187.0 cm.’ The frequency v4, we have seen, is 


7F. B. Stitt and D. M. Yost, J. Chem. Phys. 4, 82 
(1936). 


about 914 cm. The combination frequency 
vitv, should then be about 3100 cm which is 
nearly 60 cm less than the measured value. On 
the other hand v2 and »; have, respectively, been 
found to be about 975 cm and 2191 cm~. Their 
combined frequency would then be about 3166 
cm which is in satisfactory agreement with ex- 
periment. If this band is to be explained in this 
manner we should be able to find, as in the case of 
v2, two sets of independent rotation lines, one set 
associated with what has been termed (v2+v3)* 
and another with (v2+v3)~. This does, indeed, 
appear to be the case and on this basis we have 
been able to identify about seven lines on either 
side of the center of each component band. 
Beyond this, overlapping with the other band 
made the identification impossible. The identifi- 
cations are indicated in the third column of 
Table III and in Table V may be found the values 
of (B’’—¢Bo) and (B’—¢Bo) obtained when the 
combination relations are applied to the band. It 
is obvious that the ¢ appropriate for (v2+~y3) is 
equal to ¢3. 


The Frequency v;+v, 


The first-order Coriolis interaction causes the 
level v3+ v4 to split up into nine component levels. 
The amount of the splitting depends in a compli- 
cated manner upon the values of £3 and ¢4. In our 
analysis of the data on v3+¥4 we shall anticipate 
the values of ¢3; and {4 which may be obtained 
from the information which goes before. It will be 
shown in the next section that ¢;=0.046 and that 
¢4=0.454. With these values we may calculate 
the nine values of ¢€ required in (3) to give the 
energy values of the molecule for a given value of 
J. Because ¢3 is so nearly equal to zero the 
pattern becomes that of three closely spaced 
triplet components. 

For the sake of simplicity we have taken ¢3 
exactly equal to zero and {4 exactly equal to }. 
This is equivalent to letting each of the three 
triplet components for a given J value degenerate 
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into a single level. There will then be three 
different values of €;, €,, €o, and e_, each of which 
is repeated three times. These are e,=2J¢, 
—2¢, —2(J+1)¢. While the e’s are iden- 
tical with those for the state v4, the selection rules 
for J wili be more complicated. When the dis- 
cussion of Shaffer, Nielsen, and Thomas is applied 
to this case one finds, in fact, that transitions of 
the type AJ = +1, 0 may take place to each of the 
above three component levels €,, €9, and e_ of any 
value J. These selection rules yield, in effect, 
three different bands, (v3+yv4)4, (vs+v4)o, and 
(vst+v4)_, which will overlap with each other. 
With the aid of the above selection rules we 
obtain the following relations for the positions of 
the lines in the P, Q, and R branches of each 
band: 


P(J)(+0_) = (Ty — 2¢Bo) +J(J—1)B’ 
—J(J+1)B"+2a(J)Bo, (5) 


where a(J) takes, respectively, the values Jf, 0, 
and —(J—1)¢; 


Q(J)(+0_) = (Ty — 2¢Bo) 
+J(J+1)(B’—B”)+2a(J+1)Bo, (6) 
and 


R(J)(*+0_) = Ty — 2¢Bo) + (J+1)(J+2)B’ 
—J(J+1)B"+2a(J+2)Bo. (7) 


These relations will in turn yield the following 
combination relations: 


R(J—1)(*0_) — P(J+1)(+0_) 
=2(2J+1)B”, (8a) 
R(J)(*+0_) — P(J)(*0_) 
= (8b) 


The band »3+»,4 overlaps rather badly with the 
band »2+¥73, but the low frequency side is quite 
unobstructed. From the relation (5) it is clear 
that the lines of the component band (v3+»,)_ will 
spread to lower frequencies than the lines of the 
other two component bands. In fact it will be 
seen that two lines in the P branch of this com- 
ponent band will be separated by an interval of 
about 2(B’+¢Bo) which is roughly 3/2 the 
normal rotation spacing. This appears consistent 
with the frequency intervals of the most promi- 
nent lines in the low frequency side of this band. 
It has been possible to account for most of the 
prominent peaks in »3+ 4 on this basis and their 


identification is given in Table III in the third 
column. The values of the rotational constants 
from this band from the combination relations 
are set down in Table VI. No great accuracy for 
these values can be claimed, however, because of 
the approximation which was made that {,=¢ 
=0.5 exactly. 


IV 


In the tables of Section III are summarized the 
data arrived at from a remeasurement of the 
bands in the spectrum of SiH,. From our analysis 
of these data it has been possible to arrive at 
values for the rotational constants (B’’—¢,Bo) 
and (B’—¢;Bo) characteristic of the normal state 
and the excited vibration states, respectively. 
From these values and from the knowledge that 
€s+¢4=0.50 and that ¢2=0 it is possible to arrive 
at values for the quantities By, B’, and B’(»,), as 
well as the values of £3; and ¢4. The values ob- 
tained are set down in Table VI. 

From the value of Bo given in the table above 
it is possible to compute the equilibrium value of 
the moment of inertia J) of the SiH, molecule. 
Taking h to be 6.624 X 10-*’ one obtains for J» the 
value 9.33 X10-*” g cm?. From this value it is in 
turn possible to determine the equilibrium values 
of the distances between the silicon atom and the 
hydrogen atoms and the distances between two 
hydrogen atoms. It is readily verified that 
Io= (8/3)mu(rsi_n)?=mu(ru_n)?. This leads to 
the values (rsi-n)=1.55X10-* cm and (ry_n) 
= 2.58 10-* cm. 

It is of interest to note in connection with the 
values of the rotational constants that the value 
of [B” —(B’(v2))» ]= —0.08 where we mean 
by (B’(v2))w the average value of B,’(v2) and 
B_'(v2), is very nearly equal to the negative of 
[ — B’(v,) ]}=0.07 cm~. This suggests that the 
variation of B’(v2) and B’(v4) from B” is essen- 
tially due to the Coriolis interaction between the 
two frequencies and that the other contributions 
are small. It is further of interest to note that the 
values of the ¢’s obtained from the data are in 
good agreement with those predicted by Jahn on 
the basis of valence force fields. 

With a knowledge of the rotational constants 
at hand it is possible at once to infer what are the 
centers of the bands. This information is gathered 
together in Table VII. 
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TABLE VI. Rotational constants of the SiH, molecule. 


Bo =3.00 B’(v3) =2.95 =3.05 
B” =2.96 B’(v4) =2.89 Bo'(v3+v4) =2.97 
=3.06 =3.05 B_’(v3 +4) =2.94 cm=! 
B'(v2) =3.02 B_’(v2+v3) =2.98 {4 =} =0.454 


TABLE VII. Positions of the centers of the vibration- 
rotation bands in SiH,. 


v2 ¥3 “4 


974.6cm™ 2190.6cm™ 914.2cm™! 3148.8 3100.2 


The potential energy functions of the tetra- 
hedrally symmetric XY, have been discussed by 
several authors.’ The most general form may be 


written 


3 3 3 
V=}{hki 2k nikiths ni? 


3 3 
+k, coi}, (9) 


where the k’s are complicated relations involving 
nine arbitrary and independent potential energy 
constants K,---Ky. We do not have available 
enough data to evaluate all of these. It is possible, 
however, to evaluate the five constants k,- --ks 
since five independent pieces of information are 
available. Four of these are supplied by the values 
of the normal frequencies and the fifth is avail- 
able when we make use of the fact that £3 and ¢4 
are functions also of the potential energy con- 
stants. The values of the centers of the funda- 
mental bands are, of course, different from the 
actual values of the fundamental normal fre- 
quencies because of the anharmonicity. For this 
effect we have no means at our disposal to make 
corrections. We may, however, to a fair degree of 
approximation simply use for the normal fre- 
quencies the centers of the fundamental bands. 
The values of the constants k;- - -ks arrived at on 
this basis are the following: k:=1.177X105, 

8D. M. Dennison, Astrophys. J. 62, 84 (1925); J. E. 


Rosenthal, Phys. Rev. 45, 538 (1934); W. H. Shaffer, 
H. H. Nielsen, and L. H. Thomas, reference 5. 


1.4891 X 105, k3 = 3.963 10°, kg = 1.315 X 105, 
ks=0.754X105 dyne cm. It is of interest also to 
evaluate the five constants which will occur in a 
generalized valence force potential energy func- 
tion. Such a generalized valence potential energy 
function will take the form: 


4 


i=1 
+Koro ba;;(67;+6r;) 
4 
+2K, bax jx} ’ (10) 


i=1 

j= 

where 6r; is the change in the distance between 
the central atom and the ith hydrogen atom and 
da;; is the variation in the H;—Si—H, valence 
band angle. In an unpublished work Shaffer® has 
shown that the following relationships hold be- 


tween the k’s and the K’s. 


Ki Rit2kot+kst (ki t2ks) ], 
|, 
|, (11) 
], 
]. 


From these we may evaluate the K’s and their 
values are given below: 


K,=2.735 X10° dyne cm, 

K.=0.2262 X 10° dyne cm, 
K;=0.0283 X 10° dyne cm, 
K,=0.0295 X 10° dyne cm, 
K;=0.0196 X 10° dyne cm. 


We observe that the constants K3, Ky, and Ksare 
all of the order of 75 of Ko; i.e., they are of a 
smaller order of magnitude than either K, or Kz. 
This fact is consistent with the point of view that 
the valence force field concept represents a legiti- 
mate approximation to the actual force fields in 
the SiH, molecule. 


®W. H. Shaffer, unpublished work. We are grateful to 
Dr. Shaffer for making this work available to us. 
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Two of the infra-red bands of SiH,, originally measured by Steward and Nielsen, have been 
remeasured at much increased resolving power. In addition a new band near 13.0u has been 
located and measured. The bands have been identified as follows: their centers having been 
determined to be v2=930.9 cm™, v3=2113.6 cm™, and »,=819.3 cm™. From a rotational 
analysis of these bands the equilibrium value of B and the values of B for the vibration states 
v2, v3, and vy have been determined, B being h/8x*Jc. These have, respectively, the values 
Bo=2.90 cm™, B” =2.87 cm™, B’(v2) = 2.96 cm™, B’(vs) = 2.79 cm™. The values of ¢ are found 
to be the following: {4= }—{3;=0.476. From the above data four independent relations existing 
between the generalized valence force constant are computed. 


I. INTRODUCTION 


N a recent paper’ we have reported on our 
remeasurement of the infra-red absorption 
bands due to SiH,. The results were so gratifying 
that it was decided to remeasure also the absorp- 
tion bands in the spectrum of the GeH, molecule. 
Such an investigation would moreover present an 
opportunity to hunt for the forbidden vibration 
transition v2 which by analogy to the SiH, case 
might be expected to be intense enough to be 
observed. 

Such an investigation could easily be carried 
out. No major adjustments of the spectrometer 
were required inasmuch as the absorption bands 
in the spectrum of GeH, lie in the same general 
regions as do the bands in the SiH, spectrum. 
The same diffraction gratings were used, re- 
spectively, as in the preceding case; the slits 
were maintained at the same sizes; in fact 
experimentally all conditions remained the same. 
The absorption cells containing the GeH, gas 
were simply made to replace the cells which con- 
tained the SiH, gas. Nothing further need there- 
fore be said concerning the experimental pro- 


cedure. 
II. THE ABSORPTION BANDS OF GeH, 
The 4.5u Band 


For measurements on this region a cell 6 cm 
long filled with GeH, gas to a pressure of 10 cm 


* Now at the University of Toledo, Toledo, Ohio. 

** Now at Mt. Holyoke College, South Hadley, Massa- 
chusetts. 

1C. H. Tindal, J. W. Straley, and H. H. Nielsen, 
Phys. Rev. (this issue). 


Hg was used. The data are represented by the 
upper curve in Fig. 1. In it one can readily detect 
the peaks referred to by Lee and Sutherland? and 
which also were observed by Steward and 
Nielsen,* but with much less relative intensities. 
When, however, enough of the gas in the cell is 
removed so that the pressure is reduced to 6 cm 
Hg the data are represented by the lower curve 
which is quite consistent with the earlier data of 
S and N when one takes account of the fact that 
for the present curve, galvanometer deflections 
were: taken at intervals of about § of those of 
S and N. We believe the lower curve to be a far 
more accurate picture of the absorption region 
than the upper curve. The principal lines have 
approximately the same heights as before, but 
are much sharper. The subsidiary peaks of L and 
S we believe are in reality quite weak, but can 
be brought out quite strongly when enough gas 
is used so that the principal lines are strongly 
“overexposed.” 

In the first column of Table I are given the 
frequency positions of the more important lines 
in the band and their identifications are set down 
in the second column of the same table. The 
identifications have been made particularly with 
reference to the lower curve in Fig. 1. We have 
identified this absorption region as due prin- 
cipally to the fundamental vibration v3. For these 
we have 

(BY £3Bo) = 2.80 


2E. Lee and G. B. B. M. Sutherland, Proc. Camb. 
Phil. Soc. 35, 341 (1939). 

3W. B. Steward and H. H. Nielsen, Phys. Rev. 48, 
861 (1935). 
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2000 cm** 2050 cm” 200 

50) 

+50 
2100 cmt "2250 
Fic. 1. 4.7-u region of GeH,. 
TABLE I. Frequency positions and identifications of lines in the GeH, spectrum at 4.7y. 
Freq. Freq. Freq. Freq. Freq. 
(cm~!) Ident Ident. (cm~) Ident. | (cm™!) Ident. | Ident. (cm~) Ident. 
2000.8 P(19) | 2031.5 P(14) | 2061.2 P(9) | 2102.5 P(2) | 2149.0 2188.5 R(13) 
2004.6 2034.0 2062.8 2107.9 P(1) | 2151.4 R(6) | 2191.0 
2005.8 2035.8 2065.0 2116.2 2154.0 2193.6 R(14) 
2007.5 P(18) | 2037.5 P(13) | 2067.3 P(8) | 2119.3 R(O) | 2156.6 R(7) | 2198.7 R(15) 
2011.4 2039.5 2071.3 2121.0 2159.5 2203.2 R(16) 
2012.2 2041.5 2073.5 P(7) | 2122.1 2162.1 R(8) | 2208.0 R(17) 
2014.3 P(17) | 2044.0 P(12) | 2076.5 2124.5 R(1) | 2165.1 2212.5 R(18) 
2016.2 2045.4 2079.1 P(6) | 2130.3 R(2) | 2167.8 R(9) | 2217.8 R(19) 
2017.5 2046.8 2082.8 2132.5 2170.5 2222.6 R(20) 
2019.6 P(16) | 2049.4 P(i1) | 2085.3 P(5) | 2135.6 R(3) | 2173.0 R(i0) | 2227.0 R(21) 
2024.0 2051.2 2090.9 ° P(4) | 2138.0 2175.8 2232.0 R(22) 
2025.7 P(15) | 2053.2 2093.0 2141.0 R(4) | 2178.2 R(i1) | 2236.6 R(23) 
2027.8 2055.3 P(10) | 2096.6 P(3) | 2143.4 2180.5 2241.2 R(24) 
2029.5 2057.0 2146.3 R(5) | 2183.6 R(12) 
2059.2 2186.0 


and 
(B’ (v3) €3Bo) =2.78 cm". 

There are five isotopes of germanium and one 
might expect, therefore, as in the case of SiH, 
to find a band »; for each isotopic molecule. Un- 
fortunately, however, the centers of these 
isotopic lines may be expected to lie so nearly 


superimposed that they cannot be resolved 
Experimentally no real evidence of these separa- 
tions is to be found. The above rotational 
constants must therefore be regarded as average 
values for the several isotopic molecules. 

No explanation for the weaker structure in 
between the principal lines is very obvious. Lee 
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750 775 accent" 


om 


900 


875 
Fic. 2. The region 10u to 14 of GeH,. 
TABLE II. Frequency positions and identifications of lines in absorption bands near 12.0,. 
Freq. Freq. Freq. Freq. Freq. Freq. 
(em~!) Ident. Ident (cm~!) Ident. (cm~!) Ident. | (cm™) Ident. (em™~!) Ident. 
749.5 795.7 P(6) | 837.9 R(8) 890.7 P’(8) | 940.0 975.0 | R'(6) 
752.3 799.4 P(5) | 839.0 893.6 941.5 977.6 
754.2 803.3 P(4) | 840.4 R(9) 894.9 P’(7) | 943.0 R’(1) 981.3 R'(7) 
758.0 805.5 845.0 R(i1) | 898.4 945.2 984.3 
759.8 P(15) | 807.5 P(3) | 847.3 R(12) | 899.0 P’'(6) | 947.0 987.2 R'(8) 
763.5 P(14) | 810.5 P(2) | 849.6 R(13) | 904.3 P'(5) | 949.3 R’(2) 991.0 
766.0 814.1 P(i) | 851.4 R(14) | 909.2 P’(4) | 950.6 993.3 R’(9) 
768.0 P(13) | 819.7 R(0) | 863.0 914.2 P’(3) | 952.0 998.6 R’(10) 
770.5 820.5 865.0 917.5 954.0 1004.7 R’(11) 
772.7 } P(12) 822.0 R(1) | 866.6 919.0 P’(2) | 956.0 R’(3) | 1011.3 R’(12) 
773.6 822.8 869.7 P’(13) | 922.0 958.0 1015.6 
776.8 P(i1) | 824.3 R(2) | 874.0 P’(12) | 924.1 P’(1) | 960.3 1019.3 R’(13) 
781.1 P(10) | 826.4 R(3) | 878.0 P’(11) | 932.5 962.6 R’(4) | 1022.5 
785.2 P(9) 828.9 R(4) | 881.7 P’(10) | 933.5 964.5 1026.8 R’(14) 
788.2 P(8) 831.2 R(5) | 883.2 934.5 969.5 R’(5) | 1028.4 
792.5 P(7) | 833.0 R(6) | 886.1 P’'(9) | 936.3 R'(0) | 971.0 1033.5  R’(15) 
793.1 835.5 R(7) | 888.3 938.3 973.7 1039.5 R’(16) 


and Sutherland have suggested that these lines 
arise in some manner from interactions between 
vibration and rotation. This view receives little 
support from the theory in its present state 
although until the selection rules have been 


investigated to higher orders of approximation 
this question can probably not definitely be 
settled. It is suggested here that these lines might 
just conceivably be due to the overtone band 
3v4. 
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The Region 10u to 


The region reported by Steward and Nielsen 
near 10.54 was first remeasured. The data were 
here obtained by using an absorption cell 6 cm 
long containing GeH, gas at a pressure of 10 cm 
Hg. These measurements repeated the earlier 
data quite satisfactorily, especially on the high 
frequency side, when one takes into account that 
the galvanometer deflections were taken at 
intervals one-sixth as large as in the earlier work. 
It was soon discovered, however, that on the 
long wave-length side, about where the data of 
S and N stopped, a new and much more intense 
band set in. This region extended out to beyond 
14.0u and for measurements on this band enough 
gas was removed from the cell to reduce the 
pressure to 6 cm Hg. The data taken are shown 
graphically in Fig. 2 and the frequency positions 
of the more intense lines throughout the region 
are set down in Table II. 

Steward and Nielsen identified the band at 
10.54 as the fundamental frequency v4 and never 
looked farther for another band. In the light of 
this investigation it becomes clear that their 
identification was wrong and that this band was 
in reality the forbidden band v2 and that the 
much more intense band at 820 cm~ must be the 
fundamental frequency v4. 

The rotation lines will here as in the spectra 
of SiH, and CH, be split into fine structure 
because of the resonance Coriolis interaction 
between v2 and vs. The splitting is, however, not 
very great, especially for the smaller quantum 
numbers. In the same manner as we preceded 
with the SiH, bands, so here also we have taken 
a line to mean a group of component lines too 
closely spaced to be resolved spectroscopically. 
On this basis we have identified a considerable 
number of the rotation lines in the P and R 
branches in both of the bands v2 and v4. Our 
identification of the lines in these two bands is 
set down in the second column of Table II. 

In the case of the SiH, spectrum it appeared 
that the twofold degeneracy of the frequency v2 
was removed. This does not appear to be the case 
here. It does not seem possible to find more than 
one set of lines in the band v2. There are certain 
lines in v2 which remain unaccounted for, but 
they are probably to be ascribed to K splitting 
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of the principal lines. The B’(v) values seem to 
decrease with increasing quantum numbers, par- 
ticularly beyond J =6, but this effect may be due 
to centrifugal stretching. With the aid of the 
combination relations we obtain the following 
rotational constants from the band »4: 


(B” — =1.49 
(B’ (v4) — &4Bo) =1.41 


and from the band v2 we get, remembering that 
f2=0, BY’ =2.87 and B’(v2) =2.96 cm. 


The Combination Frequencies near 3.3u 


It was found that there was not enough GeH, 
gas available to fill a cell long enough to re- 
measure these bands. We refer, therefore, to the 
earlier measurements reported by S and N. We 
alter their identification of the band at 3031.0 
cm~!, however, from 7,4 to but retain 
their identification of the band at 2922.7 cm~! as 
v3+v4. In support of this identification we shall 
point out that for v2+ v3 the appropriate ¢ value 
will be ¢3. This should therefore lead to a spacing 
between rotation lines, particularly near the 
center, of about 5.6 cm~. This value agrees very 
well with the data of Steward and Nielsen. The 
band at 2936 cm™ overlaps badly with that at 
3031 cm~!, but the low frequency side remains 
quite unobscured. Our discussion of the fre- 
quency v3+v,4 for the SiH, molecule is equally 
applicable for the frequency v3+ 4 for the GeH, 
molecule. We should therefore expect to find a 
separation between lines in the low frequency end 
of the P branch of about 8.7 cm~'. This seems in 
excellent agreement vith the spacings between 
the first six lines shown in the diagram of Sand N. 


III. DISCUSSION OF THE RESULTS 


From the rotational constants derived from an 
analysis of the bands ve, v3, and v4 one may, with 
the aid of the condition ¢;+{,=4, obtain the 
values for B’’, .”’(v), and ¢4 chemselves. The 
values obtained are given in Table III. It is 
worthy of note that, here also as in the case of 
SiH, spectrum, one finds [B’’—B’(v4)] is very 
nearly equal to the negative of [B’’—B’(v2) ] 
which indicates that the principal contribution 
to the change in B for these frequencies comes 
from the resonance Coriolis interaction. 
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From the value of Bo we can compute the 
equilibrium value of the moment of inertia Jo. 
We find this to be 9.65 KX 10~*° g-cm?. From J» we 
may then deduce the equilibrium values of the 
H—H distance and the Si—H distance. We find 
these, respectively, to be ra_n=2.41X10-* cm 
and rsi-n = 1.48 X 107-8 cm. 

With a knowledge of the values of the rota- 
tional constants we may infer at once where the 
centers of the aforementioned bands lie. The 


TABLE III. Numerical values for the rotational constants 
of the GeH, molecule. 


Bo B” B’(v2) B’(vs) B’(vs) v4 
2.90 2.87 cm™! 2.96cm™! 2.85 2.79 cm™ 0.024 0.476 


values we find are the following: v2=930.9 cm~, 
v3 = 2113.6 cm™, and 819.3 cm. No data on 
the Raman spectrum of GeH, seem to be avail- 
able. The value of »; is therefore not known. 
There are consequently only four independent 
bits of knowledge from which to deduce the 
values of the potential energy constants. These 
are the values of the band centers v2, v3, v4, and 
the ¢’s. As a consequence we cannot here com- 
pute more than four of the independent relations 
among the constant k; * which would occur in the 


‘See either Shaffer, Nielsen, and Thomas, Phys. Rev- 
56, 895 (1939) or our paper dealing with SiH,. 


general potential energy function for such a 
molecule. If we take the values of the funda- 
mental frequencies w, to be equal, respectively, 
to the values of the band centers »,, these we 
find are: ki; =1.117X10° dyne cm, ke=+1.471 
X105 dyne cm, k3=3.709X10° dyne cm, ks 
=5.116X10°5 dyne cm. From these we may in 
turn calculate four independent relations be- 
tween the constants K, which will be present in 
the generalized valence potential energy func- 
tion. The relations are K,;—$K,4=2.589X10° 
dyne cm; K2=0.191X10° dyne cm; K;=0.0011 
10° dyne cm; K;=0.0104X10° dyne cm. It 
should be noted that similar to the case of the 
SiH, molecule the constants K; and Ks are small 
compared to Kz, and K,—3K,. It seems therefore 
quite likely also that K, (which is the coefficient 
of the cross product terms 


i=1 

will be small compared to these, say of the order 
of 1 percent of K,. The value of K, is probably 
not very much different from K,—4$K,4. We shall 
here assume that K, is of the order of 1 percent 
of K, and then calculate the value of ». For »; 
calculated on this basis we obtain 2090 cm™ 
which seems a quite reasonable value. 
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II. Oscillations in a Split-Anode Magnetron 
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A theoretical study is made of the possible oscillations in split-anode magnetrons. It is 
found that large oscillations can occur in certain frequency ranges, these ranges depending 
upon the number of segments. In Fig. 6 these frequency ranges are indicated with y = 21310/AH 
as a function of the number of segments (2m). The large oscillations occur in the first and 
fourth region as tabulated at the top of the figure. 


1. SPLIT-ANODE MAGNETRON— 
GENERAL REMARKS 


HE behavior of a magnetron operated 
under steady conditions was studied in a 
previous paper! where the conditions of oscilla- 
tion were also discussed for the case of a magne- 
tron with one anode of cylindrical shape. This 
first paper will be referred to under the title 
“‘Magnetron I,”’ and a numeral I will be added 
to its equation numbers to avoid confusion with 
the present paper. 

The aim of the present work is a study of 
oscillation conditions in a split-anode magnetron. 
Such magnetrons are built with an anodic 
structure consisting of an even number of 
shells equally spaced so as to form an almost 
complete cylinder about the filament. The most 
popular type has two half-cylindrical anodes, 
but other models with as many as twelve or 
more anodes have also been successfully used. 
The oscillating circuit is connected by one of its 
terminals to the odd-numbered shells, the second 
terminal being connected to the even-numbered 
shells, as shown schematically in Fig. 1. 

The exact field distribution in such magnetrons 
would be very difficult to compute, but a 
reasonable theoretical attempt is possible under 
the assumption of very small sinusoidal oscilla- 
tions. This is illustrated in Fig. 2. Let us call } 
the radius of the anodic structure and 2m the 
number of anodes. The position of the anodes as 
a function of @ is indicated along the @ axis. 
Let V. be the average constant potential of the 


* Now at Brown University, Providence, Rhode Island. 
1L. Brillouin, Phys. Rev. 60, 385 (1941); also Elec. 
Communication 20, 112 (1941). 


anodes and +U,e! the additional alternating 
potential, the + sign corresponding to the odd 
anodes and the — sign to the even ones. The 
potential distribution on the anode circle is 
represented by the curve in Fig. 2. Such a 
broken curve can easily be expanded in a 
Fourier series. Let a be the angle subtended by 
one anode while 8 is the angle between two 
successive anodes. Then a+8=2/n. The po- 
tential on the radius r=) is given by the series 


8 
V(b) = + U,e*'—— 


1 
Xd sin sin (pn6). podd (1) 


The first term p=1 is the most important and 
the next term p=3 can be made zero by taking 


sin (3/2)n8=0, B=22r/3n, a=n/3n, 


which means that the distance between the 
anodes is twice the width of the anodes. 


Fic. 1. Connections of terminals to oscillator. 
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The calculations will be made under the 
assumption that the alternating term U, is 
very small with an e coefficient the highest 
power of which will be neglected. Furthermore, 
we shall keep only the first term of the Fourier 
expansion. 


V(b) = V.(b) sin (n8), 
= sin (38). 


This is the assumed distribution of potential on 
the circle r=}, but one should keep in mind the 
fact that a magnetron built with 2” anodes may 
eventually work on one of the higher modes 
corresponding to 2’ anodes with n’ = 3n, 5n, ---. 
It is easily seen that these higher modes of 
oscillation will increase in probability if the 
interval 6 between the anodes is very small. 
At the limit B=0, a=2/n, the amplitude of the 
pth term becomes 


(2) 


8 1 4 
lim —— — sin (png) =—. (3) 


rp 


2. FUNDAMENTAL EQUATIONS 


The fundamental equations are Maxwell's 
relations which will be written in cylindrical 
coordinates 7, 6, z, where z is the axis taken 
along the filament. It will be assumed that all 
quantities are independent of z. The current 
density components J,, J», including Maxwell’s 
displacement current, are 


(4) 


v, and ve are the components of the velocities of 
the electrons, p the space charge density, EZ, and 
Es the components of the electric field, and 
E.=0. The coefficients €, wo represent the 
dielectric power and the magnetic permeability 
in vacuum. Electrostatic c.g.s. units will be 
consistently used, making 


€9= 1, (5) 


J,rd@ is the radial current flowing through rd@, 
per unit of length z, while Jedr is the angular 
component of the current through dr, per unit 


z. The total radial current per unit of length of 
zis 


Qr 
in J (6) 


It is well known that the introduction of Max- 
well’s displacement terms secures the following 
relation: 


10 10 
div J =- —(rJ,) +— —(Je) =0. (7) 
r or r 00 


-_ 3 


Fic. 2. Potential distribution on the anode circle. 


This may easily be verified, with the conservation 
equation for electric charge 


10a 10 
dp/dt= —div (pv) = —- —(rpv,)—— —(pve), (8) 
r or r 00 


together with one of Maxwell’s equations: 


10 1a 
div —(rE,) —E,=4mp/eo. (9) 
r or r 00 


Let us call H,, He, H, the components of the 
oscillating magnetic field. The symmetry of the 
problem shows that 7, and H, are zero. Conse- 


quently 
div H=0H,/dz=0, (10) 


which proves that H, is independent of z. 
Maxwell’s equations yield 


0H,/rd0=4nJ,, 
rot H=4nJ (11) 
rot E= —y,.dH/dt, (12) 
/rdr —dE,/1d0 = at. 


This system reduces to three relations: (11.1), 
(11.2), and (12). Equation (7) shows that the 


= 
Vv 
ing 
‘he 
a 
a 
by 
0- 
he at’ = 
€0 
Jo=pve+— —. 
4r dt 


— 


relations (11.1) and (11.2) are compatible as it 
gives 
,/drd0= 0°H ,/d60r. 


Hence we can write 


and we are left with one relation 


dE, 
—— = —4Arpy— (14) 
ror ot 


which is to be used together with (9) to define 
the electric field. It is well known that any field 
E can be represented as the sum of an irrota- 
tional field (grad V) and a field E’ with no 
divergence. 


V 
—-—+E’,, 
or (15) 
10V 
—+E’s, 
r 00 
0(rE’,) OE", 
div E’= + =0. (16) 
ror 


Substituting in Eqs. (9) and (14) we obtain 
100V 4rp 


AV=- —r—+-— —= -—— (17) 
ror or 00? €0 
0(rE’s) OE’, 0 
f rJ.d0. (18) 
ror ot 


The assumption will now be made that the 
second field E’ is very small compared to grad V 
and can be consistently neglected. This will be 
proved later on to be a reasonable hypothesis in 
most cases, provided the dimensions of the 
magnetron are small compared with the wave- 
length A in vacuum. 


(19) 


Some special cases where this assumption would 
lead to error will be discussed. Most of them 
lie outside the actual operation conditions of 
the magnetron. Equation (18) may be reduced 
by the use of (16). One differentiates (18) with 
respect to @ and replaces 0E’,/d@ by its value 
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—0drE’,/dr, which yields 


100 
A(rE’,) =- —r—rE’, 
r or Or 
1 


r 06? 


This equation will be used later in the discussion 
of the order of magnitude of E’, and the condi- 
tions for its negligibility. 

The preceding equations must be combined 
with the equations of motion for the electrons 
[Eqs. I.1) to (1.8) ]. 


(20) 


d*r/d?= —2rwn6, 
(21) 
d(r?6) /dt= 
in which 


wu = —dpuo(e/m)H 
Larmor’s angular velocity. (22) 


The second equation in (21) may also be written 
d e 
—[r?(6—wy) |=—rEp. (23) 
dt m 


According to the assumption about the field, 
only grad V will be used in these equations. 
The derivative of a function f(r, 6, t) is often 
to be taken following the motion of one electron. 
The symbol d/dt will be used on such occasions: 


(24) 


Applying this definition, using Eqs. (4) and (17) 
and neglecting E’, one obtains, since €9= 1 (e.s.u.), 


d 
—(rdV/dr) =rd?V/drat 
dt 
(rd V/dr) +190? V/drde 
r 


= —4arJ,—v,0° (25) 


In a similar way 


dav 
“(—) =a°V 
dt\ a0 


V/rde* 


0 
= V/odr). (26) 
r 


— —" 
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r 
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These two important formulae are the general- 
ization of Eq. (1.41). They are very useful 
because they no longer contain the space charge 
density p. Our fundamental equations thus 
reduce to Eqs. (7), (21), (23), (25), and (26). 


3. SMALL OSCILLATIONS IN A SPLIT-ANODE 
MAGNETRON 


As in ‘‘Magnetron I”’ (Section 6) it is assumed 
that very small oscillations are superimposed on 
the static potential and current distributions. 
The smallness of these oscillations is secured by 
a parameter e, the higher powers of which will 
be consistently neglected. We thus assume 


Vir, 0, t)= V(r) +eV.(r, 8, | 
J Ar, 0, t)=I./2ar+eJa,(r, 6, t), (27) 
Jo(r, 0, 9, t). 


The potential V, and the current J, are those 
defined in ‘‘Magnetron I”’ for the magnetron 
operated under steady conditions. 


6.(r) =wy(1—a?/r’). (28) 


In the new field of force electron trajectories 
are slightly modified. Let us follow one particular 
electron leaving the filament at a time fo; its 
position at time ¢ is 


r-, 9. in the unperturbed motion. 
Oc, t), (29) 
in the perturbed motion. 


In the unperturbed motion the electron describes 
the static trajectory calculated in ‘‘Magnetron 
I,” while in the perturbed motion its position 
at time ¢ is at a distance era, €0, from the unper- 
turbed location. It is necessary to calculate the 
value of quantities such as V at the point where 
the electron is actually located. Let f be such a 
quantity. Then Eq. (27) would read 


f(r, 0, )=f-(r, O+.(r, 4, t). (30) 


If this quantity be measured at the point 
ret+era, 0.+6€8, where the electron is located, 
we must write 


(ret a, t) =f.(re, 6.) 
(31) 


terms in é, e*, --- being consistently neglected. 
The fundamental equations can now be 
written with the expansions (27) to (31). 


Equation (7) first yields a very simple relation 
O(rJar)/Or= 0/08. (32) 

Equation (21) gives 

d*r./dt?+¢d*r,/dt* 


= 
m 


2wu + € + rab.) 
276-64]. (33) 


Terms independent of € correspond to the static 
formula (1.19). In ¢ terms we use (28) for 6, and 
get finally 


m 


m 


Equation (23) is treated in the same way. 
Writing only ¢ terms one finds 


d 

dt 
d a? 
dt Ye m 


Turning to Eq. (25) and applying (31), we get 
from the e terms 


d 
V./ar?+r.0 V,/dr | 
t 


= 4ar.J ar - r.6.0° V,/0r00. (36) 


Equation (26) yields 


d 0 
t 


Ye 
+er.J, «| 


+{ (r-+€Fa) 


X +d v./ar)| (37) 


| 
22) 
| 
(35) 
ns: 
17) 
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Terms independent on ¢ just cancel due to (1.21) 


V./dr) =0, (38) 
r 


and ¢ terms result in 
d 
dt Or. 
V./dr 


0 
(rei V./dr), (39) 
r 


where the last term according to (38) is equal to 


4. TYPE OF ELEMENTARY SOLUTIONS 


The new system of Eqs. (32) to (39) is linear 
with respect to all of the alternating quantities, 
a condition which enables us to superimpose 
elementary solutions and to make use of imagi- 
nary exponentials. Each alternating quantity is 
supposed to depend on @ and ¢ in the following 


way 
(40) 


falr, 0, t) = 
fa(r) is an imaginary expression including the 
phase angle factor. 

A solution of type (40) represents one ele- 
mentary solution, and the general solution of any 
particular problem may involve superposition 
of such elementary solutions in order to satisfy 
boundary conditions. One elementary solution 
corresponds to a sort of wave rotating about 
the filament with angular velocity w/n. In the 
example shown in Fig. 2, the potential distribu- 
tion on the anodic cylinder should be represented 
by Eq. (2). This result will be obtained by the 
superposition of two solutions +m and —n of 
type (40), provided these two solutions have 
the same amplitude factor f, for r=b. These 
two +z solutions will, however, depend on r in 
different ways so that the superposition will 
yield a potential (2) only on the cylinder r=0. 
This will be discussed later on. 

An elementary solution (40) must be studied 
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for 
w>O0, (41) 


One must remember, however, that a simultane- 
ous change of the signs on w and m does not 
affect the exponential in formula (40). Hence 


n positive or negative. 


w,n and —w, —n 


should yield the same results; 


w,—n and should also be equivalent. 
Every elementary solution must satisfy the 
boundary conditions on the filament 


V.(a)=0, for r=a. (42) 
From Eq. (40) one readily gets 
0/dt=iw, 8/00=—in. (43) 


Our equation system contains derivatives d/dt 
taken along the path of one electron, as expressed 
in Eq. (24). On the other hand, making use of 
Eq. (31), we have reduced any quantity f, 
observed at the point where an electron is 
located in the unperturbed motion, to a function 
of the coordinates r,., 6, defining the position of 
the same electron in the unperturbed motion. 
The result is that in order to follow the motion 
of the perturbed electron, we only need 7.6, to 
describe an unperturbed trajectory ; this yields 


=i(w—n6,)+7,0/dr.. (44) 


Calculations appear to be very complicated in 
the general case, so we make some simplifying 
assumptions : 


A. very small filament a=0, 
B. magnetron in the critical state J,=0. 


Conditions A and B were discussed in ‘“‘Magne- 
tron I,’’ Section 3, in which we found that the 


? Referring to Eqs. (29) and (31), one should notice that 
they correspond to the use of Lagrange’s definitions in 
hydrodynamics. The same problem could be discussed on 
a different basis with Euler's definitions. In (29) we 
compare the position of one particular electron in the 
unperturbed motion with the position of the same electron 
in the perturbed motion. Euler’s method consists in 
comparing a certain electron located at r, @ (velocities 
v,, 4) in the unperturbed motion with a different electron 
which, in the perturbed motion happens to come to the 
same position r, 6, but with different velocities v,+eu,, 
ve+eug. Euler's method seems to be the one chosen by 
Blewett and Ramo [Phys. Rev. 57, 635 (1940) ]. With 
such definitions Eq. (31) must be cancelled, while Eq. (44) 
should be completed with additional terms because 
following the motion of one unperturbed electron means 
comparing different perturbed electrons. 


(45) 
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Fic. 3. The roots of Eq. (56). 


electrons build a cloud with constant space 
charge density po and rotating about the filament 
with a constant angular velocity wy. 


po= /27e, 6.=wn, =0. (46) 
Hence Eq. (44) reduces to 
d/dt=ia, a=w—Nwn. (47) 


These simple relations must be introduced into 
the general equations of Section 3. The problem 
to be discussed here is very similar to the one 
treated in ‘“‘Magnetron I,’’ Section 6. There we 
had to discuss the role which some damping 
terms should play in a more complete theory, 
and the same remarks apply here. The use of 
Eqs. (43) and (47) is equivalent to taking account 
only of forced vibrations and neglecting proper 
vibrations of the electronic system, as such 
proper vibrations must be damped and practi- 
cally disappear after a short time. 
Equation (34) yields 


e 
=—-——dV,/dr. (48) 
m 


V.. is given by Eq. (1.23) 


m 
Vo= (49) 
2e 
Hence 
(50) 


Turning now to Eq. (35), one finds 


6.=enV,/mr 2a. (51) 


By using (49) and (36) the following relation is 
obtained : 


m 
2—awy V,/or tAnr.J ar, (52) 


while Eq. (32) results in 
—(rJar) = (53) 
or 


Equations (50) and (52) give 
(54) 


Equation (39) may now be greatly simplified by 
using (47) and (53): 


naV,= 


+ (4x1 /n)r.-—(reJ ar) 


le 


0 2m 
or e 


m 
(55) 
e 


Referring now to (50), (51), and (54), one finds 
a homogeneous differential equation for Vq: 


OV. OV. 


Ar—r +Br 
or or or 


+CV.=0, 


1 
A (56) 
n 


B=2wy'/a’, 
C=n(a—2wy"/a). 


Let us introduce a variable y = 
and new coefficients 


A'=(y—n)?A/wu = 

B’ =(y—n)*B/wn =2, 

¢=(y—n)*—2(y—n). 


A solution may be found by taking 


V.=Kr-, 


=0. (57) 
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We thus obtain two solutions: 


x1 B’ 1 
)- ———+—(B”-4A’C’)}. (58) 
Xe 2A’ 2A’ 

It should be immediately noticed that A’, B’, 
and C’ are functions of y and n, so that the 
solutions x}, x2 will depend only upon these two 
variables. Furthermore we should remember the 
conditions (42) stating that V, and dV,/dr have 
to be zero on the filament. This results in the 


2 


Fic. 4. ¢ as a function of y—n. 


necessity of choosing the root with a real part 


greater than one. 
Re(x) 21. (59) 


For each elementary solution (57) in V, we 
shall obtain the corresponding current densities 
from Eqs. (53) and (54). 

V.=Kr'*, 


2wy?/(w—nwy) —w |, 
ix —w | (60) 


Ja= —— —(rJar) 
n Or 


= —w J. 


5. DISCUSSION OF THE RESULTS 


A discussion of the practical results must be 
based on a numerical calculation of the two roots 
x, and x2 of Eq. (56). These roots, x; and Xx», 


must be calculated as functions of y=w/wy for 
different values of n. The solutions corresponding 
to +w and +m may also be obtained with +n 
and +w according to a remark on Eq. (41). 
Figure 3 shows a typical graph of x as a function 
of y as found for n=2, 3, ---. The case n=1 is 
exceptional. Figure 3 is self-explanatory and 
shows the behavior of the roots as functions of y. 
There are two regions, yry2 and ysys4, which give 
imaginary solutions x=x,+ix;. Curves repre- 
senting x and x,x; have been drawn separately. 
The limits of these regions are given by 


=0, 


¢=—n+(n?—1)!. 


Inside these regions of imaginary solutions we 
are interested in defining the point ys for which 
x,=1: 


x,= —B’/2A’=1/(2+¢/n) =1, 


Other points of interest are yo, where A’=0 
(¢{=—2n) and the points (y=n or n+v2) 
where the x roots are 1 and 0. All these solutions 
are easily obtained by the use of Fig. 4 relating 
¢ to y—n. The case n=1 is exceptional since the 
regions yiv2 and y3y, both disappear and the x 
roots are always real, one of them being 1 for 
any y value (Fig. 5). 

We now wish to discuss these results in relation 
to condition (59). The following cases are to be 
distinguished : 

I. Two real roots less than 1. None of these 
roots can be used for our problem. This occurs 
for 


yo<y<n—Vv2, n<y<ys, ys<y<n+v2. (62) 


II. Two real roots, one greater and the other 


4 


Fic. 5. Roots of Eq. 56 for n=1. 
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less than 1. The first root can be used. This 
occurs for 


y<yo, n—V20<y<n, n+v2<y. 


III. Two real roots greater than 1, both 
yielding possible solutions. This occurs for 


<1. 


IV. Two imaginary roots with x,<1. No 
useful root. This occurs for 


Ms. 


V. Two imaginary roots with x,>1, yielding 
two possible solutions. This occurs for 


Vi 


These results are summarized on Fig. 6, which 
shows the distribution of the different regions 
corresponding to cases I, II, III, IV, and V in 
the y,” plane. Curves have been drawn con- 
sidering as a continuous variable but the 
results have physical meaning only for integral 
values of n. In addition to the previously 
defined curves, the curve corresponding to 
x=n-—2 is also plotted. This is interesting in 
view of the approximation introduced in Section 
2, neglecting the additional field E’. (See 
Section 7.) ' 

In case II we obtain one possible solution (57) 
and there will be no difficulty in solving com- 
pletely our problem. We shall prove that the 
magnetron operated under such conditions 
behaves like a pure imaginary impedance and 
should be unable to sustain oscillations. 

In cases III and V we obtain two distinct 
solutions, yielding a potential 


Kor™. (63) 


This leaves an arbitrary constant, A,/Ke for 
instance, which we are unable to determine. 
It is thus impossible to draw any precise con- 
clusions in these cases. This difficulty seems to 
be connected with the simplifying assumption 
made in Section 4 [Eq. (45.A)] about the 
filament having a radius zero. A special discus- 
sion of a magnetron with finite filament radius 
should be made to answer the question. 

In cases I or IV no solution is obtained. Both 
roots are smaller than 1, which means that 
condition (59) cannot be fulfilled. If we try any 


one of the roots, we obtain an infinite oscillating 
electric field —d@V/dr on the filament (r=0). 
Instead of a magnetron working far from 
saturation and showing small oscillations around 
its steady conditions, we obtain large oscillations 
where the current emitted by the filament is 
zero (half of the time) or the saturation current 
for the other half periods. This shows that cases 
I and IV should represent the conditions under 
which a magnetron is able to oscillate and to 


yl 2 real roots<! / 

2real roots 
EJimag. roots<i 


imag. roots >i 


T 


Fic. 6. Summary of results. 


sustain high frequency oscillations in an outer 
circuit. The frequency is given by the corre- 
sponding y values: 


¥s<y<n—V2, 
n<y<n+v2. (64) 


We must now see how to build up the complete 
solution for the potential distribution inside a 
magnetron with 2m anodes, operated at the 
frequency w= ywy. According to the remarks at 
the beginning of Section 4 this complete solution 
will result from the superposition of elementary 
solutions corresponding to m and —n. Let us 
first suppose the corresponding roots x, and x_ 
both to be real and greater than 1 (case II). 


V.(r, 6, t) t—n@) + K_rz-eiot+ (65) 
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where the K, and K_ constants must be adjusted 
so as to reproduce, for r=6, the assumed po- 
tential distribution (2) on the anodes 


0, t) = Va(b)e! sin n8. (66) 


This gives 


1 1 
K,=——V,(b)b-*, K.=—V,.(b)b., 
2 
and 


1 
V.(r, 0, t) 


which gives the complete potential distribution 
inside the magnetron. 


6. INTERNAL RESISTANCE OF THE MAGNETRON 


Having obtained the distribution of potential 
in the magnetron, we are now in a position to 
compute the current densities and to calculate 
the current reaching the anodes. This might 
involve some difficulties when the anodes are 
narrow and far apart from each other, but if the 
magnetron is built with wide anodes and small 
intervals between them (a>), we may assume 
that each anode collects all the current flowing 
radially throughout an angle w/n. Thus the 
total current flowing to the first anode will be 


r=b, (68) 
0 


but there are m anodes (numbers 1, 3, ---, 2n—1) 
connected in parallel and the total current is n 
times greater 


ain 
Tor =en f (69) 
0 


J, will consist of two terms corresponding to 
+n, each of which is related to the corresponding 
voltage term by Eq. (60). Let us first consider 
the case of two real roots x, and x_ with a 
potential (67). The relation (60) has to be 
applied separately on the two terms of (67) 
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which gives 


2wy" 
| —x(r/b)** ( 


+x_(r/b)*- ( (70) 
wtnwy 


(71) 


Now 
f +2i/n. 
0 


Making r=), and using (69) to (71), we obtain 


2wx? 


+x-( -«)| (72) 
wt+tnwn 
And according to Eq. (2), the potential of the 


odd anodes is 
U,e'*'=eV,(b) 


piwt 
La 


8 sin (nB/2) 


me BKa. (73) 


In addition to radial currents, we have found 
that there are also angular components of the 
current density inside the magnetron. Thus we 
may have to add a correction term representing 
the current flowing to the edges of the anode; 
this would be something like 


en[ 0) a/n) jor, 


for r=b, (74) 


ér representing the apparent thickness of the 
electrode and the m factor taking care of the 
fact that there are m anodes connected in parallel. 
The edges of the first anode correspond to 6=0 
and x/n. Jas is known from Eq. (60) which 
must be applied separately to each elementary 
V, solution, yielding for each solution an equation 
of the form 


Jag = —— —(rJar) = 
n Or 


x? 
( -«) V.(r, 0, t). 
\w— 


Hence, for the superposition of the +2 solutions, 


71) 
ain 


/2) 


3) 
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Fic. 7. Static characteristics of a two-segment magnetron. 


as given by (67) 


| x4? 2wy? 
VG) 
r\ 
a 
n \w+nwu b 


Exponentials are 1 for @=0 and —1 for @=x/n 


so that the two terms in (74) are equal and 
one obtains 


+x Vor. (76) 
wtnwy 


x, and x_ are of the order of magnitude of n, so 
that the additional contribution is of the order 


or 
law, (76a) 


and will probably be small if the apparent 
thickness of the anodes is small. The important 
fact is that it comes out with the same phase 
angle as J,, when both roots x, and x_ are real 
(case II). In such cases the currents /,, and Jae 
retain pure imaginary amplitudes, while the 


potential keeps a real amplitude (73). This 
means that magnetrons operated under such 
conditions are equivalent to pure impedances 
without any (either positive or negative) resist- 
ance and seem to be unable to sustain oscillations 
in an outer circuit. 

It should be noticed that in the present paper 
the currents have been written with a sign con- 
vention opposite to the usual one in circuit 
theory. Hence the internal impedance Z, of the 
split anode magnetron should be defined as 


— 
+x( -«)| (77) 
w+n 


from Eqs. (72) and (73). This impedance is 
observed in the devices shown on Fig. 1 or 
Fig. 7. 


7. VALIDITY OF THE ASSUMPTION MADE IN 
SECTION 2, NEGLECTING THE 
ADDITIONAL FIELD E’ 


As noticed in Section 2, just after Eq. (18), 
the assumption has been made that the addi- 
tional electric field E’, induced by the alternating 
magnetic fields, could be neglected. It is now 
possible to compute the order of magnitude of 
this E’ field, which is determined by Eq. (20). 
Because the e*“'-"® factor appears in this field, 
as in all alternating quantities, Eq. (20) reads 


100 n* 
r or Or r 


where J, is the sum of terms (60) corresponding 
to each elementary term r* in V,, as may be 
seen from Eqs. (67) and (70). We may, for 
instance, write (70) this way 


Seah +n) =F n@) 
Sar | (79) 


r 2 
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And our Eq. (78) will read 
E,!=E,'(+n)+£,'(—n), 
100 n* 


—-—1r—rE,' (+n) ——E,’(+n) 
r or Or r 


(80) 


The solution is obvious since J,, depends upon 
r by rz. 
4ripowr? J, ar( +n) 
E,'(+n) 
r 


rJar( +n) 
(x(-+n)+2)?—n? 


as po=c and w/c=22/\. This shows that the 
additional field E’ is very small because of the 
factor (1/c)(r/X), so long as the dimensions of 
the magnetron are small compared with the 
wave-length, which corresponds to statement 
(19). The only case for trouble would be when 
one of the roots x should come near to +n—2, 
a condition to be discussed in a moment. 

The additional potential V’ on the anode r=b 
corresponding to E,’ can be computed. It should 
be reminded that r?J,, depends on r as r*tt, 
This gives 


V’=V'(+n)+V"(—n), 


(81) 


V'(+n)= -f E,’(+n)dr 
+2)[ (x4 +2)?—n?] 


This must be calculated for the middle of the 
first anode, i.e., for r=b, 9=2/2n. Formula (72) 
may be written 


Tor = €2b[ Jar( +0) + —n)], 


(82) 


(83) 
for r=b, @=2/2n, while 
J ar( +n) 
eV’= 
ch 
Jar(—n) 


This shows that the additional voltage V’ is out 
of phase with the current (7 factor) and very 


| r=b, 0=/2n. (84) 
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small except for the cases 


(85) 


x=+n—2. 


An order of magnitude can be obtained, if one 
remembers that the x values are of the order of n. 
This gives 


cd(n+2)(n+1) 


The addition of a voltage eV’ to the voltage 
¢V.(b) already assumed on the anodes is equiva- 
lent to adding a series impedance Z’= —eV’/J,,, 
the order of magnitude of which would thus be 


ab? 


eV’= 


(86) 


(87) 


Let us take, for instance, b>=75\ and A\=1 cm 
which will be a very unfavorable situation. 


0.32? 
(n+2)(n+1) 


which is of the order of magnitude of a fraction 
of an ohm, a very small quantity when compared 
to the thousands of ohms which electron tubes 
yield as internal resistance. 

Let us now discuss the conditions (85) for 
which the additional electric field takes a really 
important part. A glance at Fig. 3 shows that 
this condition can be fulfilled only for x(+n), 
between y=n—Vv2 and nm (where x, is slightly 
above 1) or above y=n+v2. The first case is 
for n=3, x=1, y=3—V2, 3, 3+v2; then n=4, 
etc. These points are plotted on Fig. 6. They 
all lie in the region where magnetrons have 
been proved to yield pure imaginary internal 
impedance and to be unfit for sustaining 
oscillations. 

It thus appears legitimate to neglect the role 
played by the additional field E’ induced by 
alternating internal magnetic fields. 


8. CONCLUSIONS 


The results of this study are summarized in 
Fig. 6 which shows the type of solutions obtained 
in different regions of the y,m plane. These 
results will perhaps be easier to understand for 
practical applications if we notice that Larmor’s 
angular velocity wy is proportional to the 
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magnetic field, which means 


wy =0.884 10°H, H in Gauss 
y=/wn =21310/AH. Ain cm (88) 


This enables a comparison with experimental 
results, which are usually given in terms of \/7 
corresponding to different types of oscillations. 

On the map drawn in Fig. 6, the different 
regions corresponding to the cases I to V defined 
in Section 6 [Eq. (62)] have been represented. 
The physical meaning of the various cases 
seems to be: 

I and IV. Large oscillations. 


ys<y<n—v2, n<y<n+vi. (64) 


II. No oscillations, the magnetron being 
equivalent to a pure imaginary impedance. 

III and V. The solution still contains an 
arbitrary constant; probably no oscillations. 
The different regions corresponding to these 
cases I to V are neatly divided for magnetrons 
corresponding to m=2, 3, 4 --- (number of 
anodes 2n=4, 6, 8 ---); one always finds two 
different regions, defined by conditions (64), 
where oscillations could possibly take place. 

For the usual split-anode magnetron (n=1; 
2 anodes) it is hard to foresee the conditions of 
oscillations, as the whole diagram changes just 
on the line »=1. These very typical circum- 
stances should make this magnetron very 
sensitive to all sorts of perturbations such as 
increase in the diameter of the filament, effect of 
large oscillations, etc. Experimental values on 
the usual split-anode magnetron are found in a 
paper by G. R. Kilgore,* with an attempt at 
theoretical explanation, which unfortunately 
does not take account of space charge effects, 


3G. R. Kilgore, Proc. I. R. E. (August, 1936). Reprinted 
by RCA in Radio at Ultra-High Frequencies (1940), p. 360. 


and therefore remains far away from actual 
conditions in magnetrons. Figure 7 is a repro- 
duction of Kilgore’s Fig. 3, showing the static 
characteristic for a split-anode magnetron. This 
means w=0, y=0. The magnetic field applied 
was 1.5 times the critical field. The average 
anode potential was 500 volts, and variations as 
large as +400 volts (giving E,4— Eg =800 volts) 
were applied. This means very large perturba- 
tions, for which our theory would certainly be 
only a rough first approximation. It appears 
from the curves that the negative resistance is 
zero for small AV and reaches a maximum of 
about —1500 ohms for AV=400 volts. Effi- 
ciencies of magnetrons operated under different 
conditions are given by Kilgore on page 372 and 
can be summarized as follows (H=1.5H,., H, 
critical field) : 


Efficiency 10% 2% 30% 40% 50% 
\H710- 5.45 6.7 8.1 9.8 12 (82) 
y=w/oy,. 0.39 0.32 0.265 0.218 0.178. 


If these y values are taken with a negative sign, 
they lie just between the horizontal m axis and 
the line y=n—v2. Another type of oscillation 
has been found on split-anode magnetrons with 
\H values around 12,000, which gives 


y~1.8. 


This value of y seems to correspond to the 
upper band of our diagram, as it lies between 1 
and 1+v2. 

A further prediction of the theory concerns” 
the possibility for a magnetron of type nm to 
work as magnetrons of types 3m, 5n, --- but 
with lower efficiency. 

Magnetrons n=3 (6 anodes) should behave 
differently, as the perturbation by the oscillating 
magnetic field is especially large in that case. 
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On the Decay of Slow Mesons 


M. D. bE Souza SANTOS 
Departamento de Physica, Universidade de Sao Paulo, Brasil 
June 20, 1942 


N order to measure accurately the mean life of slow 
mesons, an experiment was designed to measure their 
decay curve. The disintegration electron intensity was 
measured at different time intervals, after the mesons 
came to rest, through the use of a variable time delay 


coincidence set.! 
The counter arrangement used is illustrated in Fig. 1. 


Fic. 1. Counter arrangement. 


Counters (A, B, C) are connected to a threefold coincidence 
system with a coincidence time 7’ =2X10~* sec. The 
counters D are connected to a double channel mixer 
system? and are in anticoincidence with the vertical tele- 
scope (A, B, C). The counters E are connected to a double 
channel mixer system and the coincidences between the 
pulses from (A, B, C), delayed by a time interval T+7", 
and the pulses from E are measured for different values 
of T. T’ is a constant of the apparatus (coincidence time 
of the second mixer: 7’ = 1.2 X sec.). 

The spurious coincidences are due mainly to showers 
produced either in the lead block above the counters, 
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Fic. 2. Coincidences for different time delays. 


which absorbs the soft component, or to knock on showers 
produced in the iron block J. Both phenomena are elimi- 
nated by the anticoincidence counters D. Besides, since 
the pulses from (A, B, C) are delayed by a time interval 
T+T’, relatively to those from counters E, the chance 
coincidences can be shown to be about 10~ of the observed 
frequency (the true coincident phenomena such as showers 
being eliminated by the time delay system). 

When 7ST’, the observed frequency corresponds both 
to the decay of mesons which occurred in J and to the 
scattered mesons, but for 7>7”’ the observed frequency 
corresponds to the number of mesons disintegrated be- 
tween T—T’ and T7+7”. 

In an experiment such as this one, the results obtained 
can only be accepted provided that the breakdown time‘ 
of the counters used is smaller than the coincidence time 
used; that characteristic can only be obtained through the 
use of a multivibrator circuit’ and fast alcohol counters 
with a discharge started by a Penning effect. In this way, 
the time lags observed are smaller than the coincidence 
time used. The obtained efficiency of both apparatus and 
counters was nearly 100 percent and a complete account 
of the experimental set used together with the apparatus 
used to measure small time intervals will be published 
elsewhere. 

The preliminary results obtained are illustrated in Fig. 2 
and it is seen that the shape of the curve resembles that 
of a Gauss distribution around 7)>=5 X 10~ sec. the curve 
being symmetrical about the axis 4=5X10~® sec. which 
was considered as a measure of the “‘life’’ of the meson 
at rest. 

The important result to be emphasized here is that the 
shape of the curve is altogether different from what one 
would expect from an analogy with the laws of disintegra- 
tion of radioactive substances. If the time distribution of 
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the decayed mesons would be the same as observed for, 
say, beta-rays, the shape of the decay curve would be an 
exponential one.* We can therefore conclude that the laws 
responsible for both beta-ray phenomena are altogether 
different, and that the law of disintegration of mesons at 
rest is different from that of beta-rays. 

A complete account of the experiment together with 
new results that are being obtained will be published 
elsewhere. 

It is a pleasure to thank the British Council for the 
award of a Research Scholarship which has made possible 
the construction of the apparatus used and to Professors 
G. Wataghin and G. Occhialini for many helpful discus- 
sions. 

1 Prior to the author's measurements, the only coincidence experi- 
ments on the detection of delayed electrons were those of Montgomery, 
Ramsey, Cowie, and Montgomery, with negative results, and the experi- 
ments of Rasetti, with a similar geometrical arrangement, but with an 
altogether different experimental technique. 

2M. D. de Souza Santos, An. Acad. Bras. Sci. (in press). 

3G. Occhialini and M. D. de Sousa Santos, Symposium on Cosmic 
Rays in Rio de Janeiro (1941) to be published shortly. 

4 That is, the time interval between the moment in which an ionizing 


particle passes through the counter and the first electron avalanche 


reaches the wire. 
5M. D. de Souza Soman, An. Acad. Bras. Sci. 12, 183 (1940). 
6A. Moraes, An. Acad. Bras. Sci. (in the press). 


Temperature Diffuse Diffraction Maxima 
and Crystal Perfection 


DonaLp S. TEAGUE 
Department of Physics, Stanford University, Stanford University, 
alifornia 
July 31, 1942 

HE non-Laue x-ray diffraction maxima called by 
Raman “‘modified,”” and by Zachariasen and others 
“diffuse’’ have been observed by many investigators, 
using a continually growing list of crystals. Existing dif- 
ferences of opinion as to the detailed nature of these effects 
can be reconciled only on the basis of additional and 
quantitative observations. For these purposes no observa- 
tion is more critical than the observation of the angle of 
deviation of the diffraction maximum. The mere existence 
of a maximum of a monochromatic line off of coincidence 
with the Laue spot of the continuum proves nothing. 
Measurement of the intensity or breadth of the line maxi- 
mum may assist in classifying the maximum as something 
other than a Laue or Bragg reflection, though the test is 
rarely decisive because the theory of the width and in- 
tensity of the latter in an imperfect crystal specimen is 
impossibly complex, and because the maximum is often a 
mixture of ordinary and temperature diffuse effects. On 
the other hand, an angle of deviation varying with the 
angle of incidence is rather decisively indicative of a non- 
Laue maximum. Demonstration of a deviation which is 
constant with angle of incidence removes a crystal from 
consideration in any discussion of temperature diffuse 
reflection, since it indicates a structure so imperfect that 
the diffuse reflection, if any, is swamped by the ordinary 

reflection from misaligned blocks. 
A study of six crystalline substances has been conducted 
with a photographic spectrograph in an attempt to deter- 
mine the variation of the angle of deviation of Mo Ka 
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with angle of incidence. More than 200 exposures were 
made with the crystal faces inclined to a well-collimated 
incident beam at angles near the Bragg angle. The results 
are given in Table I. The order in which the crystals are 

TABLE I. Variation of angle of deviation of Mo Ka. A =reflecting 
planes used. B =range of angles of incidence over which line reflections 
(Mo Ka or a2) were observed. C=A deviation angle/A ae me 


angle observed. D =A deviation angle/A incidence angle, theoreti 
E =character of reflected lines. 


Crystal A B Cc D E 
Iron 100(?) largest broad, 
unsharp 
Zinc blende 110 1 »>4° Appr. zero — fairly 
sharp 
Barite 100 54’ Kai 0.00469 0.022 
Ka: 0.0203 024 Sharp 
Calcite 100 17’ Ka 0.26 .0295 cher 
Kaz 0.218 028 P 
Rochelle salt ? 6’ Ka 1.95(approx.)  .0266 sharp 


Kaz 2.4( 0264 
Quartz 110 negligible (approx) 


*S. Siegel and W. H. Zachariasen, Phys, Rev. 57, 795 (1940). 


listed is the order of increasing structural perfection.’ 
Apparently neither very perfect nor very imperfect crystals 
are favorable to studies of the angular shift of the reflected 
maximum, though the reasons differ in the two cases. 
With imperfect crystals the situation is as stated in the 
first paragraph; with a crystal as perfect as the quartz 
sample the reflection had an observable intensity only at 
or extremely near the Bragg angle of incidence and no 
systematic trend could be established even though five- 
hour exposures were used in an attempt to extend the 
range of the observations. With the crystals of intermediate 
perfection there is a qualitative confirmation of theory as 
shown in columns C and D but a satisfactory quantitative 
agreement occurs only in the case of barite. The tendency 
of the observed shift to exceed the theoretical, observed 
with calcite, has been noted before.? One might be inclined 
to regard it as an instrumental effect except that it does 
not appear with all crystals. — 


1R. M. Bozorth and F. E. Haworth, Phys. Rev. 45 821 (1934). 
Additional crystals have been arranged from evidence obtained in 
connection with the research herein described. 

2? P. Kirkpatrick Phys. Rev. 59, 452 (1941). 


The Spectrum of Chlorine I 


J. B. GREEN AND J. T. LYNN 
Mendenhall Laboratory of Physics, Ohio State University, Columbus, Ohio 
July 14, 1942 


EASUREMENTS on the Zeeman effect in the spec- 

trum of Cl I in the region Ad 3000-9200 have resulted 

in the identification of several new levels and some lines 

not listed by Kiess' as well as classification of many lines 
listed in his tables. 

The g values and g sums for a few configurations have 
been determined and show sharing of g sums for the 
3s*3p'4p and 3s*3p'5p. The discrepancies in the g sums can 
be practically completely accounted for on the basis of the 
overlapping of the 4p('D) and 5p(*P) parts of these con- 
figurations. 

A complete report will be published later. 


1 Kiess, Bur. Stand. J. Research 10, 827 (1933). 
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180 LETTERS TO 


Meson Theory and the Magnetic Moments 
of Protons and Neutrons 
H. FROHLICH 


Hi. H. Wills Physical Laboratory, University of Bristol, England 
July 27, 1942 


N various recent papers in this Journal! it was shown 

that an analysis of cosmic-ray data leads to the sugges- 
tion that the spin of cosmic-ray mesons is zero, and that 
they are described by a pseudoscalar field. The question 
now arises whether all mesons follow such a pseudoscalar 
theory, or whether there are also mesons with spin 1h 
(described by a vector field) which would have escaped 
detection in cosmic rays under the assumption that their 
lifetime is much shorter than that of pseudoscalar mesons. 
Such a combined theory was in fact proposed by Mller 
and Rosenfeld? in order to escape theoretical difficulties 
which are met in a pure pseudoscalar or vector theory. 
In this letter I should like to forward an argument in favor 
of the Mgller-Rosenfeld theory. Some time ago it was 
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shown that on the basis of the vector theory a qualitative 
understanding of the anomalous magnetic moments of 
protons and of neutrons could be reached.’ I have now 
carried out an analogous calculation for the pseudoscalar 
theory‘ and found that there are no contributions to the 
anomalous magnetic moments. A pure pseudoscalar meson 
theory, therefore, would not allow an understanding of the 
anomalous magnetic moments, whereas the Mgller-Rosen- 
feld theory does so and at the same time remains in agree- 
ment with cosmic-ray experiments. It may be hoped that 
further experiments, especially a measurement of the de- 
pendence of proper lifetime on height, would yield the 
required existence of two lifetimes. 

1R. F. Christy and S. Kusaka, Phys. Rev. 59, 405 and 414 (1941); 
J. R. Oppenheimer, Phys. Rev. 59, 462 (1941); J. R. Oppenheimer and 
J. Schwinger, Phys. Rev. 60, 150 (1941). 

2 Moller and Rosenfeld, Math. Phys. Comm. Acad. Sci. Copenhagen, 
17, 8 (1940); Rozental, Phys. Rev. 60, 612 (1941). 

3 Frdhlich, Heitler, and Kemmer, Proc. Roy. Soc. 166, 154 (1938). 

4 In the course of the calculations it turned out to be wrong to neglect 
in the matrix elements the recoil of the proton (neutron) as was done 


peeey. In the vector theory this changes only the numerical 
actor, but it is of great importance in the pesudoscalar theory. 
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